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ABSTRACT 


The  Final  Report  of  ONR  Grant  N00014-98-1-0617  -  “Design  of  a  Switched 
Reluctance  Motor  Drive  for  Electric  Propulsion”  -  has  been  divided  into  two  parts.  Part 
1  is  entitled  “Aspects  of  Switched  Reluctance  Motor  Drive  Application  for  Electric 
Vehicle  Propulsion”.  Part  2  is  entitled  “Operation,  System  Design  and  Control  of 
Switched  Reluctance  Motor  with  Multiphase  Excitation”. 

Electric  motor  drives  applied  to  propulsion  of  Electric  Vehicle  (EV)  have  to  typically 
traverse  a  wide  range  of  speeds  and  power  levels  during  their  operation.  It  is  then 
essential  to  determine  and  verify  the  adequacy  of  the  motor  design  during  each  assumed 
operating  condition.  Also,  it  is  necessary  to  determine  and  state  the  required  SRM 
continuous  power  rating  for  the  intended  EV  propulsion  application  before  starting  the 
design.  Towards  this  end,  the  Part  1  of  the  Final  Report  has  developed  a  new  sizing 
procedure  to  determine  the  SRM  peak  and  continuous  power  rating.  New  equations  were 
developed  to  estimate  the  SRM  losses  at  any  operating  condition  based  on  assumed 
losses  at  one  condition. 

Recent  literature  has  identified  that  multi-phase  operation  of  SRM  can  lead  to 
advantages  in  terms  of  torque  density,  efficiency,  torque  ripple  and  acoustic  noise.  The 
Part  1  of  the  Final  Report  introduces  new  coefficients  that  were  derived  with  the  help  of 
ANSYS™  based  FEA  studies  and  take  the  operation  of  each  multi-phase  configuration 
into  account.  These  coefficients  are  useful  when  designing  SRM  under  multi-phase 
operation.  The  design  ideas  are  then  verified  by  designing,  building  and  testing  the 
prototype  of  a  1-hp  four-phase  8/6  SRM  for  EV  propulsion. 


Torque  feedback  improves  the  performance  of  the  SRM  drive  and  aids  in  aspects 
such  as  torque  ripple  minimization.  Also,  position  sensorless  SRM  drives  are  currently 
being  widely  investigated  due  to  their  increased  reliability.  Under  multi-phase  operation, 
it  is  essential  to  account  for  the  mutual  flux  interactions  between  phases  in  order  to  get 
high  accuracy  in  these  tasks.  The  Part  1  of  the  Final  Report  proposes  a  new  Feed- 
Forward  Artificial  Neural  Network  (ANN)  based  system  to  provide  mapping  between  the 
SRM  terminal  variables  and  the  SRM  mutual  interaction  function  and  torque. 

Building  upon  the  previous  modeling  research  work,  the  Part  1  of  this  Final  Report 
proposes  a  simple  ANN  based  on-line  torque  estimator  that  eliminates  the  need  for  a 
torque  transducer,  adapts  to  the  characteristics  of  the  individual  SRM  and  can  be  easily 
extended  to  the  multi-phase  case. 

Finally,  the  Part  1  of  the  Final  Report  concludes  by  considering  an  implementation  of 
the  ANN  on  the  Texas  Instruments  expressDSP™  TMS320C6701.  The  pipelining 
features  of  the  TMS320C6701  and  the  Code  Composer  Studio  features  are  used  to  obtain 
low  execution  times  for  the  ANN.  The  implementation  verifies  the  feasibility  of  the 
proposed  ANN-based  methods. 


Chapter  One 
INTRODUCTION 


1.1  General  Overview 

The  main  advantages  of  the  Switched  Reluctance  Motor  (SRM)  in  terms  of  a  simple 
construction  (hence  higher  reliability),  high  torque  densities  and  high  fault  tolerance  have  been 
recognized  for  a  long  time.  In  addition,  the  required  power  electronics  converters  may  have 
simple  unipolar  configurations  that  are  simpler  to  design  and  operate  [1].  SRM  drives  have  been 
gaining  popularity  ever  since  the  advances  in  power  electronics  technologies  have  made  it 
feasible  to  utilize  their  advantages.  These  drives  are  now  considered  as  having  great  potential  for 
application  as  variable  speed  drives. 

Section  1.2  presents  the  basic  operating  principle  of  a  SRM  and  introduces  the  related 
terminology.  One  of  the  potential  applications  for  SRM  being  widely  investigated  is  for  the 
propulsion  of  Electric  Vehicles  (EV).  Towards  this  end.  Section  1.3  will  include  a  discussion  on 
Electric  Vehicle  (EV)  propulsion  applications.  Finally,  Section  1.4  will  identify  some  of  the 
current  applications,  design  and  operation  issues  that  need  better  understanding,  and  highlight  the 
research  efforts  in  these  areas. 

1.2  Basic  Operating  Principle  of  the  SRM  and  Associated  Terminology 

The  SRM  like  any  other  electric  motor  has  a  stator  (the  stationary  part)  and  a  rotor  (the 
rotating  part)  that  are  illustrated  in  Figure  1.1.  Both  the  stator  and  the  rotor  have  a  salient 
structure.  The  windings  are  present  only  in  the  stator  and  their  basic  unit  is  a  coil.  Each  coil  has 
two  sides.  There  can  be  several  winding  arrangements  determined  by  the  pitch  of  the  coil  sides 
forming  a  coil.  The  most  conunonly  used  arrangement  is  the  short-pitch  winding  in  which  both 
coil  sides  are  present  over  the  same  pole  (that  is,  a  coil  is  wound  around  a  pole). 
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Figure  1.1  SRM  cross-section  and  phase  winding. 

In  this  arrangement,  coils  on  physically  opposing  poles  are  connected  together  to  form  a 
phase  such  that  their  magneto-motive  force  (mmf)  aid  in  setting  up  the  flux  throughout  these  two 
poles.  The  other  winding  arrangements  that  have  been  reported  are  the  Jull  and  fractional  pitches 
[2-  3].  Due  to  the  structure  of  the  SRM  cross-section,  the  rotor  movement  causes  a  change  in  the 
reluctance  of  the  magnetic  circuit,  and  hence,  a  change  in  the  inductance  seen  by  the  stator  phase. 
To  understand  this  change  better,  we  now  consider  two  extreme  positions  of  the  rotor  with  respect 
to  a  stator  phase.  Figure  1.2  shows  the  first  position,  which  is  the  unaligned  position  in  which  two 
rotor  poles  have  equal  angular  displacement  with  respect  to  the  center  of  a  stator  phase;  in 
particular,  the  one  indicated  in  Figure  1.1.  In  the  ideal  case  of  a  perfectly  balanced  construction, 
exactly  equal  and  opposite  tangential  force  will  act  on  the  two  rotor  poles  when  this  stator  phase 
is  excited.  As  a  result,  the  net  shaft  torque  is  zero.  A  small  perturbation  in  any  of  the  two  possible 
directions  results  in  increased  torque  in  that  direction  and  one  of  the  two  rotor  poles  will  move 
towards  the  stator  pole  of  the  excited  phase  until  it  is  aligned  with  it  as  shown  in  Figure  1.3. 
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A  Finite-Element  Analysis  (FEA)  derived  plot  of  the  flux  distribution  shows  that  it  is 
symmetric  about  the  center  of  the  excited  stator  phase  (see  Figure  1.2).  Figure  1.3  shows  the  other 
extreme  position,  which  is  the  aligned  position  in  which  a  rotor  pole  is  aligned  with  a  stator  pole. 
When  the  stator  phase  is  excited  in  this  position,  the  tangential  force  acting  on  the  rotor  poles  is 
zero.  So,  the  shaft  torque  is  again  zero. 


Figure  1.2  SRM  unaligned  position. 


Figure  1.3  SRM  aligned  position. 
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A  FEA  plot  of  the  flux  distribution  in  this  case  shows  that  it  is  again  symmetric  about  the 
center  of  the  excited  stator  phase  and  the  number  of  flux  lines  are  much  higher  in  value  than  those 
at  the  unaligned  position.  Hence,  the  flux  linkage  (the  product  of  the  flux  and  the  number  of 
circuit  turns  linked  by  it)  is  much  higher  in  the  latter  case.  These  two  positions  are  equilibrium 
points  for  the  rotor  with  respect  to  an  excited  stator  phase.  The  SRM  thus  uses  the  simple 
principle  that  forces  in  any  magnetic  circuit  forces  act  to  align  the  system  in  the  position  of 
minimum  reluctance.  The  corresponding  forces  are  due  to  the  change  in  reluctance  of  the 
magnetic  circuit  and  the  associated  torque  acting  on  the  poles  is  known  as  reluctance  torque. 

The  flux  linkage  of  a  phase  at  a  constant  phase  current  increases  from  the  minimum  value  at 
the  unaligned  position  to  the  maximum  value  at  the  aligned  position.  At  a  given  rotor  position,  the 
flux-linkage  increases  with  increasing  phase  current.  Figure  1.4(a)  shows  the  plot  of  the  flux- 
linkage  characteristics  as  function  of  the  rotor  position  for  a  phase  of  the  four-phase,  8/6  (8  stator 
and  6  rotor  poles)  SRM,  considered  in  the  previous  illustrations.  In  these  cases,  the  unaligned 
position  is  at  -30°  and  the  aligned  position  is  at  0°  with  respect  to  the  excited  phase.  The  plot 
shows  that  the  rate  of  increase  of  the  flux  linkage  with  increasing  phase  current  decreases  as  the 
rotor  approaches  the  aligned  position  (0°).  This  is  due  to  the  saturation  of  the  magnetic  material 
used  to  build  the  SRM.  This  saturation  is  essential  for  maximizing  the  efficiency  of  the  energy- 
conversion  process  in  SRM  [1].  However,  it  leads  to  a  non-linear  system  that  is  difficult  to 
control. 

Let  Ns  denote  the  number  of  stator  poles  and  Nr  the  number  of  rotor  poles;  then,  the  angular 
displacement  of  the  rotor  pole  center  with  respect  to  the  stator  pole  center  at  the  unaligned 

position  is  e  =  °.The  angular  displacement  at  the  aligned  position  is  always  0  =  0°.  The 

flux-linkage  characteristics  such  as  the  one  in  Figure  1.4(a)  form  the  basis  for  all  the  simulation 
models  and  computations  that  are  done  for  the  SRM. 
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An  important  characteristic  that  can  be  derived  from  the  flux-linkage  data  is  the  torque 
characteristic  that  gives  the  variation  of  the  developed  electromagnetic  torque  with  rotor  angular 
position  at  different  phase  current  levels.  This  torque  characteristic  is  derived  using  the  result  that 
the  torque  in  a  rotary  electromagnetic  energy  conversion  system  is  given  by: 

where  Wc  is  the  co-energy  of  the  electromagnetic  system  and  is  given  by  the  area  under  the  flux 
linkage  characteristic  at  a  considered  rotor  position  till  the  considered  phase  current.  By  using  this 
calculation  for  the  flux-linkage  characteristics  of  the  SRM  in  Figure  1.4(a),  we  get  the  torque 
characteristics  shown  in  Figure  1.4(b)  for  the  case  when  only  one  phase  of  the  SRM  is  carrying 
current.  The  torque  characteristics  show  the  trend  that  was  described  before;  i.e.,  the  developed 
torque  is  zero  at  the  aligned  (0°)  and  unaligned  (-30°)  positions  for  any  value  of  the  phase  current. 
Thus,  the  torque  at  the  shaft  of  the  SRM  would  have  a  large  difference  between  the  maximum 
and  minimum  values  when  a  single  phase  is  only  excited;  this  difference  is  known  as  torque 
ripple.  This  would  lead  to  shaft  oscillations  at  low  speeds  and  increased  acoustic  noise  [4].  In 
order  to  minimize  this  torque  ripple  and  also  fully  utilize  the  SRM  power  output  capability, 
overlap  in  the  phase  conduction  periods  of  a  SRM  (and  hence  the  phase  current  waveforms)  is 
essential  when  it  is  applied  in  a  practical  motor  drive  application. 

1.3  EV  Propulsion  Application 

An  electric  motor  drive  for  EV  propulsion  has  to  satisfy  several  operating  conditions 
determined  by  the  performance  requirements  listed  as  part  of  a  specification  for  the  EV.  In  this 
application,  the  motor  drive  has  an  almost  continuously  varying  duty  in  terms  of  the  load  torque 
and  speed  requirements.  The  EV  America  Technical  Specifications  [5],  the  ANSI-IEEE  Standard  for 
Rotating  Electrical  Machineiy  for  Rail  and  Road  Vehicles  [6],  and  the  NEMA  Standard  for  Motors  and 
Generators  [7]  address  the  rating  of  electric  motors  used  in  EVs. 


6 


These  standards  require  that  the  peak  and  continuous  power  ratings  be  identified  separately  for  a 
complete  specification  of  a  motor  used  for  electric  propulsion  of  an  EV.  For  a  particular  electric 
propulsion  application,  the  EV  main  specifications  are  vehicle  maximum  mass  (niev),  maximum 
and  rated  speeds  (Vev,niax  and  Vev,r).  acceleration  and  deceleration  times  (ta  and  tj),  maximum  slope 
or  gradient  (a),  rolling  coefficient  (kr),  EV  frontal  area  (Af),  dragging  coefficient  (kd),  headwind 
speed  (Vhw)«  wheel  radius  (r*)  and  gear  box  ratio(s).  In  addition,  the  considered  performance 
specifications  must  also  include  the  driving  schedules  (city  and  highway)  and  other  particular 
operating  conditions.  The  specifications  determine  the  road-load  characteristics  (see  Fev  Figure 
1.5)  that  the  EV  has  to  overcome,  calculated  using  the  following  equation  [8]  (ad  -  density  of  air): 

FevM  =  K^eyg  +  +  m,,gsin(a)  (1.1) 

Figure  1.5  also  shows  the  envelope  Fem  of  the  EV  motor  drive  characteristic  that  depends  on 
the  motor  drive  peak  power  rating,  rated  speed  and  maximum  speed. 


Figure  1 .5  The  road-load  and  motor-drive  forces  as  function  of  vehicle  velocity. 
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The  motor  drive  has  to  traverse  its  entire  operating  characteristic  (that  is,  constant-torque  and 
constant-power  regions).  The  envelope  determines  the  performance  obtained  from  the  EV.  The 
shaded  area  in  Figure  1.5  represents  the  integral  that  determines  the  acceleration  performance  to 
rated  speed  that  can  be  obtained  from  the  EV.  The  road-load  characteristic  for  a  given  EV  is  not 
constant,  but  changes  with  the  considered  condition.  When  designing  the  motor  drive 
characteristics  for  EV  propulsion,  one  has  to  meet  the  considered  performance  specifications, 
while  achieving  a  high  average  drive  efficiency  and  low  drive  weight  and  cost  [9].  The  vehicle 
speed  is  related  to  the  electric  motor  speed  through  the  gear  ratio  and  wheel  radius.  In  the  case  of 
EV,  there  is  usually  a  single  gear  ratio,  and  the  desired  characteristic  at  the  drive  shaft  is  obtained 
through  electronic  control  (i.e.,  electronic  gearing).  Hence,  we  can  carry  out  all  the  studies  in 
terms  of  the  vehicle  speed,  thereby  eliminating  the  need  to  consider  the  values  of  the  gear  ratio 
and  wheel  radius. 

1.4  SRM  Application  Issues  and  Contributions  of  this  Report 

For  EV  applications,  it  is  essential  to  identify  and  state  both  the  peak  and  continuous  power 
ratings  of  the  SRM  [6].  Towards  this  end,  there  is  a  need  at  the  design  stage  to  formulate  a 
methodology  to  calculate  the  continuous  power  rating  of  the  SRM  by  predicting  the  efficiency 
contours  of  the  SRM  based  on  the  assumed  efficiency  at  one  operating  point.  Thus,  it  is  useful  to 
have  scaling  equations  relating  losses  at  one  (known)  operating  point  of  the  motor  to  those 
(unknown)  losses  at  any  other  operating  point.  It  is  required  to  assume  the  losses  at  one  operating 
point  because  the  SRM  is  yet  to  be  designed  and  built  at  this  stage.  This  would  allow  the  designer 
to  estimate  the  motor  losses  under  each  operating  condition  in  order  to  check  the  adequacy  of  the 
motor  thermal  design  for  each  operating  condition.  An  estimation  of  the  thermal  duty  would  also 
allow  proper  choice  of  the  safety  margins  used  in  the  design  which  would  then  determine  the  life 
of  the  SRM  drive  depending  on  the  temperature  attained  by  the  winding  under  different  operating 
conditions  [6-7]. 
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Chapter  Two  first  presents  a  detailed  SRM  sizing  procedure  [12]  for  the  propulsion  of  EV 
along  with  the  development  of  the  equations  that  allow  relating  losses  at  one  operating  point  to 
another  one.  Using  this  approach,  a  code  was  developed  in  C++  for  sizing  the  SRM. 

In  recent  literature,  the  potential  for  multi-phase  operation  (more  than  one  phase  conducting 
current  simultaneously)  has  been  identified  to  obtain  higher  torque  densities  and  higher 
efficiencies  [2-3, 10-11].  These  literatures  only  provide  the  means  for  understanding  the  operation 
under  these  conditions.  There  is  however  a  need  for  a  methodology  to  design  the  SRM  under 
multi-phase  operating  conditions,  while  taking  into  account  the  loading  constraints  of  the 
magnetic  material  used  to  build  the  SRM  and  the  desired  output  torque.  Chapter  Three  includes 
an  approach  for  incorporating  the  multi-phase  operation  into  the  SRM  design  equations  [13].  The 
theory  for  this  approach  was  developed  jointly  by  Mr.  Edilberto  Hall  Mitre  and  myself.  I 
suggested  the  idea  of  introducing  coefficients  to  account  for  the  multi-phase  operation  into  the 
already  existing  design  equations  [1]  for  SRM  under  single-phase  operation.  Mr.  Hall  Mitre  and  I 
determined  the  values  of  these  coefficients  based  on  the  results  of  the  FEA  studies.  I  performed 
all  the  FEA  studies.  Mr.  Hall  Mitre  then  developed  the  MATLAB™  code  for  designing  the  SRM 
using  this  approach.  Using  this  sizing  and  design  approach  as  an  aid,  a  1  hp  four-phase  8/6  SRM 
prototype  was  designed  by  me  and  built  by  Baldor  Motors  and  Drives  (Fort  Fort  Smith,  AR).  Mr. 
Roberto  Marcelo  Schupbach  and  I  then  carried  out  static  measurements  on  the  SRM  to  quantify 
the  performance.  A  description  of  this  design  and  the  obtained  performance  are  included  in 
Chapter  Three. 

Since  the  SRM  is  designed  to  be  non-linear  and  has  a  coupled  magnetic  circuit,  it  is  expected 
that  superposition  principles  cannot  be  applied  to  predict  its  behavior  under  multi-phase 
operation.  However,  it  is  essential  to  identify  and  predict  the  non-linear  torque  characteristics  of 
the  SRM  under  multi-phase  operating  conditions  for  tasks  such  as  torque  feedback  or  torque 
ripple  minimization. 
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In  addition,  it  is  essential  to  identify  and  adapt  to  the  characteristics  of  each  individual  SRM 
due  to  variations  in  their  characteristics  due  to  manufacturing  tolerances  and  variance  in  the 
magnetic  material  characteristics.  For  these  tasks,  it  is  necessary  to  develop  feasible  solutions  that 
provide  adequate  accuracy.  To  this  end.  Chapter  Four  proposes  Feed-Forward  (FF)  Artificial 
Neural  Network  (ANN)  based  solutions  to  achieve  this  non-linear  mapping  between  the  SRM 
terminal  variables  and  its  performance  [14].  The  idea  of  using  ANN  for  this  task  was  initially 
proposed  by  me.  Mr.  Schupbach  and  I  then  tried  out  several  cases  of  using  ANN  for  the  torque¬ 
mapping  task.  The  simulation  studies  by  incorporating  ANN  into  the  SRM  model  were  done  by 
Mr.  Schupbach.  A  new  FF  ANN  based  torque  estimator  is  proposed  in  Chapter  Five  that  can  learn 
and  adapt  to  the  characteristics  of  the  individual  SRM  based  on  on-line  measurements  [15]. 

Finally,  Chapter  Six  considers  an  implementation  of  the  considered  FF  ANN  based  solutions 
on  the  TI  TMS320C6701  DSP  and  reports  on  the  performance  obtained.  A  scheme  is  also 
suggested  for  communicating  the  torque  feedback  fi’om  the  TMS320C6701  to  the  TMS320F240 
EVM  (the  SRM  PWM  generator)  [16]. 

Throughout  this  report,  references  are  given  at  the  end  of  each  chapter. 
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Chapter  Two 


SIZING  A  SWITCHED  RELUCTANCE  MOTOR 
FOR  ELECTRIC  VEHICLES 


2.1  Introduction 

The  EV  America  Technical  Specifications  [1],  the  ANSI-IEBE  Standard  for  Rotating  Electrical 
Machinery  for  Rail  and  Road  Vehicles  [2],  and  the  NEMA  Standard  for  Motors  and  Generators  [3]  address 
the  rating  of  electric  motors  used  in  EVs.  Ifence,  these  are  used  when  appropriate  throughout  this  chapter. 

These  standards  require  that  the  peak  and  continuous  power  ratings  be  identified  separately  for  a 
complete  specification  of  a  SRM  used  for  electric  propulsion  of  an  EV.  A  detailed  procedure  is  presented  in 
this  chapter  to  determine  these  power  ratings  given  the  vehicle  specifications.  To  determine  the  continuous 
power  rating  of  the  SRM,  it  is  necessary  to  know  the  variations  of  the  SRM  losses  with  vehicle  speed  and 
power  requirement.  To  this  end,  the  report  work  has  led  to  the  development  of  new  equations  that  are  used 
in  the  sizing  calculations  to  estimate  the  variations  of  the  SRM  losses  with  speed  and  power  output.  The 
equation  derivations  are  included  in  Appendix  A. 

For  a  particular  electric  propulsion  application,  the  EV  main  specifications  are  vehicle  maximum  mass 
(mev),  maximum  and  rated  speeds  (Vevjnax  and  Vevj),  acceleration  and  deceleration  times  (t*  and  td), 
maximum  slope  or  gradient  (a),  rolling  coefficient  (1q),  EV  frontal  area  (Af),  dragging  coefficient  (kd), 
headwind  speed  (Vhw),  wheel  radius  frw)  and  gear  box  ratio(s).  The  EV  specifications  considered  in  this 
study,  listed  in  Table  2.1,  are  the  minimum  specifications  for  an  EV  qualifying  as  an  EV  America 
“production-level”  vehicle  [1]. 


Table  2.1  The  EV  Specifications. 


Vehicle  mass 

1000  kg  (GM  Impact) 

Mass  considered  in 
calculations  {mj 

1 172  kg  (vehicle  with  2  166- 
lb  occupants) 

Rolling  Coefficient  (k) 

0.013 

Acceleration  duty 
(O 

0  to  50  mph  in  13.5  s 

Aerodynamic  Drag  Coefficient 

0.29 

Maximum  speed 

70  mph  (drive  should  be  able 
to  operate  continuously 
under  this  condition) 

Frontal  Area  (A,) 

2.13 

Gradient  (a) 

0%,  3%  at  55  mph,  6%  at  45 
mph  and  capable  of 
ascending  a  25%  gradient 

Wind  Velocity 

5  mph 
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The  force  vs.  speed  characteristic  (or  the  so-called  “road-load”  characteristic)  is  determined  from  the 
EV  main  specifications  (see  Figure  1.5).  The  considered  performance  specifications  also  include  the 
driving  schedules  (city  and  highway)  to  be  satisfied  and  other  operating  conditions.  The  driving  schedule 
data  are  given  in  the  form  of  vehicle  velocity  at  time  instants  equally  spaced  (see  Figure  2. 1).  The  sizing 
calculations  presented  here  consider  the  Federal  Test  Procedure  (FTP)  and  the  Highway  Fuel  Economy 
Driving  Schedule  (HWFET)  representing  city  and  highway  driving  conditions,  respectively  [4].  The  other 
considered  operating  conditions  are  based  on  the  minimum  performance  specifications  for  a  commercial 
EV  given  in  [1]. 

The  vehicle  speed  is  related  to  the  electric  motor  speed  through  the  gear  ratio  and  the  wheel  radius.  All 
the  results  are  reported  in  terms  of  the  vehicle  speed,  thereby  eliminating  the  need  to  consider  the  values  of 
the  gear  ratio  and  the  wheel  radius  in  this  study. 

In  the  analysis,  the  SRM  total  losses  are  equal  to  the  sum  of  copper  and  total-iron  losses  (i.e.. 


Figure  2. 1  The  FTP  is  composed  of  the  UDDS  (Urban  Dynamometer  Driving  Schedule) 
followed  by  the  first  505  seconds  of  the  UDDS. 
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WsRM  =  Wcu  +  Wpe)-  Rotational  losses  are  not  considered  although  they  could  be  included  as 
suggested  in  [2].  In  the  following  sections,  we  explain  the  proposed  procedure  for  determining 
(sizing)  the  ratings  of  a  SRM  for  EV  propulsion.  This  explanation  is  followed  by  the  results  of 
applying  this  procedure  to  sizing  a  SRM  for  the  EV  specified  in  Table  2.1.  The  Chapter  also 
includes  a  discussion  on  the  continuous  power  ratings  and  the  thermal  duty  of  other  electric 
motors  as  compared  to  the  SRM  for  the  same  EV  propulsion  application. 

2.2  Procedure  for  Sizing  the  SRM 

For  a  complete  specification  of  the  electric  motor,  (i)  the  peak  and  continuous  power  ratings  should  be 
identified  separately,  and  ( ii )  the  peak  power  (or  shcat-time  capability)  should  be  corrected  according  to  the 
requirements  imposed  by  all  considered  driving  schedules  and  operating  conditions  [1-3].  In  [5],  the  peak 
power  rating  is  also  called  the  rated  or  continuous  power  of  the  electric  motor;  this  might  lead  to  over¬ 
dimensioning  the  electric  motor.  We  can  determine  the  peak  power  rating  and  the  continuous  power 

rating  _ by  taking  into  account  the  vehicle  specifications,  considering  the  principles  of  operation  of 

SRMs,  and  assuming  (or  knowing)  the  motor  loss  values  at  only  one  operating  point. 

Initially,  the  peak  power  rating  Pc^  of  the  SRM  (i.e.,  the  electric  motor)  is  calculated  from 
the  specified  acceleration  duty  as  shown  in  [5].  This  value  is  then  compared  with  (a)  Pemmax,  the 
power  for  cruising  at  the  vehicle  maximum  velocity  [5],  as  well  as  the  powers  Pemg  required  when 
the  EV  is  moving  at  (b)  55  mph  on  a  3%  gradient,  (c)  45  mph  on  a  6%  gradient,  and  (d) 
ascending  a  25%  gradient  at  low  speed.  The  value  of  Pemp  is  changed  to  Pemmax  or  the  highest  Pemg 
if  any  of  these  values  is  greater  than  the  initial  value  of  Pemp-  To  determine  the  continuous  power 
rating  Pemcont  of  the  SRM,  we  need  to  identify  that  particular  operating  condition  imposing  the 
most  severe  thermal  performance  requirement  on  the  SRM  (i.e.,  the  condition  with  the  highest 
losses). 

This  is  done  by  determining  an  equivalent  continuous  power  rating  Pcont.i  at  a  chosen  speed  (we  chose 
the  rated  speed  Vem,r)  for  each  operating  condition  i  (for  example,  like  the  driving  schedule  in  Figure  2.1). 
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Figure  2.2  illustrates  that  the  power  losses  change  as  function  of  the  vehicle  (motor)  velocity.  These  losses 
increase  the  tenperature  of  the  electric  motor.  By  definition,  the  losses  dissipated  at  Pconu  the  chosen 
speed  are  identical  to  those  dissipated  throughout  the  operating  condition  (where  the  speed  changes). 
Hence,  the  ten^jerature  rise  is  the  same.  This  equivalent  continuous  power  rating  is  calculated  at  the  same 
chosen  speed  for  all  conditions  to  identify  that  condition  resulting  in  the  highest  PcantM  which  then 
determines  the  SRM  continuous  power  rating  Pemcom  valid  only  at  the  chosen  speed,  hi  other  words,  the 
SRM  is  capable  of  delivering  this  Pemcom  at  the  chosen  speed.The  SRM  temperature  rise  under  this 
continuous  power  rating  should  be  within  the  limits  specified  in  [2]  for  each  class  of  insulation.  In  other 
words,  the  SRM  will  have  the  same  losses  at  the  continuous  power  rating  Pmcom  and  the  rotor  speed 
corresponding  to  Vem,r  as  the  average  losses  corresponding  to  the  most  severe  operating  condition. 

Due  to  the  variation  of  SRM  losses  with  vehicle  speed  (see  Figures  2.2  to  2.4),  a  certain  condition  (e.g., 
cruising  at  maximum  vehicle  speed)  may  produce  more  losses  than  when  the  electric  motor  is  delivering 
peak  power  at  the  chosen  speed.  To  simulate  this  most-severe  condition  and  test  thermally  the  electric 
motor,  the  motor  output  would  have  to  be  increased  above  P^p  resulting  in  Pmcom  being  greater  than  Pemp, 
at  the  chosen  speed.  However,  when  Pemcom  is  specified  at  another  speed  at  which  the  SRM  efficiency  is 
lower,  it  should  be  less  than  Pemp- 

By  using  the  sizing  procedure  described  in  this  report,  the  SRM  does  not  need  to  be  capable  of 
delivering  Pemp  continuously  at  all  the  speeds  in  the  operation  range  of  the  drive  as  in  [5],  but  it  has  to  be 
designed  to  dissipate  the  losses  occurring  during  all  the  specified  conditions  of  EV  operation. 

In  order  to  calculate  P^^.  it  is  necessary  to  know  the  SRM  losses  at  each  power  requirement  imposed 
by  the  operating  condition;  unfortunately,  this  is  not  possible  since  the  SRM  (or  the  electric  motor)  has  yet 
to  be  designed.  To  relate  assumed  (or  known)  motor  losses  at  one  particular  operating  point  (for  example, 
we  select  P^  and  rated  speed  v^,,)  to  the  equivalent  continuous  power  rating  P^  (also  at  rated  speed  v^,) , 
we  define  the  ratio  s,  at  the  SRM  rated  speed  between  P^  and  P^^. 

=  Pern/ Pemp-  (2-1) 
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For  each  specified  operating  condition  i  (or  driving  schedule),  is  calculated  by  solving  for  s.  in  an 
iterative  manner  using  the  following  motor-loss  equation  (refer  to  the  Appendix  A  for  its  derivation): 

{(1-c)  +  cxf}  +  5,  X c  X  (1-f)  =  (1-c)  xK,+  c  xfx  cx(l-f)x  (2.2) 

The  Newton’s  algorithm  is  used.  A  first  estimate  of  s^  is  obtained  ignoring  the  exponent  0.8.  The 
coefficient  c  represents  the  ratio  between  the  total-iron  loss  and  total  losses  of  the  electric  motor,  and  the 
coefficient  /  represents  the  ratio  between  the  eddy  current  and  total-iron  losses.  These  coefficients 
correspond  to  the  peak  output  power  at  any  arbitrary  speed  (we  selected  the  SRM  rated  speed);  that  is, 
c  =  and/=  W^/  We  used  typical  values  for  these  two  coefficients  based  on  SRM  designs 

presented  in  the  literature  [6-8].  The  accuracy  of  the  calculated  or  estimated  P^  depends  on  the 
“closeness”  between  the  assumed  (or  known)  loss  values  and  those  of  the  designed  and  built  SRM.  In 
(2.2),  the  copper-loss  coefficient  K,,  the  eddy-current-loss  coefficient  and  the  hysteresis-loss  coefficient 
corresponding  to  the  equivalent  continuous  power  rating  F„^,  are  given  by  (see  Appendix  A  for  their 


derivations): 

^  1  ^em,r  ^  Pdi 

P 

^  ^emp  i=l  ^evj 

(2.3) 

N  P  V 

(2.4) 

■j-  _  1  1  yr 

W  nO.8  0,2 

*  emp^em,r  »=1 

(2.5) 

With 

C,  =  PT^Z  for 

or 

C,  =  X  for  >  \n.r 

N  represents  the  number  of  time  instants  in  the  considered  driving  schedule  (or  operating  condition).  For 
example,  the  FTP  driving  schedule  has  1874  points. 

The  power  P^  in  the  above  equations  is  the  power  required  from  the  SRM  (electric  motor)  at  the  center 
of  each  time  interval  in  a  drivmg  schedule  assuming  a  linear  velocity  profile  between  any  two  consecutive 
time  instants.  This  power  is  given  by: 
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(2.6) 


The  required  force  components  are  determined  as  follows: 

1)  is  the  acceleration  component  to  accelerate  the  vehicle  from  velocity  at  time  instant  t,_, 
to  velocity  v^,at  time  instant  t..  This  force  is  given  by: 


2)  is  the  force  required  to  overcome  the  road  load  [5]  at  the  velocity  in  the  center  of  the  considered 

time  interval. 

A  negative  value  of  means  that  the  vehicle  requires  braking  power,  that  is,  the  electric  motOT 
participates  in  regenerative  braking.  In  this  case,  the  calculated  is  inverted  in  sign  for  the  purpose  of  only 
estimating  the  motor  losses.  Any  additional  force  required  during  brakmg  over  the  capability  of  the  electric 
motor  at  a  given  vehicle  speed  is  assumed  to  be  provided  by  a  mechanical  braking  system. 

The  constant  takes  into  account  the  rotating  masses  within  the  EV.  The  total  force  ( F„+  FJ  during 
acceleration  must  be  checked  if  it  is  less  than  or  equal  to  the  SRM  peak  capability  based  on  the  value  of  P^ 
calculated  earlier.  If  not,  the  peak  SRM  power  P^  must  be  made  equal  to  the  value  corresponding  to  this 
total  force.  The  C-h-  program  for  sizing  the  SRM  is  included  in  Appendix  B. 

2.3  Results  of  the  Sizing  Calculations 

Assuming  a  ratio  of  5:1  between  the  maximum  and  rated  rotor  speeds  [5],  the  sizing  procedure 
produces  the  following  results  for  the  EV  specifications  in  Table  2.1: 


A.  SRM  peak  ( output)  power  rating 

(1)  Accelerating  from  0  to  50  mph  in  13.5  sec  =  28  kW. 

(2)  Ascending  a  25%  gradient  (at  20  mph)  =  28  kW. 

(3)  Maintaining  55  mph  in  a  3%  gradient  =20kW. 

(4)  Maintaining  45  mph  in  a  6%  gradient  =  21  kW. 
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(5)  Cruising  at  the  maximum  vehicle  velocity 


=  19  kW  at  the  maximum  motor  speed. 


Thus,  the  peak  output  power  rating  of  the  SRM  is  specified  as  28  kW. 


B.  Continuous  (output)  power  rating 

The  values  of  K^,  K^,  and  s,  were  determined  using  the  driving  schedules  FTP  and  HWFET  as  well 
as  the  power  requirements  at  the  maximum  vehicle  velocity  of  70  mph  and  for  sustained  speeds  on 
gradients.  Values  of/=  0.333  and  c  =  0.4  (considered  typical  for  a  SRM)  were  assumed  [6-8].  The 
calculated  results  are  listed  in  Table  2.2. 

Using  (2.3),  (2.4),  (2.5)  and  finally  (2.2)  with  N=\,P^=  P^,  v^,=  ,  P^^  of  19  kW  at  the 

maximum  speed  corresponds  to  of  18  kW  at  the  SRM  rated  speed  with  the  K^,  and  s.  values 
listed  in  Table  2.2.  This  duty  of  cruising  at  the  maximum  vehicle  speed  results  in  the  highest  value  of  P^. 


Table  2.2  The  Continuous  Power  Ratings  at  Rated  Speed  for  the  SRM 


Condition 

K,(3) 

KJ4) 

KJ5) 

s,(2.2) 

FTP 

0.121 

0.381 

0.252 

0.172 

5 

HWFET 

0.103 

1.222 

0.455 

0.324 

9 

Maintaining  55  mph  in  a  3  %  gradient 

0.182 

2.806 

0.507 

0.601 

17 

Maintaining  45  mph  in  a  6  %  gradient 

0.233 

2.411 

0.558 

0.592 

17 

Cruising  at  maximum  vehicle  speed 

0.133 

3.333 

0.446 

0.626 

18 

The  continuous  power  rating  P^^  of  the  SRM  is  specified  as  18  kW  at  the  SRM  rated  speed  (see  Table 

2.2). 

Assuming  a  safety  factor  of  1.2,  the  output  power  ratings  of  the  SRM  to  be  designed  and  built  are: 

Peak  power  =  34  kW  at  and  above  the  rated  speed  of  the  SRM. 

Continuous  power  =  23  kW  at  the  maximum  speed  of  the  SRM. 

Figure  2.2  shows  the  ideal  (torque  or  force)  envelope  of  the  SRM  drive  (referred  to  the  EV  wheels),  the 
road-load  characteristic  (as  per  (1.1)  for  the  EV  specifications  in  Table  2.1)  and  the  SRM  loss  variation 
(when  the  SRM  is  developing  the  power  required  to  satisfy  the  road-load  characteristic).  Here,  Area  I  is  the 
specified  operating  region  (when  the  vehicle  speed  is  less  than  the  maximum  specified  value  of  70  mph) 
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and  Area  II  is  the  region  where  the  EV  is  still  capable  of  operating  due  to  the  SRM  peak  power  output 
capability  Pemp-  However,  the  SRM  drive  is  not  capable  of  operating  continuously  in  this  region  due  to 
thermal  considerations.  It  is  clear  from  the  figure  that  this  EV  is  not  capable  of  travelling  at  velocities 
greater  than  about  90  tt^ih. 


2.4  SRM  Loss  Trends 


Figure  2.2  The  electric  propulsion  drive  envelope,  EV  road-load  characteristics,  SRM  loss 

variations  and  the  operating  areas. 

hi  order  to  achieve  a  low  loss  dissipation  or  low  thermal  duty  for  the  electric  motor  and  also  have  the 
optimum  efficiency  for  the  EV  drive,  it  is  essential  to  understand  the  behavior  of  each  motor  loss 
component  as  a  frinction  of  the  EV  speed  and  power  requirement.  To  this  end.  Figure  2.3  shows  the 
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variations  of  the  SRM  losses  with  vehicle  velocity  when  the  SRM  is  developing  a  constant  torque  equal  to 
that  value  required  at  the  vehicle  maximum  velocity.  Figure  2.4  shows  the  variations  of  SRM  losses  with 
vehicle  velocity  when  the  SRM  is  developing  the  power  required  to  overcome  the  specified  road  load,  hi 
these  figures,  each  loss  conponent  is  normalized  (in  per-unit)  with  respect  to  the  total  loss  at  the  peak 
power  output  and  rated  speed.  The  total  losses  are  plotted  using  the  assumed  values  of  /and  c.  For  a 
constant  torque  demand,  the  SRM  losses  increase  with  increasing  vehicle  velocity  (and  drive  speed).  This 
is  mainly  due  to  the  increase  in  the  iron  losses  w  vehicle  speed;  in  particular,  the  eddy-current  loss 
corcqionent  exhibits  the  most  rapid  increase.  Finally,  it  is  clear  fi-om  Figure  2.4  that  the  losses  are  the 
highest  when  the  vehicle  is  cruising  at  the  maximum  vehicle  speed. 

Next,  we  analyzed  the  effects  of  the  contributions  of  the  loss  components  on  the  continuous  power 
rating  of  the  SRM.  For  c  (=  vaiying  between  0.3  to  0.7  with/=  0.333,  the  continuous  power 

rating  at  rated  speed  varies  between  0.5  to  1 .03  times  the  peak  power  rating  P^.  On  the  other  hand,  varying 
fi=W^/  )  between  0.3  to  0.7  with  c  =  0.4  results  into  the  continuous  power  rating  at  rated  speed 
vaiying  between  0.6  to  1.07  times  of  P^.  Since  the  iron  loss  shows  an  increase  with  the  drive  speed,  it  is 
desirable  to  keep  the  iron  losses  at  a  lower  fiiaction  of  the  total  losses  by  choosing  c  low  and  to  design  the 
SRM  with  low  eddy-current  losses  by  choosing/low.  This  reduces  the  losses  at  the  maximum  cruising 
velocity,  and  hence,  it  reduces  the  loss-dissipation  requirements  for  the  SRM.  Thus,  the  continuous  power 
rating  and  the  loss-dissipation  duty  imposed  on  the  SRM  are  also  fimctions  of  the  contributions  of  the 
individual  loss  components  to  the  total  motor  losses;  these  are  determined  to  a  certain  degree  during  the 
design  stage.  The  development  of  a  design  methodology  using  these  ideas  is  not  within  the  scope  of  this 
report 
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(nd)  iio’i 


Figure  2.3  Variation  of  loss  components  and  total  loss  with  vehicle  velocity  at  constant  torque. 


Figure  2.4  Variation  of  loss  components  and  total  loss  with  vehicle  velocity. 


2^  Comparison  with  Other  EV  Electric  Motors 

To  get  a  better  insight  into  the  thermal  duty  of  the  SRM  as  compared  to  other  electric  motors  used  for 
this  EV  application,  a  similar  sizing  study  for  other  popular  drives  was  carried  out.  The  considered  electric 
motors  are  the  Permanent  Magnet  Synchronous  Motor  (PMSM),  the  Induction  Motor  (IM)  and  the  DC 
Motor.  The  performance  of  each  motor  is  described  below  and  the  sizing  results  are  presented  in  Table  2.3. 
The  results  are  presented  for  a  drive  with  =  2  and  for  a  drive  with  =  5  (long  constant- 

power  region).  See  Appendix  B  for  the  MATLAB™  code  for  sizing  these  motors. 


Table  2.3  Results  of  Sizing  Calculations  for  Different  Motor  Types. 


Ratio 

V  /v 

«n/nax  tmf 

2 

P^=  59  kW 

5 

Pe^  =  34  kW 

s 

P 

Bl/ilW 

S 

P 

SRM 

0.36 

0.34 

20 

0.64 

0.63 

22 

PMSM 

0.84 

0.85 

50 

1.15 

1.12 

38 

IM 

0.47 

0.35 

21 

0.60 

0.58 

20 

DC 

0.18 

0.43 

26 

0.56 

0.75 

26 

W-  losses  to  be  dissipated  by  the  motor  under  the  most  severe  operating  condition. 
W  -  losses  at  peak  drive  power  output  and  rated  speed. 


The  PMSM  is  considered  widely  considered  as  the  ideal  candidate  for  EV  propulsion  due  to 
its  high  efficiency  in  the  constant  torque  region  [9].  It  is  clear  from  our  analysis  that  the  PMSM 
has  more  losses  as  compared  to  other  motors  when  applied  for  this  EV  propulsion  application, 
since  it  has  to  traverse  a  large  constant-power  region  during  operation.  The  losses  increase  due  to 
the  high  currents  required  for  flux  weakening  in  the  constant  power  region.  The  studies 
sununarized  in  Table  2.3  indicate  that  the  average  efficiency  of  the  PMSM  (as  indicated  by  the 
losses  W)  in  operation  is  much  lower  than  the  other  motors.  When  the  SRM  loss  contributions  are 
chosen  appropriately,  it  becomes  comparable  to  the  conventional  motors  (IM  and  DC  Motor). 
The  behavior  of  each  motor  is  discussed  in  detail  below. 


A.  PMSM  Discussion 
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The  PMSM  has  a  fixed  excitation  determined  by  the  permanent  magnets.  During  operation,  the 


PMSM  needs  to  satisfy  the  following  equations  at  a  given  power  output  and  speed  [10]; 

(2.7) 


(2.8) 


where  is  the  back-emf  induced  due  to  the  fixed  excitation  of  the  permanent  magnets,  and  denote 
respectively  the  direct  and  quadrature  axis  reactance  values  of  a  phase,  and  denote  respectively  the 
direct  and  quadrature  axis  components  of  the  phase  current,  and  V  is  the  terminal  voltage  .  To  achieve  a 
long  constant-power  region,  it  is  necessary  to  weaken  the  excitation  using  the  d-axis  component  4  of  the 
armature  current  This  leads  to  the  armature  current  increasing  above  the  rated  speed  with  an  increase  in 
speed.  The  winding  current  at  the  maximum  motor  speed  rises  close  to  or  greater  than  the  current  at  the 
peak  power  output  and  rated  speed  irrespective  of  the  motor  power  output  In  turn,  the  copper  losses 
increase  at  high  speeds  at  all  power  output  levels.  The  required  q-axis  current  slightly  decreases  due  to  the 
contribution  by  the  reluctance  torque.  The  core  losses  increase  with  speed  because  flux  weakening  in 
proportion  to  the  motor  speed  is  not  exactly  achieved. 

Die  result  is  that  the  copper  and  core  losses  at  the  maximum  drive  speed  are  high  for  very  high  ratios 
of  maximum  to  rated  speed  (long  constant-power  region).  In  the  PMSM,  it  is  difficult  to  reduce  the  loss 
dissipation  duty  substantially  by  the  choice  of  the  individual  components  of  the  total  loss.  This  is  because 
the  iron  and  copper  losses  have  similar  increasing  contributions  as  a  function  of  the  vehicle  speed.  In 
addition,  efforts  to  reduce  the  d-axis  component  of  current  through  control  result  in  increased  fluxing  in  the 
core  leading  to  an  increase  in  core  losses. 

We  have  considered  a  PMSM  with  interior  permanent  magnet  rotor  having  = 0.80  pu,  X^ = 0.50  pu 
and  X^ = 0.66  pu,  with  c  =  0.4  and /=  0.333. 

The  controller  is  implemented  following  the  concepts  of  [10],  in  order  to  determine  the  winding 
currents  as  function  of  the  vehicle  speed  and  motor  power  output.  Till  rated  speed,  the  controller  sets  the 
terminal  voltage  proportional  to  the  frequency.  Above  this  speed,  the  terminal  speed  is  held  constant.  Die 
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continuous  power  rating  is  determined  by  the  duty  of  cruising  at  the  maximum  vehicle  velocity.  For  this 
duty,  the  winding  current  is  close  to  the  winding  current  corresponding  to  peak  power  output  at  rated  speed 
for  =  2  and  is  higher  than  this  current  at  =  5  even  though  the  motor  power  output  is 

only  19  kW.  The  iron  losses  are  also  higher  than  those  at  rated  speed  (at  which  the  continuous  power  rating 
is  specified). 

B.  IM  Discussion 

The  IM  is  assumed  to  operate  under  simple  V^control.  Hence,  it  operates  at  substantially  constant  flux 
upto  rated  speed.  The  flux  is  then  inversely  proportional  to  speed  in  the  constant  power  region.  At  a  given 
power  level  and  above  the  rated  speed  of  the  motor,  the  rms  armature  current  remains  almost  constant  with 
changes  in  speed.  Also,  the  core  losses  reduce  after  the  rated  speed  of  the  motor  because  the  flux  reduction 
plays  a  greater  role  than  the  fiequency  increase.  The  loss-dissipation  duty  is  reduced  if  the  core  loss 
component  is  minimized  in  the  design. 

We  have  considered  an  IM  with  c  =  0.4  and/=  0.333.  The  continuous  power  rating  is  determined  by 
the  duty  of  maintaining  a  speed  of  45  mph  in  a  6  %  gradient. 

Since  this  speed  is  greater  than  the  rated  speed  of  the  motor  at  both  speed  ratios,  the  core  losses  are 
lower  than  those  at  rated  speed.  The  winding  current  is  approximately  equal  to  the  ratio  of  the  power 
required  at  the  above  duty  over  the  peak  power  output  for  both  speed  ratios.  The  choice  of  the  individual 
components  has  a  weak  influence  on  the  continuous  power  rating  of  the  motor. 

C.  DC  Motor  Discussion 

The  armature  current  in  the  dc  machine  is  constant  with  speed  at  a  given  power  level  after  the  rated 
speed  of  the  motor.  The  continuous  power  rating  is  again  determined  by  the  duty  of  maintaining  a  speed  of 
45  mph  on  a  6  %  gradient.  The  winding  current  is  approximately  equal  to  the  ratio  of  the  power  required  at 
this  duty  to  the  peak  motor  power  output  for  both  speed  ratios. 

References 


24 


[1]  "1999  EV  America  Technical  Specifications",  Department  of  Energy  EV  web-site  - 

http://ev.inel.gov/fop/eva/99EVATECHSPEC.PDF 

[2]  "IEEE  Standards  Collection  -  Electric  Machinery",  1993  Edition. 

[3]  “NEMA  Standard  for  Motors  and  Generators”,  NEMA  MG  1-1993. 

[4]  “EPA  Dynamometer  Driving  Schedules”,  Environmental  Protection  Agency  web-site  - 

http://www.epa.gov/OMSWWW/labmthod.htm. 

[5]  M.  Ehsani,  K.  M.  Rahman,  H.  A.  Toliyat,  “Propulsion  System  Design  of  Electric  and  Hybrid 
Vehicles”,  IEEE  Transactions  on  Industrial  Electronics,  Vol.  44,  No.  1,  pp.  19-27,  February  1997. 

[6]  T.  J.  E.  Miller,  Switched  Reluctance  Motors  and  their  Control  Oxford  Science  Publications,  1993. 

[7]  T.  J.  E.  Miller,  Y.  Hayashi,  “A  New  Approach  to  Calculating  Core  Losses  in  the  SRM”,  IEEE 
Transactions  on  Industry  Applications,  Vol.  31,  No.  5,  pp.  1039-1046,  September/October  1995. 

[8]  P.  N.  Materu,  R.  Krishnan,  “Estimation  of  Switched  Reluctance  Motor  Losses”,  IEEE  Transactions  on 
Industry  Applications,  Vol.  28,  No.  3,  pp.  668-679,  May/June  1992. 

[9]  D.  J.  Patterson,  “High  efficiency  permanent  magnet  drive  systems  for  electric  vehicles”.  Proceedings 
of  the  1997  23rd  Annual  International  Conference  on  Industrial  Electronics,  Control,  and 
Instrumentation,  pp  391-396,  Nov  9-14  1997  New  Orleans,  LA. 

[10]  P.Pillay  (editor),  Performance  and  design  of  permanent  magnet  ac  motor  drives.  IEEE  tutorial  course 
presented  at  the  IEEE  Industrial  Applications  Society  Annual  Meeting,  October  1989  -  Chapters  6,  8. 


25 


Chapter  Three 

DESIGN  APPROACH  UNDER  MULTI-PHASE  EXCITATION 


3.1  Introduction 

Recent  literature  has  identified  the  potential  for  operating  the  SRM  advantageously  under  multi-phase 
excitation  [1-  4].  The  operation  of  SRM  configurations  working  under  multi-phase  excitation  is  described 
in  these  literatures  in  terms  of  the  phase  inductances.  This  description,  though  useful  for  understanding  the 
operations  of  these  configurations,  cannot  be  used  when  designing  the  SRM,  given  a  specificaticm  set  to  be 
satisfied.  When  designing  electric  machines,  one  has  to  take  into  account  the  constraints  inposed  by  the 
electric  and  magnetic  loadings  of  the  machine  (i.e.,  current  and  flux  densities,  respectively).  To  this  end, 
the  Chapter  describes  the  work  carried  out  as  part  of  this  Report  that  attempts  to  obtain  a  better 
understanding  of  these  SRM  configurations  when  operated  in  the  multi-phase  excitation  mode.  FEA 
studies  based  cm  ANSYS™  [5]  were  the  means  to  obtain  and  understand  the  flux  distribution  in  each 
configuration.  Based  on  these  studies,  the  design  procedure  of  [6]  was  modified  by  incorporating 
coefficients  that  account  for  each  configuration  and  mode  of  excitation.  The  following  sections  introduce 
these  coefficients  and  briefly  describe  the  design  procedure,  which  uses  these  coefficients.  A  1-hp 
prototype  of  a  SRM  for  EV  propulsion  was  then  designed  and  built  based  on  this  procedure.  The  details  of 
this  prototype  are  included  together  with  a  discussion  of  the  obtained  measurement  results. 


3.2  Nomenclature 

Bs  Saturation  flux  density 

Lsnc  Stack  length 

yr  Rotor  yoke  thickness 

Ks  Slot-fill  factor 

Tavg .  peak  Average  peak  torque 
developed  by  the  SRM 


Dr,  Ds,  Dshefi 


tR,  ts 
kd 

TpeiA 

^peak 


Rotor,  stator  and  shaft  diameters 
Rotor  and  stator-pole  widths 
Ratio  of  rotor  to  stator  diameter 
Peak  torque  produced  by  a  pole  pair 
Peak  power 
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N 

Number  of  tums/pole 

Nr. Ns 

Number  of  rotor  and  stator  poles 

m 

Number  of  phases 

P 

Number  of  pole  pairs/phase = Ns  /m 

rpm,  (Or 

Rated  speed  in  rpm  and  rad/sec 

I.J 

Current  and  current  density 

y,E 

Voltage  and  electromotive  force 

A .  Ai »  At 

Electrical,  mechanical  and  stroke 

angles 

A?.  As 

Rotor-  and  stator-pole  arcs 

TRV 

Torque  per  unit  rotor  volume 

Subscripts  avg,  mech,  cont,  mis,  I,  DC,  T,  ST,  min,  STK  denote  average,  mechanical, 
continuous,  effective  value,  cunent,  direct  current,  torque,  stroke,  minimum,  and  stack, 
respectively. 

33  SRM  M^;netic  Circuit  and  Excitation  Modes 

When  exciting  a  phase  of  the  SRM,  the  flux  flows  through  the  stator  poles,  the  rotor  poles,  the  rotw 
body  and  traverses  the  air  gap.  This  flux  closes  its  path  by  returning  through  the  depth  of  the  stator  and 
rotor  bodies  behind  the  poles  known  as  the  back  iron  (see  Figure  3.1).  When  multiple  phases  are  excited,  it 
is  the  level  of  saturation  of  the  back  iron  that  primarily  determines  the  extent  to  which  the  multi-phase 
excitation  can  be  supported  by  the  SRM  while  still  getting  advantage  due  to  such  an  operation.  Based  on 
the  current  polarity  with  which  adjacent  phases  are  excited,  the  pattern  of  flux  flow  in  the  SRM  magnetic 
circuit  changes  and  determines  the  saturation  level  of  the  back-iron.  FEA-based  flux  plots  obtained  under 
the  two  main  excitation  modes  are  displayed  in  Figures  3.2  and  3.3  to  illustrate  this  better.  These  are  the 
long-flux-path  mode,  in  which  the  two  adjacent  excited  phases  produce  the  same  direction  of  flux  into  the 
back  iron  resulting  into  the  same  fluxing  pattern  as  the  single-phase  excitation  case  (see  Figure  3.2).  The 
flux  flowing  through  the  back  iron  is  hence  increased.  The  rotor  and  stator  back-iron  width  should  be 
increased  adequately  to  carry  this  additional  flux  without  reaching  deep  saturation  in  order  to  get  advantage 
in  terms  of  torque  density  when  using  multi-phase  excitation  in  this  mode. 
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Figure  3.1  Flux  distribution  for  a  8/6  four-phase  SRM  in  the  single-phase  excitation  case 

for  a  rotor  position  of  0  =  -15°. 


Figure  3.2  Flux  distribution  for  a  8/6  four-phase  SRM  in  the  multi-phase  excitation  case 
and  long-flux-path  mode  for  a  rotor  position  of  0  =  -15°. 


By  deep  saturation,  we  mean  that  the  back  iron  is  carrying  so  much  flux  that  any  increment  in  phase 
current  causes  a  negligible  further  increment  in  flux.  The  other  excitation  mode  is  the  short-flux-path  mode 
in  which  the  two  excited  phases  have  opposite  direction  of  flux  flow  into  the  back  iron.  The  flux  flowing  in 
the  back  iron  is  not  increased  substantially.  Thus,  this  mode  requires  no  increase  in  the  depth  of  the  back 
iron  and  at  the  same  time  it  supports  the  advantages  that  can  be  obtained  due  to  multi-phase  operation  (see 
Figure  3.3).  These  two  excitation  modes  increase  the  torque  density  of  the  SRM  by  increasing  the  number 
of  pole-pairs  participating  simultaneously  in  torque  production.  However,  the  long-flux-path  mode  requires 
a  corresponding  increase  in  the  back-iron  depth.  When  the  number  of  phases  is  even,  the  short-flux-path 
mode  requires  that  the  converter  have  the  capability  of  supplying  bipolar  currents  to  maintain  short-flux- 
path  mode  continuously.  When  the  number  of  phases  is  odd,  only  unipolar  currents  are  required  and  the 
short-flux-path  mode  operation  can  be  obtained  through  appropriate  winding  connections.  These  points 
will  be  explained  in  more  detail  later  in  this  chapter. 


Figure  3.3  The  flux  pattern  in  a  8/6  four-phase  SRM  in  the  multi-phase  excitation  case  and 
short-flux-path  mode  for  a  rotor  position  of  6  =  -15°. 


There  are  other  excitation  modes  (like  the  full  and  fractional  pitch  (FRSRM)  winding 
configurations)  that  attempt  to  increase  the  torque  density  of  the  SRM  by  increasing  the  number 
of  phases  acting  simultaneously  on  a  pole  pair.  As  a  result,  the  peak  mmf  to  be  produced  by  a 
phase  reduces.  However,  the  time  interval  for  which  each  phase  is  carrying  a  current  is  increased. 
So,  the  rms  phase  current  remains  almost  the  same.  This  diminishes  the  advantages  that  can  be 
obtained  by  these  modes  of  operation;  namely,  the  full  and  fractional  pitch  (PTISRM)  winding 
configurations.  Figures  3.4  and  3.5  (a)  and  (b)  show  the  mmf  and  flux  distribution  in  a  FRSRM 
carrying  square  wave  currents  [3].  In  this  case,  the  mmf  acting  on  a  pair  of  poles  is  a  function  of 
currents  flowing  in  two  phases.  In  an  equivalent  short-pitch  winding  wound  over  pole  1,  a  current 
having  a  magnitude  of  (4  -  h)  flows  and  determines  the  torque  produced  by  the  equivalent  short- 
pitch  winding. 

In  case  of  pole  2,  the  current  flowing  is  (4  -  4)-  For  pole  3,  the  current  flowing  is  (4  -  4)- 
Figures  3.5  (a)  and  3.5(b)  shows  the  magnetic  flux  distribution  in  the  FRSRM  for  square-wave 
currents  when  6  =  0°  and  0=15®  respectively.  At  0  =  0°,  pole  1  produces  positive  torque  and  pole 
3  produces  a  torque  equal  in  magnitude  to  that  of  pole  1,  but  negative  in  sign.  Pole  2  does  not 


Figure  3.4.  The  mmf  distribution  in  a  section  of  the  FRSRM. 


produce  any  torque. 
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Figure  3.5(b).  Flux  distribution  in  the  FRSRM  for  square-wave  currents  and  6  =  15®. 


At  0  =  15®,  the  currents  in  the  equivalent  short-pitch  windings  remain  the  same.  Pole  3 
produces  positive  torque,  pole  1  and  2  do  not  produce  any  torque.  The  factor  relating  the  rms 
phase  current  to  the  peak  phase  current  will  be  greater  in  this  case  since  all  the  three  phases  are 
carrying  current  at  all  times.  The  work  of  this  report  has  modified  the  design  equations  of  [6]  to 
account  for  the  various  configurations  so  that  for  a  particular  SRM  configuration  operated  under 
multi-phase  excitation,  the  design  procedure  would  account  for: 

1 )  The  increased  number  of  pole-pair  participation. 

2)  The  increased  mmf  acting  on  a  pole-pair. 

3)  The  torque  waveform  obtained  using  a  particular  configuration. 

4)  The  factor  relating  the  average  torque  to  the  peak  torque. 

As  mentioned  earlier,  the  short-flux-path  mode  of  operation  would  give  the  most  benefit  in 
terms  of  torque  density.  Thus,  the  design  procedure  has  been  developed  to  consider  only  this 
mode  of  operation.  For  instance,  the  design  procedure  does  not  consider  the  required  increment  in 
the  depth  of  the  back  iron,  which  is  essential  in  the  long-flux-path  mode  of  operation.  FEA 
analyses  were  carried  out  for  several  important  SRM  configurations  to  identify  the  major  flux 
paths  in  each  case.  Sample  ANSYS™  FEA  code  is  included  in  Appendix  C. 

Hgure  3.6  (a)  shows  that  in  the  case  of  the  four-phase  8/6  SRM  (even  number  of  phases),  at  the  end 
of  4  strokes,  phase  current  reversal  is  required  to  continue  operation  in  the  short-flux  path  mode.  However, 
Figure  3.6(b)  illustrates  that  in  the  case  of  the  five-^haselO/8  SRM  (odd  number  of  phases),  the  original 
current  direction  in  the  first  phase  can  still  sustain  operation  in  the  short-flux-path  mode  at  the  end  of  5 
strokes.  In  this  case,  appropriate  connection  of  the  phases  is  only  required  to  sustain  short-flux-path 
operation.  If  operation  is  not  sustained  in  the  short-flux-path  mode,  the  flux  would  take  the  long  path  every 
m  strokes  resulting  into  a  reduction  in  the  instantaneous  torque,  increased  losses,  noise  and  vibrations. 
Thus,  the  converter  should  have  the  capability  of  supplying  bipolar  currents  when  the  number  of  phases  is 
even  with  the  short-flux-path  mode  of  operation  (see  Figure  3.6).  Tte  flux  in  any  SRM  operating  under  the 
long-flux-path  mode  will  take  the  short  path  mode  every  m  strokes. 
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Direction  of  rotor 
revolution 


Due  to  this,  current  reversal  is 


required  in  all  the  phases  after 


ever^  m  strokes  (m  =  number  of 
phases)  to  sustain  operation  in  the 
short  ftuK  path  mode,  mhen 
number  of  phases  is  even. 


Direction  of  rotor 
revolution 


Figure  3.6  Winding  connections  for  the  short-flux-path  mode  over  one  revolution:  (a)  Four-phase  8/6 
SRM  (change  in  the  flux  pattern),  and  (b)  Five-phase  10/8  SRM  (no  change  in  the  flux  pattern). 


The  results  are  summarized  in  Tables  3.1.  It  is  clear  that  some  configurations  increase  the 
number  of  pole  pairs  participating  in  the  torque  production  while  others  increase  the  number  of 


phases  acting  on  a  pole  pair. 

Table  3J  FEA  of  SRM  ConHgurations  -  Part  1. 


SRM  Configuration 
and  Operation  Mode 


Three-phase  6/4 


Four-phase  8/6  in 
short-flux-path 
mode 


Five-phase  10/8  in 
short-flux-path 
mode 


FEA  Analysis  and  Major  Flux  Path 


XT, 


Ideal  Torque  Waveform 


^  peak 


rirLFiiu' 

kpeak 


■  ►1^  p.  p,  -►!  e 


2T, 


peak 


■*  peak 


1 

3.4  Design  G)efficients  for  Multi-phase  Operation 

This  Section  introduces  the  coefficients  enabling  the  design  of  the  SRM  under  multi-phase  excitation. 
These  coefficients  are  useful  for  extending  established  design  procedures  on  SRMs  operating  under  single- 
phase  excitation  (such  as  the  one  described  in  [6]).  It  must  be  noted  that  some  of  the  equations  listed  in  [6] 
are  only  valid  for  the  specific  case  of  the  three-phase  6/4  SRM  under  single-phase  excitation. 

It  is  well  know  that  the  SRM  operates  in  a  sequence  of  strokes.  Under  single-phase  excitation, 

each  phase  may  conduct  for  a  stroke  angle  given  by:  0  =  (3.1) 

In  some  SRM  configurations  (e.g.,  the  8/6  and  10/8  considered  in  this  Report),  the  stator- 
pole  arc  ps  is  greater  than  the  stroke  angle.  As  a  result,  only  the  fraction  Osj/Psof  the  stator-pole 
width  is  actually  utilized  in  torque-production  under  single-phase  operation.  To  account  for  this,  we 
introduce  the  stroke  angle  coefficient  that  is  given  by: 

K  =  Aj—  (3.2) 

^  ST  a 

u  57. 

This  can  be  easily  understood  by  considering  the  flux-linkage  variation  in  Figure  1.4(a)  along  the  12  A 
line.  In  this  case,  one  can  realize  that  limiting  A6  to  values  smaller  than  30°  will  yield  smaller  flux-linkages 
and  energy  conversion.  With  this  background,  we  present  the  equations  for  calculating  the  number  of  turns 
per  pole,  N,  the  peak  torque  per  pole  pair,  Tpeak  (which  is  used  for  sizing  the  rotor  dimensions  [6]),  the  peak 
winding  current,  /peat  and  the  rms  winding  current,  While  doing  this,  we  will  introduce  other  design 
coefficients. 

A.  Number  of  tums/poleN 

The  magnetic  flux  per  pole  is  given  by 

fs  bsTK  .  (3.3) 

^  ST 
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The  coefficient  in  the  denominator  of  (3.3)  accounts  for  the  fact  that  only  the  fraction  ( 1/  Ksr)  of  the 
total  stator  pole  area  (fsX  Lsnd  is  actually  used  in  some  configurations  because  the  stroke  angle  is  smaller 
than  the  stator-pole  arc.  From  (3.3)  the  flux  linkage  for  two  coils  in  series  is  given  by 

X  =  1N^  .  (3.4) 

The  DC  voltage  Vdc required  to  produce  the  magnetic  flux  in  the  core  is  calculated  as  follows  [6]: 


Vpc  =  {dXldO)  0)  = 


^  A  ^ 


to. 


(3.5) 


Substituting  (3.3)  and  (3.4)  into  (3.5)  and  solving  for  the  number  of  tums/pole  yields: 


N  =  Kst  X 


n  xy 


DC 


CO  X  m  X  1 J  X  Lsj-k  X  N  Bg 


(3.6) 


where  all  linear  dimensions  are  in  m,  Vdc  is  in  volts,  Bj  is  in  weber/m^  and  (O  is  in  rad/sec. 


N  =  IC  X _ 46500  xVpc _ .  (3.7) 

^  rpm  X  mxtg  X  x  N  gX  Bg 

When  the  coils  on  nrore  than  one  pole  pair  are  connected  in  series  to  form  a  phase,  (3.4)  must  be 
divided  by  (P/2).  The  constant  46500  is  a  conversion  factor  that  arises  when  all  linear  dimensions  are 
expressed  in  inches  and  the  speed  value  in  rpm  as  used  in  [6]. 


B.  Peak  torque  per  pole  pair 

The  average  peak  torque  is  related  to  the  peak  power  and  the  rated  speed  a  at  which  the  torque  has  to 
be  developed  as  follows: 

T  =  (3.8) 

avg  ,peak  ^ 

In  the  previous  Section  (see  Tables  3. 1),  we  noted  that: 

•  Only  one  pair  of  poles  develops  torque  in  the  three-phase  6/4  SRM  at  any  given  time,  and 

•  The  average  torque  developed  in  the  five-phase  10/8  SRM  is  twice  of  that  developed  by  that  one 
pole  pair. 

Thus,  the  peak  torque  to  be  developed  by  a  pole  pair  in  the  10/8  SRM  for  a  given  average  peak  torque 
is  half  of  that  to  be  developed  by  one  pole  pair  in  the  6/4  SRM. 
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To  account  for  this  in  the  design  equations,  we  introduce  the  torque  coefficient  that  accounts  for  the 
reduction  in  the  developed  torque  by  only  one  pole  pair  in  a  particular  configuration.  As  per  our  definition, 
Kt  =  1  for  the  6/4  SRM  and  /sTj-  =  2  for  the  10/8  SRM. 

Up  to  this  point,  we  only  considered  the  unsaturated  operation  of  the  SRM.  However,  the  actual  SRM 
operates  under  saturation.  Thus,  we  introduce  the  coefficient  Ki^eak  accounting  for  the  effects  of  saturation 
and  given  as  the  ratio  of  the  peak  torque  to  the  average  value  obtained  from  FEA  analysis;  that  is  Kjpeak  = 


Tpeak^mg- 

hi  general,  the  average  peak  torque  developed  by  the  SRM  is  given  by 


Kr  X 


avg  ,peak 


^  Tpeak  ^  ^  srr 


-X  T 


peak 


(3.9) 


From  (3,8)  and  (3.9),  the  peak  torque  to  be  developed  by  only  one  pole  pair  is  given  by: 


^  ^  Tpeak  ^  ^  ST  ^ peak 

~~  <0 
2 


(3.10) 


C.  Peak  winding  current  Ipeak 

In  a  three-phase  6/4  SRM  with  short-pitch  windings  under  single-phase  excitation  (that  is, 
with  only  one  pole  pair  producing  torque),  the  peak  winding  current  for  developing  the  torque 
Tpeaki&  given  by  [6]: 

I  _ _ _ .  (3.11) 

N  X  Bs  X  DgX  Lstk 

In  SRM  with  full-  and  fractional-pitch  windings,  the  winding  arrangement  results  in  increased 
mmf  acting  on  the  flux  paths.  As  a  result,  the  winding  current  can  be  reduced  for  a  desired  flux 
density  to  be  developed  in  the  air-gap.  To  account  for  this,  we  introduce  the  coefficient  Kj  which 
is  the  amount  by  which  the  winding  current  can  be  reduced  as  compared  to  the  short-pitch 
winding  case.  By  this  definition,  K/  =1  in  the  6/4  and  8/6  SRM,  and  Ki  =  2  in  the  12/8  SRM  with 
full-  or  fractional-pitch  windings.  Thus,  the  peak  winding  current  can  be  calculated  as: 
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(3.12) 


^  _  1  W _ ^peak _ 

K,  NxBsXDpXLsTK 

This  equation  applies  when  the  coils  on  the  pole  pairs  of  a  phase  are  connected  in  series.  For 
the  coil  parallel  connection,  (3.12)  should  be  multiplied  by  (P/2). 

In  summary,  we  can  calculate  the  tums/pole,  N,  the  peak  torque  developed  by  only  one  pole 
pair,  Tpeak  and  the  peak  winding  current,  Ipeak  using  the  new  coefficients  to  modify  the  equations  of 
the  design  procedure  in  [6].  For  sizing  the  winding  conductors,  we  need  to  calculate  the  rms 
value  of  the  winding  current  Inns  as  shown  below. 


D.  RMS  winding  current  I  rms 

The  coefficient  Ka  is  defined  as  the  ratio  between  the  rms  and  peak  winding  currents.  This 
coefficient  takes  different  values  depending  on  the  number  of  phase  currents  applied  at  any  one  time  and 
the  SRM  configuration.  For  continuous  single-phase  operation,  this  ratio  is  given  by  [6]: 

(3.13) 

since  there  is  a  current  in  only  one  phase  for  the  duration  of  the  stroke  angle  Ogr-  In  the  8/6  and  10/8  SRMs 
under  multi-phase  excitation,  the  current  pulse  exists  for  the  duration  of  the  stator-pole  arc  Ps-  This  increase 
in  conduction  (and  hence,  rms  winding  current  value)  can  be  accounted  for  by  introducing  the  coefficient 
Kgr  in  (3.12).  For  exan^le,  Ps  =  18°  and  0st  =  9°  in  a  10/8  SRM  so  the  period  of  the  current  pulse  is 
increased  by  a  factor  of  18/9  =  2.  Likewise,  Ps = 22.5°  and  0st  =  15°  in  the  case  of  the  8/6  SRM  with  short- 
pitch  windings  under  multi-phase  excitation;  the  period  of  the  current  pulse  increases  by  a  factor  of  22.5/15 
=  1.5. 

In  general,  this  coefficient  for  SRM  with  short-pitch  windings  becomes: 
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In  the  full-  and  fractional-pitch  winding  SRM,  is  given  by  the  square  root  of  the  ratio  between  the 
number  of  phases  conducting  at  a  time  and  the  total  number  of  phases.  For  the  three-phase  full-pitch 


winding  SRM  having  two  phases  conducting  at  any  dme, 


For  the  three-phase  fractional- 


pitch  windmg  SRM,  all  three  phases  conduct  at  any  given  time  and  k 


The  multi-phase  design  coefficients  introduced  in  this  section  are  listed  in  Table  3.2  for  the  considered 
SRM  configurations. 


3.5  SRM  Design  Steps  and  Prototype  Design 

Using  the  above  design  coefficients,  a  simple  method  for  designing  SRM  under  multi-phase 
excitation  was  developed  in  [7].  The  method  ideas  were  incorporated  into  a  computer  program 
that  calculates  the  SRM  main  dimensions,  peak  current,  total  weight  and  winding  losses  given  the 
input  specifications.  Figure  3.7  shows  the  program  flowchart  which  is  described  in  more  detail  in  Dr.  E. 
HaU  Mitre’s  doctoral  report  [7]. 

From  Chapter  2 ,  the  SRM  satisfying  the  assumed  performance  specifications  must  have  the  following 
ratings: 

Peak  Power  Rating  =  34  kW 

Continuous  Power  Rating  =  23  kW 
The  other  specifications  of  the  SRM  were  assumed  as: 

DC  Supply  Voltage  =  300  V 

Rated  Speed  =  1200  rpm 

Maximum  Speed  =  6000  rpm 


41 


Table  3.2.  Design  coefficients. 


Ns /Nr 

Winding  Scheme 

Kt 

Kst 

mm 

6/4  or  12/8 

Short  pitch 

1.0 

1.0 

1.0 

0.5778 

6/4  or  12/8 

Full  pitch 

1.0 

■umi 

2.0 

0.8165 

6/4  or  12/8 

Fractional  pitch 

1.0 

2.0 

1.0000 

8/6 

Short  pitch,  short  flux  path 

1.5 

MEM 

1.0 

0.6124 

10/8 

Short  pitch,  short  flux  path 

2.0 

1  2.0  1 

1.0 

0.6324 

Figure  3.7  Flowchart  of  the  computer  program  that  calculates  the  SRM  main  dimensions,  peak  current,  total  weight  and  winding  losses. 


The  maximum  speed  was  chosen  as  6000  rpm  considering  the  limitations  of  the  TMS320C240  EVM. 
The  rated  speed  was  chosen  as  1200  tpm  based  on  a  choice  of  the  speed  ratio  of  5  for  the  drive.  The 
different  SRM  candidates  are  those  ones  listed  in  Table  3.2.  For  each  SRM,  Table  3.3  lists  the  stator  outer 
radius  and  stack  length,  weight,  peak  winding  current  and  copper  losses. 
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These  results  were  obtained  using  the  same  constraints  on  maximum  flux  and  current  densities  for  all 
configurations.  It  is  clear  that  the  8/6  SRM  and  the  10/8  SRM  have  almost  identical  weights. 
However,  the  8/6  SRM  has  a  greater  angle  range  available  for  tum-on  angle  advancement  to  build 
up  the  phase  current  even  at  high  speeds.  As  a  result,  its  torque  characteristic  is  better  at  high 
speeds.  It  is  noted  that  electric  motors  for  EV  must  operate  over  a  wide  speed  range.  The  FRSRM 
and  the  full-pitch  winding  SRM  have  highly  coupled  windings,  and  hence,  the  phase  currents  are 
difficult  to  set  up  and  control.  Due  to  these  reasons,  the  8/6  SRM  configuration  was  chosen  for  this 
application  since  it  has  the  feasibility  for  continuous  operation  in  the  short-flux-path  nxxle  and  adequate 
angle  control  for  phase  advancement.  The  8/6  SRM  design  was  then  scaled  using  the  following  scaling 
principle  proposed  in  [8]: 


^  avgypeak 


V  ^  ; 


(3.15) 


FEA  studies  were  carried  out  to  refine  the  design.  Some  changes  to  the  geometry  were  made  such  as 
the  manufacturing  requirement  of  straight  poles  for  easy  insertion  of  coils  wound  outside  the  stator  on  a 
former.  The  drawings  for  this  SRM  were  prepared  so  that  it  could  be  manufactured  by  Baldor  Motor  and 
Drives,  Fort  Smith,  AR.  Table  3.4  shows  the  specifications  of  the  manufactured  SRM.  The  geometry  of  the 
stator  and  rotor  laminations  is  shown  in  Figure  3.8. 

Table  3.3.  Design  Results  [7] 


Ns/ 

Nr 

Winding 

Scheme 

Excitation 

Type 

Stator 
Outer 
Radius  (m) 

Stack 

Length  (m) 

Weight 

(kg) 

Peak 

Current  (A) 

Copper 
Losses  @ 

Front  (Watts) 

6/4 

Short 

Single  phase 

0.232 

0.298 

306 

258 

665 

12/8 

pitch 

0.188 

133 

264 

646 

12/8 

Full  pitch 

Multi  phase 

0.172 

111 

274 

773 

12/8 

Multi  phase 

0.179 

0.237 

120 

274 

802 

8/6 

Multi  phase 

0.185 

0.237 

145 

132 

563 

10/8 

(short  flux  path) 

0.186 

0.237 

139 

132 

688 
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Table  3.4  1-HP  Prototype  SRM  Speciflcations 


Rated  Speed 

1200  rpm 

Maximum  Speed 

6000  rpm 

Continuous  Power  Rating 

490  W  @  Rated  Speed 

DC  Bus  Voltage 

75  V 

Winding  Resistance/Phase  @  75  °C 

0.43  Ohms 

Stator  Outer  Radius 

81.8  mm 

Stack  Length 

72  mm 

Winding  Conductors 

68  turns  of  2  conductors  of  AWG  20  wire  on  each  pole. 

Peak  Winding  Current 

15  A 

Multiple  Excitation  Mode 

Two-Phase,  Short-Flux-Path  Mode 

Stator  lamination  Rotor  lamination 


Figure  3.8.  Stator  and  rotor  lamination  cross-section. 


3.6  Description  Of  the  Experimental  Results 

A.  Comparison  of  Measurements  vs.  FEA  Simulations 

The  prototype  SRM  was  characterized  using  2D  FEA  studies.  A  deviation  was  expected 
between  the  flux-linkage  characteristics  from  FEA  and  those  from  actual  measurements  at  the 
unaligned  position.  From  this  point  of  view,  a  20  %  margin  was  incorporated  between  the  desired 
peak  power  rating  and  the  actual  design  value.  The  single-phase  measurements  are  based  on  the 
procedure  described  in  [6],  while  the  multi-phase  measurements  are  based  on  the  new  approach 
presented  in  Chapter  4  for  the  multiply-excited  6/4  and  8/6  SRM.  The  flux  data  were  obtained  by 
integrating  the  back  emf  equation  using  the  measured  terminal  voltages  and  phase  currents. 
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The  phase  resistance  was  measured  using  a  digital  multi-meter  as  well  as  by  the  V/I  method. 
The  latter  was  also  used  for  resistance  measurements  in  order  to  determine  the  average  winding 
temperature  rise.  The  torque  was  measured  using  a  Lebow  torque  meter/transducer. 

Measurements  were  also  carried  out  for  the  torque  characteristics  under  single-phase 
excitation  as  well  as  multi-phase  excitation  (short  and  long  flux  path  modes).  The  torque 
characteristics  under  the  various  modes  are  included  and  discussed  in  detail  in  Chapter  Four. 
Figures  3.9  (a)  and  (b)  show  the  comparative  torque  and  flux-linkage  characteristics  under  single¬ 
phase  excitation  obtained  using  FEA  and  measurements.  The  observed  difference  between  the 
FEA  and  measured  characteristics  deviates  from  the  usually  expected  one  of  difference  occurring 
in  2D  FEA  results  from  the  actual  flux-linkage  characteristics  at  the  unaligned  position.  In  this 
case,  a  comparatively  greater  deviation  is  observed  at  the  aligned  position.  This  pattern  of 
deviation  could  be  attributed  only  to  a  difference  between  the  specified  and  actual  coil  turns  of 
about  20%.  Also,  a  phase  resistance  measurement  shows  that  it  is  less  than  the  calculated  value  of 
0.43  ohms  for  a  winding  conductor  of  gage  20  by  about  20  %.  The  actual  resistance  value  is  about 
0.34  ohms.  However,  this  possible  reduced  number  of  turns  in  the  prototype  could  not  be  verified 
with  the  manufacturer. 

Figure  3.10  shows  a  plot  of  the  measured  flux-linkage  characteristics  for  the  8/6  SRM  under 
single-phase  excitation.  Figure  3.11  shows  the  measured  average  winding  temperature  rise  when 
a  phase  is  carrying  the  peak  phase  current  of  15  A.  This  verifies  the  short-time  overload 
capability  of  the  designed  SRM,  since  the  temperature  rise  during  this  overload  is  well  within  the 
allowable  temperature  rise  for  the  Class  F  insulation  used  in  the  machine. 
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Figure  3.10.  Flux-linkage  surface  for  the  8/6  SRM  prototype  under  single-phase  excitation. 


Fig.  3.11  Temperature  rise  in  the  windings  at  the  peak  phase  current  of  15  A. 


Chapter  Four 

ARTIFICIAL  NEURAL  NETWORKS  FOR  MODELING 
THE  MULTIPLY  EXCITED  SRM 

4.1  Introduction 

This  chapter  presents  a  new  modeling  technique  for  multi-phase  excited  SRM  that  uses  FF  ANN  to 
account  for  the  mutual  flux  interaction  among  excited  phases  and  to  compute  the  torque,  while  requiring  a 
smaller  measured  data  set  as  compared  to  earlier  modeling  approaches. 

In  recent  years,  the  need  for  multi-phase  excitation  modeling  techniques  has  increased  greatly  [1]  since 
it  has  been  prqwsed  that  multi-phase  excitation  can  substantially  inprove  the  performance  of  SRM  in 
terms  of  torque  ripple,  torque  density,  acoustic  noise  and  efficiency.  Moreover,  studies  on  sensorless 
operation  (torque  and  position)  and  advanced  control  strategies  increase  further  the  need  for  a  simple  multi¬ 
phase  model  based  on  a  small  number  of  measurements  and  capable  of  achieving  high  accuracy. 

SRM  have  been  traditionally  driven  using  single-phase  excitation;  i.e.,  exciting  only  one  phase  at  a 
time.  However,  even  a  SRM  designed  to  operate  under  single-phase  excitation,  such  as  the  three-phase  6/4 
SRM,  operates  under  multi-phase  conditions  during  the  commutation  from  one  phase  to  the  next  one.  An 
accurate  model  of  the  SRM  should  consider  the  influence  of  the  currents  in  the  other  excited  phases.  In  the 
past,  single-phase  excitation  modeling  techniques  have  been  extended  usmg  superposition  to  model  the 
SRM  multi-phase  operation  due  to  several  reasons;  namely,  they  are  sinqiler  to  inqilement,  faster  to 
compute,  and  require  a  smaUer  number  of  data  points  than  multi-phase  excitation  modeling  techniques. 
This  type  of  approximation  is  frequently  adequate  as  a  first  approach  to  model  the  SRM  dynamic  behavior 
and  to  study  sinple  control  strategies.  However,  superposition  principles  will  not  yield  high  accuracy  due 
to  the  non-linearities  of  the  system  as  well  as  the  existence  of  mutual  interactions  between  the  phases. 
Hence,  proper  multi-phase  excitation  modeling  techniques  should  be  implemented  when  requiring  higher 
accuracy. 
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In  the  past,  many  techniques  have  been  reported  for  on-line  calculation  of  the  SRM  performance 
variables  and  parameters  as  function  of  the  terminal  variables  (i.e.,  current,  position)  under  the  single-phase 
operating  mode  for  exanple  [2].  However,  it  is  difficult  to  develop  equations  to  calculate  easily  and 
accurately  these  performance  variables  and  parameters  in  the  multi-phase  excitation  case. 

The  chapter  first  presents  the  proposed  multi-phase  model  for  the  SRM  together  with  a  discussion  on 
the  behavior  of  the  phase  flux  and  net  electromagnetic  torque  under  multi-phase  excitation.  The 
measurement  techniques  used  to  obtain  the  data  to  train  the  ANN  are  then  described.  Finally,  the  chapter 
presents  the  simulation  results  obtained  using  the  new  model  and  conpares  them  with  the  results  obtained 
using  the  single-phase  model. 


4.2  The  Multi-Phase  Modeling  of  the  SRM 


The  phase  voltage  equation  of  the  multiply  excited  SRM  is  given  as  follows: 


rfA, 

V,  =/?,i,+^ 
^  ^  ‘  dt 


(4.1) 


where  Vj,  Rj,  ij,  and  Xj  are  the  applied  voltage,  resistance,  current  and  flux  linkage  of  the  j-th  phase, 
respectively. 

The  flux  linkage  of  the  j-th  phase,  Xj,  can  be  expressed  as  follows: 


X. 

J 


ph 

=X.(i.,0.)+  S  A.,  =A.//.,6'.;+A. 

]}  J  J  Jk  JJ  J  J  jm 


k=l 


(4.2) 


where  Nph  is  the  number  of  phases  simultaneously  excited,  is  the  self-flux  linkage  of  the  y-th  phase,  Xjk  is  the 
mutual-flux  linkage  between  the  j-  and  ^-th  phases,  Oj  is  the  rotor  position  of  the  y-th  phase  and  represents  the 
total  mutual-flux  linkage  interacting  with  the  y-th  phase.  The  mutual-flux  linkages  are  in  general  non-linear 
function  of  the  current  in  all  the  excited  phases  and  the  rotor  position  with  respect  to  the  phase  y. 
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Since  all  the  other  phases  and  rotor  poles  have  fixed  geometric  relationships  with  respect  to  the  phasej,  we 
only  need  to  consider  the  rotor  positions  with  respect  to  the  f'  phase.  Thus,  the  mutual-flux  linkage  is  a 
nonlinear  function  of  the  relative  rotor  position  and  the  currents  in  all  the  excited  phases.  By  using  a  model  for 
mapping  this  function,  the  mutual-flux  linkage  portion  can  be  subtracted  from  the  total  flux  hnkage  in  order  to 
obtain  the  self-flux  linkage  portion.  From  the  self-flux  linkage,  the  phase  current  can  be  calculated  using  the 
flux-linkage  characteristic  under  single-phase  excitation.  Rnally,  from  the  phase  currents,  the  total  torque  can  be 
then  calculated  using  another  non-linear  function  representing  the  full  mapping  of  the  multi-phase  torque  with 
respect  to  currents  in  all  the  excited  phases  and  the  rotor  position.  Fig.  4. 1  shows  a  generic  phase  of  the  multi¬ 
phase  model. 

The  sub-sections  below  discuss  in  more  detail  the  behavior  of  the  torque  and  mutual  flux  interactions 
between  the  phases  when  exciting  multiple  phases  as  well  as  the  multi-phase  modeling  steps  and  concepts. 


A.  Mutual-Flux  Linkage  and  Torque  Behavior 

The  mutual  flux  interactions  between  the  phases  of  a  SRM  are  determined  by  several  factors  such  as 
the  current  level  in  the  excited  phases,  the  rotor  position  and  the  saturation  level  of  the  core,  in  particular 
that  of  the  back  iron.  Also,  the  SRM  may  operate  in  either  short  or  long  flux  path  mode  depending  on  the 
polarity  of  the  currents  in  the  excited  phases.  Suppose  that  we  consider  only  two  adjacent  phases  excited 
simultaneously.  In  the  case  of  the  short-flux-path  mode,  the  fluxes  in  the  two  adjacent  phases  have  opposite 
directions  of  flow  into  the  back  iron  (see  Fig.  4.2).  Thus,  the  majority  of  the  flux  tends  to  concentrate  in  the 
region  of  the  back  iron  between  the  two  phases  (see  Fig.  4.2).  In  the  long  flux  path  mode,  the  fluxes  have 
the  same  direction  of  flow  (see  Fig.  4.2).  Thus,  the  fluxing  pattern  remains  the  same  as  the  single-phase 
excited  case.  The  excitation  of  the  two  phases  simultaneously  drives  the  back  iron  further  into  saturation 
than  in  the  single-phase  excitation  case. 
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This  chapter  will  later  introduce  the  concept  of  the  mutual  flux  interaction  function  M(ii, 
which  is  the  non-linear  fimction  determining  the  mutual  interaction  between  the  excited  phases.  The 
function  M  may  produce  either  an  increase  or  decrease  in  the  net  phase  flux  depending  on  the  current  levels 
in  the  excited  phases  //,  ii-jNph  and  the  rotor  position  6.  In  the  case  of  the  short-flux-path  mode  at  low 
current  levels,  the  change  in  the  fluxing  pattern  under  multi-phase  excitation  results  in  increasing  the  net 
flux  in  each  phase  due  to  the  decrease  in  the  reluctance  of  the  flux  path  (due  to  a  shorter  length).  At  high 
current  levels,  however,  the  net  flux  decreases  due  to  the  saturation  of  the  back  iron  between  the  two 
phases. 

For  the  long-flux-path  mode  of  the  operation,  the  net  flux  reduces  at  all  current  levels.  The  peak  level 
of  mutual  interaction  is  lower  m  the  short-flux-path  case  as  compared  to  the  long-flux-path  case.  This 
difference  in  interaction  level  between  the  two  modes  will  be  more  pronounced  when  increasing  the 
number  of  participating  phases.  As  a  result,  the  changes  in  the  mutual  flux  and  torque  are  smaller  in  the 
shoit-flux-path  case  than  in  the  long-flux-path  case.  Fig.  4.3  illustrates  these  points  through  a  plot  of  the 
function  M  at  a  rotor  position  of  -7.5°,  which  was  derived  using  ANSYS™-FEA  simulations  of  a  6/4 
SRM.  It  can  be  observed  that  M  is  positive  for  certain  current  levels  in  the  short-flux-path  mode.  For  all 
other  conditions,  it  is  negative. 

The  torque  produced  by  each  phase  is  a  function  of  the  flux  level  in  the  air-gap  portion  associated  with 
that  phase.  The  torque  produced  by  each  of  the  conducting  phases  also  changes  since  the  flux  level  of  a 
phase  changes  with  the  presence  of  currents  in  the  other  phases.  This  torque  change  follows  the  pattern  of 
the  change  in  the  flux.  Since  the  mutual  flux  interactions  are  in  general  higher  in  the  long-flux-path  mode, 
the  deviation  in  net  torque  from  that  obtained  by  superposition  of  single-phase  values  is  greater  in  the  long- 
flux-path  excitation  mode.  Hence,  the  principle  of  superposition  should  not  be  used  to  model  the  SRM 
under  multi-phase  excitation.  Similar  behavior  of  SRM  under  multi-phase  excitation  has  been  recendy 
reported  in  [3]. 


53 


0.01 


Figure  4.3.  The  variation  of  the  mutual  flux  interaction  function  M  with  phase  winding  currents  at  0  =  -7.5°. 

In  this  report,  the  case  of  Np*  =  2  has  been  considered.  Different  non-linear  function  approximation 
techniques  can  be  inplemented  in  order  to  represent  the  mutual  interaction  M(ii,i2,0)  and  total  torque 
T^ii.izO)  functions.  Most  approximation  techniques  when  applied  to  these  functions  would  require  a  large 
number  of  data  points  to  get  the  required  fit  to  the  non-linear  behavior  of  the  SRM.  Also,  these  techniques 
would  involve  tedious  calculations  (several  corr^utations  of  trigonometric  functions  and  surface  integrals 
at  each  time-step)  [4]  for  torque  calculation  that  are  difficult  to  implement  efficiently  for  simulations  and 
on-line  estimation.  For  example,  one  can  use  a  Fourier  series  like  formulation  to  express  the  variation  of  the 
M  function  with  rotor  position.  Knowing  the  single-phase  characteristics  of  the  phases,  we  need  the  M 
values  as  a  function  of  current  in  the  2  excited  phases  (surfaces  like  that  in  Fig.  4.3)  at  a  minimum  of  4  rotor 
positions  to  confute  the  coefficients  of  this  Fourier  series. 
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At  each  rotor  position,  we  need  10  steps  of  one  phase  current  x  10  steps  of  the  other  current  or  100 
measurements;  hence,  we  require  a  total  of  400  measurements.  These  nurhbers  are  based  on  ANSYS™- 
FEA  studies  aimed  at  representing  these  interaction  functions  for  several  SRM  configurations.  In  some 
cases,  data  at  more  positions  are  required  to  be  able  to  use  a  Fourier  series  representation.  For  estimating  the 
torque,  it  is  essential  to  compute  the  total  co^nergy  of  the  system  by  consideration  of  the  self  and  mutual 
fluxes.  Other  modeling  methods  must  be  sought  to  reduce  the  large  number  of  measurements  required  to 
obtain  the  model  data  and  simplify  the  computations.  In  particular,  two  FF  ANN  have  been  applied  that  are 
expected  to  require  a  substantially  reduced  data  set  and  can  be  trained  to  generalize  and  estimate  the 
behavior  of  the  mutual  interactions  and  torque. 

The  role  of  the  flux  ANN  is  to  provide  static  mapping  between  the  currents  in  the  two  excited  phases, 
the  rotor  position  and  the  mutual  flux  interaction  function.  Likewise,  the  torque  ANN  is  used  to  provide 
static  mapping  between  the  two  currents,  the  rotor  position  and  the  electromagnetic  torque  developed  in  the 
SRM. 

B.  The  Mutual  Flux  Interaction  Function  M( i i,  h,  0) 

A  multiply  excited  magnetic  system  (in  this  case,  the  SRM  under  multi-phase  excitation)  is  governed 
by  the  following  relationship  among  the  flux  linkage  A,,  the  phase  current  the  total  system  coenergy 
and  the  total  system  field  energy 

(4-3) 

In  the  case  of  a  SRM  with  two  phases  simultaneously  excited  {N^  =  2),  (4.3)  can  be  expressed  as: 

A//,+A2i2  =  W^c+lY/  (4.4) 

The  total  phase  fluxes  A;  and  A2  are  divided  into  self-flux  linkages  (Ayy  and  A22)  and  mutud-flux 
linkages  (Ay2  and  A2y)  as  follows: 
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A/  —  A/7+A;2  3nd  A2  =  A22+A2/ 


(4.5) 


Based  on  these  equations,  a  mutual  flux  interaction  function  M  can  be  defined  as  follows: 

A/2  =  M{iu  hOlh  and  A2/  =  M{ ij,  i2,6)ii  (4.6) 

Then,  (4.4)  can  be  rewritten  as 

A;;  i/+A22  i2+2M(ii,  izdflih^  Wc+Wf  (4.7) 

The  mutual  flux  interaction  function  M(i],i2,  d)  can  be  calculated  as  follows: 

(4.8) 

2i,i2 


To  obtain  the  plot  of  Figure  4.3,  Wc  Wf,  Xu  and  A22  were  obtained  at  each  current  level  and  rotor 
position  using  ANSYS™-simulations. 


C.  ANN-based  Representation  ofM  and  Torque 

The  SRM  operates  in  strokes  over  the  conduction  angle  given  by: 

0  =i^°  (4.9) 

^  mN, 

where  m  is  the  number  of  phases  and  Nr  is  the  number  of  rotor  poles.  The  positive  torque  capability  of 
each  phase  is  from  g  _  -180  °  to  0°  where  6  is  the  rotor  position  with  respect  to  a  phase  measured 

between  the  centers  of  the  rotor  and  stator  poles,  respectively.  To  use  a  FF  ANN  to  represent  the  mutual 
flux  interactions  between  the  SRM  phases,  we  consider  the  influence  of  the  leading  phase  on  the  working 
phase  (in  the  direction  of  rotation)  and  vice-versa  for  rotor  positions  d  at  the  working  phase  varying 

between  °  to  0°.  Correspondingly,  the  rotor  position  at  the  leading  phase  varies  between 

Nr 
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For  the  three-phase  6/4  SRM  considered  here,  the  flux  and  torque  ANN  will  account  for  mutual  flux 
interactions  between  the  leading  and  working  phases  as  well  as  the  electromagnetic  torque  output  for  6 
varying  between  45°  to  0°  at  the  working  phase.  The  rotor  position  with  respect  to  the  leading  phase  will 
correspondingly  vary  between  -15°  to  30°.  This  choice  of  angles  will  not  only  allow  representing  the 
mutual  flux  interactions  and  the  electromagnetic  torque  from  the  minimum  possible  tum-on  angle  at  the 
workmg  phase  but  also  allows  for  delayed  turn-off  at  the  leading  phase.  For  a  four-phase  8/6  SRM,  the  FF 
ANN  will  account  for  the  mutual  flux  interactions  within  an  angle  span  from  0=-3O°  to  0°  when  the  rotor 
position  with  respect  to  the  leading  phase  varies  between  -15°  to  15°. 


4.3  ANN  Training  using  Measurements  of  M  and  Torque 

The  FF  ANN  were  first  trained  using  data  points  computed  from  ANSYS™-2D  simulations.  This  is 
adequate  for  the  initial  stages  of  model  development,  for  determining  the  ANN  architecture  and  for 
carrying  out  initial  studies.  However,  it  is  well  known  that  2D  FEA  cannot  be  expected  to  give  an  accurate 
magnetic  characterization  of  the  SRM  because  of  neglecting  the  end-winding  effects  [5].  An  accurate 
modeling  of  the  SRM  for  verification  and  tuning  of  the  drive  performance  should  be  based  on  measured 
data.  Static  measurements  were  carried  out  on  a  three-phase  6/4  SRM  designed  for  single-phase  operation 
as  well  as  a  1-hp  four-phase  8/6  SRM  prototype,  designed  for  multi-phase  operation  in  the  short-flux-path 
mode  for  EV  propulsion  applications. 

The  multi-phase  data  were  recorded  in  the  FEA  simulations  at  6  equally  spaced  intervals  of  angular 
position  in  the  range  45°  to  0°,  5  equally  spaced  intervals  of  current  from  3  A  to  15  A  in  the  two  phases. 
The  training  data  thus  have  (6x5x5=  150  pomts).  The  single-phase  excitation  data  are  composed  by  16 
equally-spaced  rotor  positions  and  15  equally-spaced  current  phase  resulting  in  240  more  data  points. 
Remember  that  the  Fourier  series  based  representation  requires  400  measured  points  in  addition  to  the 
single-phase  excitation  data. 
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Also,  the  multi-phase  torque  would  have  to  be  computed  using  conplex  formulations  that  are  difficult 
to  implement  either  in  the  off-line  or  on-line  mode  [4].  In  the  off-line  case,  the  torque  under  multi-phase 
excitation  would  be  pre-computed  and  stored  as  3D  lookup  tables.  At  the  phase  current  and  rotor  position 
values  measured  from  the  drive,  the  torque  would  then  have  to  be  con5)uted  by  on-line  interpolations.  On 
the  other  hand,  the  on-line  case  would  be  more  difficult  to  inclement  due  to  requiring  estimation  of  flux- 
linkage  characteristics  as  well  as  torque  conqiutation  on-line. 

The  ANN  were  trained  and  their  stracture  was  optimized  using  MATLAB™.  For  the  FF  ANN 
mapping  of  the  M  function  to  the  terminal  variables  in  the  6/4  SRM  case,  the  desired  performance  accuracy 
was  obtained  using  an  ANN  having  an  input  layer  of  3  neurons,  two  hidden  layers  having  respectively  3 
and  2  neurons  (with  log-sigmoid  activation  functions),  and  an  ouqiut  layer  of  1  neuron  (having  a  linear 
activation  function).  Figure  4.4  below  shows  the  M  function  prediction  using  the  FF  ANN  at  different 
current  levels  in  the  working  and  leading  phases  as  function  of  rotor  position.  It  is  noted  that  the  M  function 
is  always  negative  meaning  that  the  phase  flux  will  always  be  reduced  in  the  6/4  SRM  operating  under 
multi-phase  excitation  and  long-flux-path  nxxle.  Similarly,  the  torque  ANN  has  an  input  layer  of  3  neurons, 
a  hidden  layer  of  7  neurons  and  an  output  layer  of  1  neuron.  Figure  4.5  shows  the  topology  of  the  torque 
prediction  FF  ANN.  Figure  4.6  shows  the  electromagnetic  torque  prediction  using  the  FF  ANN  at  different 
current  levels  in  the  working  and  leading  phases  as  function  of  the  rotor  position.  The  torque  values  were 
obtained  using  FEA. 

In  order  to  take  the  static  measurements,  a  test/measurement  bench  was  set  up  as  illustrated  in  Figure 
4.7  that  also  illustrates  conceptually  the  electrical  circuit  for  the  measurements.  The  torque  measurenfients 
were  made  using  a  Lebow  torque  meter  at  9  rotor  positions  and  5  levels  of  phase  current  at  equally  spaced 
intervals  resulting  in  5  x  5  x  9  =  225  measurements  for  the  multi-phase  excitation  case.  The  single-phase 
excitation  case  was  characterized  using  16  x  5  =  80  measurements;  that  is  a  total  of  305  measurements. 

Rgure  4.8  shows  the  static  torque  measurements  on  the  four-phase  8/6  SRM  when  operated  in  the 
short-  and  long-flux-path  modes. 
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Rotor  position  with  respect  to  working  phase  (Degrees) 


Figure  4.4  The  output  of  the  FF  ANN  for  estimating  the  mutual  flux  interaction  function  at 
different  current  levels  as  function  of  the  rotor  position  with  respect  to  the  working  phase. 


Figure  4.5  The  FF  ANN  used  for  torque  mapping. 


Figure  4.6  The  output  of  the  FF  ANN  for  estimating  the  torque  at  different  current  levels 
as  function  of  the  rotor  position  with  respect  to  the  working  phase. 
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Figure  4.7.  Test/measurement  bench  for  torque  and  mutual-flux  linkage  measurements. 


Since  this  SRM  was  designed  to  take  advantage  of  the  shoit-flux-path  mode  concept,  the  torque  output 
is  diminished  in  the  long-flux-path  mode  due  to  deep  saturation  of  the  back  iron.  The  torque  value  that 
would  be  obtained  by  using  superposition  of  the  torque  of  the  working  and  leading  phases  operated 
independently  is  also  shown  to  substantiate  the  fact  that  superposition  cannot  be  used  for  modeling  the 
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multiply  excited  SRM.  fo  addition,  the  single-phase  excitation  characteristics  for  the  leading  and  working 
phases  are  shown  for  con^leteness.  The  measured  torque  data  were  used  to  train  the  FF  ANN  to  provide 
static  mapping  between  the  terminal  variables  and  the  developed  electromagnetic  torque.  The  desired 
performance  accuracy  of  the  mean  squared  error  of  [0.0005  x  (the  peak  expected  torque  value)^]  was 
obtained  using  a  FF  ANN  having  the  same  topology  as  that  one  in  the  6/4  SRM  case.  Fig.  4.9  shows  the 
performance  of  the  FF  ANN  for  the  8/6  SRM  prototype  in  the  torque  estimation  task.  The  single-phase 
flux-linkage  characteristics  at  different  positions  were  obtained  by  single-shot  rqrplications  of  a  dc  voltage 
source  to  the  phase. 

The  SRM  phase  terminal  voltages  and  currents,  recorded  using  a  Tektronix™  TDS  460  four-channel 
digitizing  oscilloscope,  were  used  to  conqrute  the  flux  linkages  as  function  of  current  For  conputing  the 
multi-phase  characteristics,  the  dc  voltage  is  simultaneously  applied  to  two  phases  (see  Figure  4.7).  The 
voltage  sources  can  be  turned  on  simultaneously  or  they  can  be  turned  on  one  at  a  time.  One  of  the  voltage 
sources  is  set  to  get  the  desired  current  level  in  one  of  the  phases.  Then,  the  other  voltage  source  is  turned 
on  so  that  the  current  in  that  phase  increases  fiom  zero  to  the  peak  phase  current.  Even  though  the  current  in 
the  first  phase  will  vary  slightly  due  to  the  mutual  interactions,  it  is  acceptable  since  the  idea  is  to  sanple 
the  space  of  the  function  M  to  use  these  values  for  training  of  the  ANN.  The  terminal  voltage  and  currents 
of  the  two  phases  are  recorded  and  used  to  compute  the  phase  fluxes  (see  Rgure  4.10).  Since  the  single¬ 
phase  characteristics  are  already  known,  we  can  then  compute  the  mutual  flux  and  M  value.  For  measuring 
M,  we  need  only  5  current  levels  in  one  phase  (say  leading  phase)  at  6  rotor  positions.  The  current  is  set  up 
in  the  other  phase  (say  the  working  phase)  by  switching  on  the  dc  voltage  source  whose  voltage  is  selected 
so  as  to  set  up  the  peak  phase  current  based  on  knowing  the  phase  resistance  value.  In  this  case,  each 
switching  of  the  dc  voltage  source  at  the  working  phase  will  provide  data  to  measure  the  M  function  values 
along  trajectories  parallel  to  the  working  phase  current  axis  on  surfaces  such  as  that  in  Figure  4.3.  Thus,  we 
need  5  x  6  =  30  measurements  to  characterize  M  completely. 
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In  Figure  4. 10,  the  current  in  (leading)  phase  2  was  set  at  a  fixed  level.  Then,  a  voltage  was  applied  to 
the  (working)  phase  1.  As  a  result,  the  current  in  phase  1  increases  fixim  zero  and  correspondingly,  the 
current  in  phase  2  undergoes  a  change  due  to  the  mutual  flux  interactions.  The  voltage  source  applied  to 
phase  2  changes  to  maintain  the  current  at  the  set  current  limit.  Such  waveforms  were  recorded  at  each 
measurement  point  and  used  to  calculate  the  phase  flux  linkages  under  these  conditions.  Since  the  self-flux 
linkages  have  been  already  recorded,  we  can  use  (4.8)  to  compute  the  M  function. 

These  computed  values  at  different  current  levels  and  rotor  positions  are  then  used  to  train  the  FF  ANN 
for  mapping  the  terminal  variables  to  the  M  function  values.  The  desired  performance  accuracy  was 
obtained  with  a  FF  ANN  having  an  input  layer  of  3  neurons,  two  hidden  layers  of  respectively  5  and  4 
neurons,  and  an  output  layer  of  1  neuron.  Figure  4.11  illustrates  the  performance  of  the  FF  ANN  in 
mapping  the  SRM  terminal  variables  to  the  measured  mutual  flux  interaction  functions.  It  is  noted  again 
that  M  has  negative  and  positive  values  meaning  that  the  effect  may  be  to  either  decrease  or  increase  the 
phase  flux.  As  anticipated  from  Figure  4.3,  the  mutual  interaction  functions  cannot  be  “fitted”  by  a  single 
mathematical  expression  (e.g.,  Fourier  series).  Therefore,  the  ANN  provides  better  performance  than  the 
other  methods.  Table  4. 1  summarizes  the  number  of  measured  pomts  in  the  ANN-approach  as  conqrared  to 
the  conventional  techniques  for  the  two-phase  excitation  case.  With  higher  number  of  excited  phases,  the 
number  of  measirrement  and  storage  points  required  with  the  conventional  methods  for  torque  mapping 
would  scale  up  by  one  more  order  (e.g.,  10000  points  of  measurement  and  storage  would  be  required  in  the 
three-phase  excitation  case).  In  the  case  of  the  mutual  interaction  function,  the  number  of  points  would 
scale  to  similar  values  since  a  Fourier  series  approximation  to  interpolate  between  the  function  values  at  4 
positions  would  no  longer  be  applicable. 

4.4  Simulations  And  Results 

The  described  multi-phase  model  was  applied  in  simulations  to  study  the  dynamic  performance  of  the 
SRM  drive  since  a  fully  operational  SRM  drive  was  not  available  in  the  laboratory.  The  full  model 
Including  the  torque  and  mutual  flux  interaction  function  ANN  was  inqrlemented  in  SIMULINK™. 
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Table  4.1  Summary  of  Conventional  and  ANN-approach 

Conventional  (interpolation)  ANN 


Mutual 

Interaction 

Function 

6  points  for  single  phase  plus 

40  points  with  two  phases.  Can  use 
spline  function  plus  Fourier  Series. 
400  points  storage. 

6  points  for  single  phase  plus  30 
points  for  two  phases. 

Torque 

About  1000  points  of  measurements 
(or  computation)  and  storage.  Must 
use  spline  function. 

80  points  for  single  phase  plus  225 
points  for  two  phases. 

-3 

x10 


Rotor  position  with  respect  to  working  phase  (degrees) 

Figure  4.11.  Performance  of  the  FF  ANN  for  the  mutual  flux  interaction  function 
mapping  in  the  four-phase  8/6  SRM  case. 

See  Appendix  D  for  the  full  SIMULINK™  model.  To  explore  one  of  the  possible  applications  for  this 
new  model,  a  simulation  study  was  also  implemented  for  a  position  sensorless  scheme  correcting  for  the 
mutual  flux  linkages. 

Figure  4.12  shows  the  simulation  results  for  the  torque  waveform  using  the  FF  ANN  based  on  the 
measured  torque  data  from  the  three-phase  6/4  SRM  at  a  speed  of  300  rpm.  The  torque  waveforms  clearly 
illustrate  the  difference  between  the  torque  obtained  using  a  multi-phase  model  and  the  principle  of 
superposition.  Measured  torque  are  not  shown  in  the  figure  since  the  torque  FF  ANN  was  tramed  using 
measured  results  for  a  mean  squared  error  value  of  [0.0005  x  (the  peak  torque  value)^]. 
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The  actual  average  torque  developed  by  the  SRM  is  lower  by  about  14%  when  compared  to  the 
“superposition”  results.  Thus,  the  “ANN”  torque  provides  a  far  better  estimation  of  the  actual  SRM  torque. 
The  degree  of  multi-phase  operation  is  also  illustrated  in  Figure  4.12.  It  is  observed  that  the  6/4  SRM 
operates  with  two  phases  excited  simultaneously  for  certain  periods  and  with  only  one  phase  excited  for 
other  periods. 

Figure  4.12  also  shows  the  simulation  results  for  the  position  error  in  the  position  estimation  system 
based  on  the  multi-phase  flux  linkage  model  at  300  rpm.  These  simulations  are  based  on  data  obtained 
from  FEA  simulations.  If  the  phase  flux  is  not  corrected  for  the  mutual  flux  linkage,  the  estimated  position 
may  deviate  from  the  actual  value  by  about  10°  leading  to  improper  SRM  operation. 
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Figure  4.12  Simulation  results  in  the  multi-phase  mode  using  superposition  and  the 
oroDosed  ANN-based  model. 


Chapter  Five 

ADAPTIVE  ON-LINE  TORQUE  ESTIMATION 

5.1  Introduction 

For  SRM,  the  electromagnetic  torque  has  been  normally  calculated  using  the  flux-linkage 
characteristics  of  the  motor  to  compute  the  co-energy  (as  a  function  of  rotor  position  at  a  constant 
current)  and  then  applying  the  classical  equation  of  torque  as  a  function  of  co-energy.  However, 
the  actual  flux-linkage  characteristics  for  a  particular  SRM  may  be  different  from  that  of  the 
designed  motor  or  the  initial  prototype  for  which  the  electromagnetic  torque  computations  have 
been  made  off-line.  This  is  because  of  manufacturing  tolerances  and  the  variance  in  the  behavior 
of  magnetic  materials  among  other  factors. 

Any  approach  that  computes  the  torque  based  on  standstill  measurements  of  inductance 
performed  on  one  prototype  [1-2]  is  not  adequate  when  accurate  torque  feedback  is  desired  for 
robust  dynamic  control.  The  problem  is  further  compounded  by  the  dependence  of  torque  not 
only  on  the  current  level  but  also  on  the  rotor  position.  Extensive  off-line  measurements  and 
calculations  would  be  required  to  obtain  the  electromagnetic  torque  characteristics  as  function  of 
current  and  rotor  position.  The  torque  data  to  be  stored  would  be  extremely  large  and  the 
algorithm  retrieving  the  torque  from  those  values  stored  at  discrete  current,  rotor  position  and 
speed  values  would  have  to  carry  out  extensive  on-line  interpolation  over  non-linear  surfaces. 
Moreover,  these  techniques  would  fail  with  slight  changes  in  the  “off-line”  characteristics  of  the 
machine  with  aging  (e.g.,  changes  in  the  air-gap  or  in  the  magnetic  material  properties).  An 
improved  technique  was  presented  in  [3]  to  estimate  the  flux-linkage  characteristic  of  the  SRM 
using  on-line  measurements.  In  [1-3],  the  electromagnetic  torque  developed  is  indirectly 
computed  through  complex  co-energy  calculations  by  integrating  the  flux-linkage  characteristics 
up  to  the  current  level  at  which  the  torque  is  desired  between  two  different  rotor  positions. 
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This  technique  would  also  involve  surface  integrals  to  compute  the  flux  linkage  and  co¬ 
energy  to  get  the  actual  torque  produced  in  the  multi-phase  case  [4]. 

In  the  multi-phase  case,  results  obtained  by  superposition  of  the  single-phase  torque  for  each 
conducting  phases  is  unacceptable  as  shown  in  Chapter  4  [5].  Thus,  techniques  such  as  those  in 
[1-3]  would  be  computationally  expensive,  especially  for  the  multiply  excited  SRM  case. 

ANN-based  models  providing  static  mapping  between  the  SRM  terminal  variables  and  the 
electromagnetic  torque  have  been  proposed  and  applied  in  [5-6].  To  overcome  the  limitations  of 
other  techniques,  this  chapter  presents  a  new  and  simple  technique  for  estimating  the  SRM  torque 
in  the  multi-phase  case  using  an  ANN  based  on  the  SRM  terminal  voltage  equation.  This 
technique  is  simpler  to  apply  and  can  adapt  to  changes  in  the  SRM  characteristics.  The  proposed 
technique  is  intended  to  allow  the  ANN  SRM  model  suggested  in  [5]  to  adapt  and  correct  its 
training  based  on  the  current  characteristics  of  the  individual  SRM.  In  the  proposed  technique,  an 
ANN  estimator  is  applied  to  learn  the  characteristics  of  the  machine  using  on-line  measurements 
of  the  phase  currents  and  voltages.  The  technique  does  not  require  any  prior  off-line  simulations 
or  computations;  hence,  it  is  simpler  to  implement  than  other  methods.  The  system  is  used  to 
estimate  the  speed  component  of  the  net  applied  terminal  voltage.  This  component  is  then  used  to 
directly  calculate  the  electromagnetic  torque  developed  in  the  SRM. 

An  approach  for  applying  ANN  for  on-line  torque  estimation  has  been  proposed  previously 
[7].  This  approach  requires  knowledge  of  the  torque  in  the  previous  time  instants.  When  such 
estimation  operates  in  the  absence  of  a  torque  transducer  (the  main  goal  of  the  estimation  is  to 
eliminate  the  torque  transducer),  the  results  obtained  are  inaccurate.  The  technique  presented  here 
requires  no  prior  knowledge  of  the  torque  or  the  system  characteristics.  Figure  5.1  shows  the  role 
of  torque  feedback  in  the  control  schematic  of  a  SRM.  The  entire  estimation  and  torque  mapping 
is  proposed  to  be  implemented  on  a  fast  DSP  such  as  the  C6701  expressDSP™  from  Texas 
Instruments.  This  implementation  will  be  described  in  more  detail  in  the  Chapter  Six. 
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5.2  Novel  Torque  Estimation  Technique 

The  terminal  voltage  equation  of  the  SRM  is  given  by: 


The  back  emf  e  is  then  given  by  : 


D-  dX 
V  =  Ri  + - 

dt 


dX']  (dX'\  f di']  fdX^fdd 


-  +  —  —  =e«  + 


The  ANN  topology  could  be  chosen  to  identify  the  net  back  emf  e  given  by  (5.1).  The  inputs 
and  output  of  the  ANN  are  chosen  as  samples  at  equally  spaced  time  instants  like  corresponding 

to  a  practical  digital  implementation.  To  represent  ,  we  use  (0*-  Ot-i)  and  for  ,  we  use 


(ir  ik-i).  Since  the  sampling  time  T  is  constant,  the  ANN  weights  adjust  to  take  care  of  this  value, 
allowing  the  use  of  difference  values.  In  the  first  topology,  the  speed  component  of  the  back  emf 
would  be  determined  by  setting  the  current  change  (4-  4.y)  to  zero  (see  Figures  5.1  and  5.2). 

However,  this  choice  of  ANN  topology  results  in  the  disadvantage  that  the  ANN  output  Ck  is 
presented  with  binary  values  (either  or  0)  depending  on  whether  the  power  converter 


70 


switches  are  ON  or  OFF.  Thus,  the  ANN  does  not  have  adequate  output  patterns  to  be  able  to 
generalize  the  output  for  other  cases.  As  a  result,  this  topology  does  not  perform  well.  To 


overcome  this  limitation,  we  use  the  observation  that  the  value  of 


is  rich  in  patterns  that 


vary  with  the  applied  voltage,  speed,  rotor  position  and  current  level.  Hence,  we  choose  a  second 
ANN  topology  that  has  the  form  of  Figure  5.3.  In  this  case  the  ANN  is  trained  to  identify  the 
change  in  the  phase  current  (4-  4-/)-  The  ANN  emulates  the  SRM  terminal  equation  re-written  as; 


raA) 

1  (40) 

L  dt 

L  _  ''jt  ^(Ok 

(5.2) 


di 


Ml 

3/ J* 


To  determine  the  speed  component  of  the  back  emf  we  need  to  determine  the  back  emf 


that  would  set  the 


to  zero.  To  do  this,  we  set  the  ANN  output  to  0  and  then  invert  the  ANN 


to  determine  the  corresponding  e*  (see  Figure  5.4).  This  value  is  the  speed-emf  component,  say 
gfl*.  Having  determined  the  speed-emf,  the  developed  electromagnetic  torque  at  the  k*  time 
instant  is  then  computed  using  the  following  equation; 

xi.xT 


y  _  ^(ok 


(5.3) 
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Figure  5.4.  ANN  inversion  for  determining  the  speed  emf  component. 

5.3  The  Inversion  Algorithm 

Several  techniques  have  been  proposed  in  the  literature  for  inverting  feed-forward  ANN  [8- 
9].  In  our  case,  the  task  is  simplified  because  three  of  the  four  input  variables  are  already  known 
at  the  k*  instant.  Furthermore,  we  can  use  the  property  of  electric  motors  that  the  speed  emf  has 
to  be  a  continuous  function  of  time. 

Based  on  this,  we  can  start  with  a  speed  emf  value  of  zero  and  in  each  case  search  for  the 
neighborhood  values  that  minimize  the  output  (/*-  4-/)  of  the  ANN  in  Fig.  5.4.  Thus,  we  set  values 
of  in  the  neighborhood  of  the  value  in  the  previous  time  instant  and  carry  out  forward 
propagation  to  search  for  that  value  of  Coik  minimizing  the  output  of  the  ANN.  The  algorithm  for 
the  estimation  system  was  implemented  in  MATLAB™  and  estimates  the  electromagnetic  torque 
developed  in  a  four-phase  8/6  SRM  drive  system  modeled  in  SIMULINK™  (see  Appendix  D) 
using  on-line  measurements  of  voltage,  current  and  rotor  position. 
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Figure  5.2.  ANN  topology  for  back  emf  identification  (see  (5.1)). 


Figure  5.3  ANN  topology  for  current  change  identification  (see  (5.2)). 
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The  SIMULINK™  SRM  model  was  used  to  produce  the  “on-line”  measurements  because  a 
fully  operational  SRM  drive  is  not  available  in  the  lab. 

Fig.  5.5  shows  the  flow  chart  of  the  algorithm  for  the  electromagnetic  torque  computation.  At 
the  start  of  the  algorithm,  the  speed  emf  estimate  is  set  to  zero.  At  every  time  instant  k,  the 
position  Qk-i  and  current  4-/  from  the  previous  time  instant  are  retrieved  from  the  SIMULINK™ 
SRM  model  and  used  to  compute  the  position  and  current  changes.  The  position  6k  and  current  4 
samples  then  overwrite  the  values  stored  in  the  previous  time  instant.  Based  on  these  “measured” 
and  computed  values,  the  ANN  in  Figure  5.3  is  trained  to  identify  the  change  in  current  (4  -  4-/4 
The  speed  emf  is  then  obtained  by  inverting  the  ANN  using  an  iterative  search  technique 
assuming  a  range  for  Ca*  as  shown  in  the  flow  chart.  Remember  that  e^k  cannot  experience  a  large 
variation  instantaneously,  so  we  assume  that  this  variation  is  within  ±5  %  of  the  DC  supply  V*- 


Figure  5.5.  The  flow  chart  for  the  electromagnetic  torque  computation 
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5.4  Simulation  and  Implementation  Results 

The  ANN  topology  for  the  estimation  task  was  optimized  using  MATLAB™.  The  ANN  has 
an  input  layer  of  4  neurons,  a  hidden  layer  with  3  neurons  using  tan-sigmoid  activation  functions 
and  an  output  layer  with  one  neuron  having  a  linear  activation  function. 

The  ANN  estimator  is  initiated  with  random  weights.  The  estimator  trains  to  predict  the 
change  in  current  (4  -  ik./),  for  given  back  emf  e*,  current  4,  rotor  position  0*  and  speed  cOk-  It  was 
found  that  speed  and  current  level  changes  in  the  SRM  are  essential  for  the  ANN  to  be  able  to 
train  correctly  and  assign  the  individual  components  of  the  back  emf  Ck  to  the  appropriate  causes 
(i.e.,  speed  and  current  changes).  The  data  corresponding  to  such  changes  can  be  easily  obtained 
since  speed  and  current  changes  are  inherently  present  in  a  practical  system  during  transients  due 
to  load  changes  and  startup. 

Fig.  5.6  shows  the  initial  response  of  the  estimation  system  after  startup  with  random  weights 
and  biases.  The  transients  during  changes  in  speed  and  current  levels  are  evident.  During  the 
simulations,  the  speed  and  currents  were  perturbed  by  step  changes  after  fixed  periods.  This  was 
done  to  allow  the  ANN  to  identify  the  characteristics  of  the  system. 

Fig.  5.7  shows  the  plot  of  the  estimated  electromagnetic  torque  along  with  that  obtained  from 
FEA  data  corresponding  to  the  SRM  current  waveform  at  1000  rpm  and  13  A.  The  estimation 
system  that  has  been  running  for  1  sec  (i.e.,  the  time  offset)  is  able  to  estimate  the  torque  to  a  very 
good  accuracy.  By  time  offset,  we  mean  that  the  ANN  started  with  no  prior  knowledge  of  the 
system  and  takes  1  sec  to  learn  the  system  characteristies.  Remember  that  an  advantage  of  the 
technique  proposed  here  is  that  no  off-line  measurements  are  required  as  in  other  techniques. 
Figure  5.8  shows  the  back  emf  and  its  speed  component  as  estimated  by  the  ANN  estimator. 
Preliminary  simulation  studies  for  the  multi-phase  case  indicate  that  the  influence  of  the  other 
phases  maybe  added  by  considering  additional  inputs  to  the  ANN.  In  this  study,  we  have  not 
considered  them. 
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FEA  computed  torque 


Figure  5.6.  The  ANN  estimator  response  following  startup  with  random  weights  and  biases. 


Time  (secs) 


Figure  5.7.  The  comparison  of  the  electromagnetic  torque  computed  using  the  ANN 
estimator  (solid)  and  using  FEA  (dashed)  at  1000  rpm  and  13  A.  (Time  offset  =  1  sec). 
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Figure  5.8.  The  measured  back  emf  e  and  the  estimated  speed  component  1000 
rpm  and  13  A  (Time  offset  =  1  sec). 
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Chapter  Six 

DSP  IMPLEMENTATION  OF  ANN  BASED  CONTROL  SOLUTIONS 

6.1  Introduction 

Previous  chapters  have  illustrated  that  ANN-based  solutions  can  be  applied  to  simplify  the 
control  tasks  in  SRM.  The  main  practical  problem  with  applying  ANN  till  now  has  been  that  their 
execution  speeds  are  limited  by  the  capability  of  the  available  processors.  However,  advanced 
architectures  such  as  the  expressDSP™  [1-2]  from  Texas  Instruments  (TI)  with  multiple  levels  of 
pipelining  and  higher  clock  speeds  are  now  available.  In  addition,  excellent  code  development 
tools  such  as  the  Code  Composer  Studio  (CCS)  [3]  support  these.  The  chapter  reports  on  the 
implementation  of  the  ANN-based  mapping  solutions  entirely  in  ANSI  C  using  the  features  of 
CCS  Version  1.2.  The  profiling  techniques  available  in  the  CCS  to  aid  the  rapid  code 
development  are  then  applied  to  measure  the  execution  times  of  the  different  ANN  stages.  The 
SRM  input  variables  obtained  through  simulations  at  different  operating  conditions  (in  lieu  of 
actual  measurements  being  available)  are  then  used  to  test  the  ANN  operation  and  response.  The 
chapter  also  gives  an  outline  of  the  schematic  for  interfacing  the  C6701  EVM  with  a  C240  EVM 
that  takes  care  of  the  low-level  tasks  of  generating  the  PWM  signals  for  SRM  drive  control. 
Finally,  the  chapter  concludes  by  describing  the  implementation  of  the  ANN  for  torque  mapping 
(see  Chapter  Four)  as  well  as  the  obtained  performance  results. 

6.2  Description  of  the  ANN  Implementation 

From  the  ANN  structure  presented  in  previous  chapters,  one  can  anticipate  that  an  important 
portion  of  the  ANN  computational  time  is  due  to  the  calculation  of  the  weighted  sums  at  each 
neuron  of  the  hidden  and  output  layers.  The  ANN  were  implemented  on  the  TMS320C6701  EVM 
entirely  in  software  using  the  C  language  and  the  features  available  in  the  CCS  Ver  1.2  [3j. 
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The  first  step  of  the  code  is  to  compute  the  weighted  sum  at  each  neuron  of  the  hidden  layer 
for  the  activation  functions.  The  C  language  code  implementing  these  computations  has  a  nested 
loop  structure.  The  bias  at  each  neuron  is  then  added  to  this  sum  to  obtain  the  net  activation  input. 
This  is  then  propagated  through  the  activation  function  at  the  neuron.  The  code  in  our  first 
implementation  had  the  following  form: 

void  comptoutO 

I 

int  i,j; 

for  (j=Oy<nl;j++)  netinllj]=biasllj]; 
for  (j=0;j<nl;j++) 

{ 

for  (i=0;i<n0;i++)  netinl[j]+=weightl|j][i]*netinO[i]; 
actvl  |j]=squash(netinl  [j]); 

} 

netin2[0]=bias2[0]; 

for  (i=0;i<nl;i++)  netin2[0]+=weight2[0][i]*actvl[i]; 

netout=netin2[0]; 

1 


where  float  squash(float)  computes  the  log-sigmoid  function  using  the  functions  from  the  math  library. 

Another  time-consuming  activity  is  the  computation  of  the  activation  function  using  the 
following  log-sigmoid  expression: 


fix)  = 


(6.1) 


One  can  use  the  function  “expO”  or  “expf()”  to  implement  this  equation  in  the  C  program. 
The  CCS  profiler  was  used  to  measure  the  number  of  cycles  required  by  each  function  and  the 
precision  of  the  value  returned  (by  the  exp()  and  expf()  functions).  Based  on  this,  it  was 
determined  that  the  expf()  function  requires  much  lesser  number  of  cycles  (390)  than  the  exp() 
function  (1030)  and  also  offers  adequate  precision  for  the  application.  The  log-sigmoid  function 
value,  thus  computed,  has  a  value  of  0  for  less  or  equal  than  -15  and  a  value  of  1.0  for  greater  or 
equal  than  15.  As  a  result,  it  is  necessary  to  compute  this  function  only  in  the  range  [-15  <  x  < 
15]. 


With  the  log-sigmoid  computation  using  the  expf()  function,  the  forward  propagation 
(computing  the  torque  output)  takes  about  3900  cycles  in  the  worst  case  when  [-15  <  x  <  15]. 
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This  is  because  the  usage  of  the  expf()  function  leads  to  the  disadvantage  that  the  code 
contains  function  calls  that  cannot  be  software  pipelined. 

To  be  able  to  use  pipelined  code,  we  use  the  “lookup-table”  technique  for  determining  the 
log-sigmoid  function  values  for  inputs  between  the  range  -15  to  15;  these  computed  values  were 
stored  in  a  lookup  table.  For  a  given  input  value  for  which  the  function  value  is  required,  the 
index  into  the  table  is  calculated  using: 

(6.2) 

I  30  J 

where  sig Points  is  the  number  of  points  in  the  lookup  table  for  the  log-sigmoid  function.  This 
equation  was  implemented  in  the  C  program  using  type-casting  to  integer  value.  However,  it  was 
found  that  direct  implementation  of  this  equation  in  the  C  program  leads  to  the  requirement  of  a 
large  number  of  cycles  (about  2800  cycles  in  the  worst  case  with  [—15  <  x  <  15])  because  the 
division  process  which  is  used  in  the  index  computation  is  very  time  consuming.  This  was 
overcome  by  replacing  the  fraction  sigPoints/30  by  its  value  (say  q)  in  the  program;  that  is 

m  =  round  (6.3) 

The  initial  code  using  this  formulation  required  about  540  cycles  for  the  forward  propagation 
because  the  code  uses  simple  addition  and  multiplication  that  are  the  “least  computationally 
expensive”  for  the  processing.  Also,  the  compiler  was  able  to  setup  better  pipelining  since  the 
number  of  registers  required  during  each  iteration  was  small.  The  loop  structure  and  the 
arrangement  of  the  code  was  further  improved  using  the  CCS  Ver  1.2  features  [8-11]  and  the  final 
code  obtained  required  only  217  cycles  for  the  torque  computation  by  forward  propagation. 
These  217  cycles  at  a  clock  of  133  MHz  mean  that  the  C6701  requires  only 

=  1  63  u  sec  for  the  torque  computation  by  forward  propagation. 

133x10^ 
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The  following  features  of  the  CCS  were  used  for  improving  the  code  arrangement  and 


pipelining  [3]: 

•  The  profiler  was  used  to  determine  the  number  of  cycles  taken  by  each  section  of  the 
code  and  thus  identify  the  “expensive”  sections. 

•  The  #pragma  MUST_ITERATE  directive  was  used  to  pass  the  information  about  the  trip 
count  (number  of  loop  iterations)  to  the  compiler. 

•  The  #pragma  UNROLL  directive  was  used  to  unroll  small  loops  and  increase  the  number 
of  instructions  available  for  execution  in  the  pipeline. 

•  The  -pm  and  -mt  compiler  options  were  respectively  used  to  direct  the  compiler  to  use 
the  trip  count  data  for  pipelining,  and  to  indicate  the  absence  of  memory  aliasing  [4]. 

•  The  -k  option  preserved  the  .asm  file  for  inspection.  The  compiler  includes  feedback 
about  the  pipelining  in  this  file  and  also  indicates  better  options  that  the  user  might 
choose.  The  code  was  rearranged  till  smaller  iteration  intervals  and  a  larger  number  of 
parallel  iterations  were  obtained. 

The  final  code  for  the  forward  propagation  has  the  following  form: 

void  comptoutO 
{ 

intj; 
int  m; 

#pragma  MUSTJTER  ATE(n  1 ); 
for  (j=0;j<nl;j++) 

1 

netinl[j]=biasl|J]; 

netin  1  [j  ]+=weight  1  [j  1  [0]  *netin0  [0]+weight  1  [j  ]  [  1 1 
netin0[  1  ]+weight  1  [j  ]  [2]netin0[2] ; 

1 

#pragma  MUST_ITERATE(nl); 
for(j=0y<nl;j++) 

{ 

if  ((netin  1  [j]>- 1 5)&&(netin  1  |jl<l  5)) 

1 

m=(int)((netinl  [jl+l  5.)*80); 
actvl  |J]=sigTable[m]; 

1 

if  (netinl[j]<-15)  { actvl |J]=0.0;} 
if  (netin  1  |j]>  1 5)  { actv  1  [)]= 1 .0; } 

1 

netin2[0]=bias2[0]; 

#pragma  UNROLL(nl); 

for  (i=0;j<nl;j++)  netin2[01+=weight2[0]|j]*actvl|J]; 
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netout=netin2[0]; 


} 

where  sigTablet]  contains  the  value  for  the  log-sigmoid  function  in  the  range  -15  to  15  and  has  2400  points. 

The  process  of  learning  by  error  back-propagation  involves  the  computation  of  activations  at 
the  neurons.  The  back-propagation  algorithm  was  also  implemented  in  C  and  the  code  was  then 
optimized.  The  resulting  code  executed  in  344  cycles  or  2.58  ns  per  iteration. 

For  applications  involving  electric  drives,  torque  feedback  at  a  sampling  rate  of  about  25  kHz 
maybe  considered  to  be  adequate  under  most  conditions.  Since  the  C6701  is  able  to  compute  the 
torque  output  under  multi-phase  operation  at  a  much  higher  speed,  there  is  sufficient  time 
available  between  the  samples  to  schedule  training  of  an  adaptive  ANN-based  model  (which 
changes  with  the  actual  motor  parameters)  based  on  on-line  measurements  (see  Chapter  Five). 
This  can  be  used  to  schedule  calculations  to  compute  the  new  actual  values  of  torque.  Finally,  the 
error  between  the  ANN  torque  output  and  the  new  calculated  value  can  be  used  to  schedule 
training  of  the  torque  ANN.  Thus,  the  C6701  can  be  used  to  implement  several  ANN  that  cannot 
only  estimate  the  drive  output,  but  also  adapt  to  changing  drive  parameters. 

6.3  Testing  of  the  ANN  using  Probe  Points 

The  CCS  has  the  facility  to  connect  data  from  files  on  the  PC  to  the  program  running  on  the 
target  DSP.  This  is  done  using  probe  points  [3].  The  SRM  current  and  rotor  position  data  were 
obtained  from  SIMULINK™  simulations  [5]  and  the  data  corresponding  to  two  different  speeds, 
1(X)  rpm  and  6000  rpm,  were  stored  in  ASCII  text  files  in  floating-point  format.  In  a  real 
implementation,  real  measurements  would  take  the  place  of  these  ASCII  files.  Probe  points  were 
then  used  to  couple  these  data  files  at  runtime  to  the  C6701.  This  was  used  to  test  the  ANN 
output. 
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Portions  of  the  resulting  graphs  for  the  two  winding  currents,  the  rotor  position  and  the  output 
torque  are  included  as  Figures  6.1  and  6.2,  which  are  screen  snapshots  taken  during  the 
simulations. 

With  sigPoints  equal  to  2400,  it  was  found  that  the  resulting  torque  waveform  includes 
additional  high-frequency  noise  introduced  by  the  ANN  due  to  rounding  of  the  index  value. 
Increasing  sigPoints  to  6000  mitigated  this  noise. 

6.4  Interfacing  the  C6701 EVM  to  the  C240  EVM 

The  intended  role  of  the  C6701  in  the  SRM  drive  system  is  to  provide  computational  support 
to  the  drive  controller.  The  main  reason  is  that  present-day  DSPs  for  motion  control  do  not  have 
the  required  bandwidth  to  implement  novel  control  algorithms  (e.g.,  model  reference  adaptive 
control)  when  the  electric  motor  is  running  at  high  speeds  (e.g.,  greater  than  6000  rpm).  In  our 
implementation,  the  SRM  will  be  controlled,  at  the  low-level,  through  PWM  signals  generated 
using  the  C240  EVM.  An  interface  between  the  two  DSPs  has  to  be  designed  to  transfer  the 
computed  torque  value  to  the  C240  EVM  as  feedback.  The  interface  between  the  C6701  and  the 
C240  EVM  can  be  implemented  conveniently  through  the  C6701  Host  Port  Interface  (HPI)  [6]. 
The  HPI  is  a  I6-bit-wide  parallel  port  through  which  the  C6701  memory  space  is  visible  to  the 
host  processor  (the  C240,  in  this  case).  The  C240  functions  as  the  master  to  the  interface.  The 
data  exchange  can  take  place  using  internal  or  external  memory.  The  C240  can  also  access  the 
memory-mapped  peripherals.  Connectivity  to  the  C670rs  memory  space  is  provided  through  the 
DMA  controller  in  the  CPU. 

The  HRDY  pin  can  be  used  to  insert  a  delay  during  a  transfer  and  characterize  the  data 
transfer  speed  through  the  HPI.  Since  the  C240  has  a  separate  address  bus,  the  HAS  can  be 
inactive  high  at  all  times. 
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Figure  6.2.  Test  results  of  the  ANN  processing  using  probe  points  in  CCS  Ver  1.2  at  6000  rpm.  The 
graphs  indicate  (CCW  from  the  top  left  corner)  the  rotor  position  (with  respect  to  the  working  phase),  the 
ANN  torque  output,  the  current  in  the  leading  phase,  and  the  current  in  the  working  phase. 
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HCNTL  [1:0]  values  would  be  set  by  decoding  the  address  of  the  C240.  The  base  address  is 
determined  by  the  decoding  logic  and  is  chosen  such  that  HCNTLl  =  0,  HCNTLO  =  0  and 
HHWIL  =  0. 

On  the  C240  EVM,  the  RAMOE  on  the  GAL  device  controls  the  read  from  external  memory 
operation.  On  the  other  hand,  RAMWE  controls  the  write  to  memory  operation  [7].  These 
signals  are  used  as  separate  strobes  for  the  hds  i  and  hds  2  for  read  and  write.  The  C240 
R/W  signal  is  connected  to  the  HPI  signal  input.  HWOB=l  indicates  that  the  first  write  is 
the  least  significant  half  word.  The  connection  diagram  for  the  interface  between  the  C6701  and 
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Figure  6.3  Interface  between  the  C6701  and  C240  through  the  HPI. 
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the  C240  EVM  is  shown  in  Figure  6.3. 


HBE  [1:0]  =0  since  only  a  half  word  (16  bits)  is  written  by  the  C240  during  a  write 
operation.  The  HPI  Control  register  (HPIC)  is  normally  the  first  register  accessed  to  set  the 
configuration  bits  and  initialize  the  interface  [6].  The  DSPINT  bit  can  be  used  by  the  C240  to 
interrupt  the  C6701.  The  FETCH  bit  of  the  HPIC  is  used  with  the  hrdy  during  a  read  or  write  to 
set  up  a  software  handshake. 

The  base  address  of  the  C6701  HPI  is  chosen  during  design.  The  addresses  of  the  control, 
address  and  data  registers  of  the  HPI  are  relative  displacements  to  this  address  and  are  thus 
accessible  by  programs  written  on  the  host  C240.  For  example,  if: 

“hpiBaseAddr”  denotes  the  base  address  of  the  C6701  HPI  chosen  to  access  the  control 
register  first  half  word. 

Then:“hpiBaseAddr+r’  accesses  the  second  half  word  of  the  control  register  HPIC  to 
change  the  control  settings  and  initialize  the  interface. 

“hpiBaseAddr+2”  accesses  the  first  half  word  of  the  address  register  HPIA  to  write  the 
address  of  the  C6701  address  map  to  which  access  is  required. 

“hpiBaseAddr+3”  accesses  the  second  half  word  of  the  address  register  HPIA  to  write  the 
address  of  the  C6701  address  map  to  which  access  is  required. 

“hpiBaseAddr+4”  accesses  the  first  half  word  of  the  data  corresponding  to  the  address 
placed  in  the  address  register  through  HPID. 

“hpiBaseAddr+5”  accesses  the  second  half  word  of  the  data  corresponding  to  the  address 
placed  in  the  address  register  through  HPID. 
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6.5  Preliminary  Implementation  of  the  Adaptive  Estimation  ANN  and  DSP/BIOS 
A  first  implementation  of  the  Adaptive  Torque  Estimation  ANN  indicates  that  each 
training  by  back  propagation  would  require  less  than  2  psec  and  each  iterative  search  by  inversion 
would  require  less  than  5  psec. 

Since  the  ANN  for  static  torque  mapping  requires  about  1.63  psec  for  forward 
propagation  and  2.4  psec  for  training  by  back  propagation,  when  the  complete  system  is  in  a 
training  mode,  the  estimation  and  computation  tasks  would  require  less  than  12  psec.  The 
simulations  have  assumed  a  sampling  time  of  50  psec  (sampling  frequency  of  20  kHz). 

A  DSP/BIOS  implementation  was  also  done  to  test  the  execution  of  the  torque  ANN  thread  by 
simulating  on  the  DSP  conditions  similar  to  that  under  which  the  torque  ANN  would  execute  in 
the  final  implementation  (a  thread  is  a  schedulable  unit  of  execution  in  a  multitasking  system). 
The  results  are  displayed  in  Figure  6.4.  A  hardware  interrupt  (HWI)  driven  by  the  timer  clock  at 
intervals  of  25  psec  initiates  a  dummy  data  transfer  function. 

This  then  triggers  a  software  interrupt  that  executes  the  torque  ANN  (ANN_SWI).  A  load 
loop  of  1000  instructions  is  also  executed  with  the  torque  ANN  to  allow  adequate  margin  for 
additional  operations.  After  the  execution  of  ANN_SWI,  KNL_swi  is  initiated  to  invoke  the  task 
scheduler  for  executing  other  threads  (see  Figure  6.4). 
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Figure  6.4.  The  task  scheduling  in  the  DSP/BIOS  implementation;  the  interval  between  tics 

on  the  time  axis  is  25  psec. 
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Chapter  Seven 
CONCLUSIONS 


This  Part  1  of  the  Final  Report  has  provided  the  following  contributions  to  the  fields  of  Switched 

Reluctance  Motors  and  electric  vehicle  propulsion: 

•  Presented  an  approach  for  determining  the  SRM  power  ratings  for  a  complete  specification  in  an 
electric  propulsion  application.  A  comparison  of  the  thermal  duty  of  the  SRM  with  other  electric 
motors  was  also  presented.  The  PMSM  has  considerably  higher  thermal  duty  as  compared  to  the  other 
motors.  The  results  show  that  the  thermal  duty  of  the  SRM  can  be  reduced  considerably  by  appropriate 
choice  in  the  design  process  of  the  contribution  by  the  individual  loss  components  to  the  total  machine 
losses.  The  SRM  is  comparable  to  the  traditional  motors  under  these  conditions.  Future  work  is  the 
need  to  develop  a  procedure  to  design  a  SRM  so  that  the  losses  ate  limited  to  the  desired  proportion  of 
the  total  and  the  thermal  duty  is  kept  low. 

•  Developed  an  approach  for  designing  SRM  under  multi-phase  excitation  mode  while  taking 
the  practical  loading  constraints.  The  studies  leading  to  the  development  of  this  approach  also 
provided  insight  into  the  relative  merits  and  de-merits  of  each  of  the  suggested  and  popular 
SRM  configurations  so  as  to  further  aid  in  the  choice  of  an  appropriate  configuration.  The 
design  results  show  that  an  improvement  in  torque  density  is  obtained  in  those  configurations  that 
increase  the  number  of  pole  pairs  participating  in  the  torque  production.  To  this  end,  the  full-  and 
fractional-pitch  windings  do  not  change  the  number  of  pole  pairs  participating  in  torque  production. 
Also,  it  is  not  possible  to  improve  the  torque  density  substantially  by  using  full-  and  fi^ictional-pitch 
windings  due  to  the  thermal  constraints  inqxjsed  by  the  maximum  flux  density.  The  apparent 
advantage  of  a  better  utilization  of  the  winding  area  in  the  full-  and  fiactional-pitch  SRM  is  reduced  by 
the  higher  value  of  Ka  for  these  configurations.  The  short-pitch  SRM  use  unipolar  converters  and  the 
voltage  applied  across  the  windings  is  twice  of  that  for  the  full-  and  fractional-pitch  SRM  which  may 
use  bipolar  converters  and  start-connected  windings.  The  peak  current  for  the  converter  switches  thus 
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remains  unchanged.  The  choice  of  a  particular  SRM  configuration  should  be  also  based  on  other 
factors  such  as  high-speed  operation,  use  of  readily  available  standard  inverters,  ease  of  manufacturing, 
requirement  of  low  torque  ripple  using  current  profiling  between  phases,  noise-reduction  needs  and 
system  losses.  In  case  of  SRM  having  odd  number  of  phases,  continuous  operation  in  the  short-flux- 
path  mode  can  be  sustained  by  sinple  re-connection  of  phase  leads  to  get  the  appropriate  current 
directions.  Thus  in  the  case  of  the  1(V8  five-phase  SRM,  simple  uni-polar  converters  can  be  einiloyed 
while  still  getting  all  the  benefits  of  the  short-flux-path  mode  such  as  lower  acoustic  noise,  vibrations 
and  core  losses.  All  other  suggested  solutions  for  obtaining  this  mode  of  operation  in  the  10/8  SRM 
require  bipolar  converters.  Even  the  shared  bridge  configuration  suggested  in  reference  [1]  of  Chapter 
Three  will  not  lead  to  continuous  short-flux-path-mode  operation. 

•  Developed  a  simple  FF  ANN  model  for  the  multiply  excited  SRM  that  can  be  applied  in 
simulations  as  well  as  for  real-time  implementations.  The  proposed  techniques  require  small 
measured  data  sets  and  simple  ANN  topologies.  The  feasibility  of  implementing  these 
techniques  was  demonstrated  through  implementation  on  the  TI  TMS320C6701 
eXpressDSP™.  The  implementation  shows  that  these  techniques  can  be  applied  for  the  real¬ 
time  estimation  tasks,  while  leaving  adequate  room  for  other  on-line  adaptive  ANN  based 
solutions. 

•  A  simple  method  that  can  be  applied  for  direct  on-line  estimation  of  the  electromagnetic 
torque  developed  in  the  SRM  was  also  proposed.  The  technique  can  be  scheduled  in  a 
training  mode  during  commissioning  of  the  drive  or  occasionally  during  the  drive  life-time  to 
identify  the  current  electromagnetic  torque  characteristics  of  the  SRM  and  train  another  ANN 
that  would  then  provide  torque  feedback.  The  technique  is  simple  and  feasible  to  implement 
and  can  be  easily  extended  to  account  for  interactions  in  the  multiply  excited  case.  This 
method  does  not  work  well  at  high  speeds  (greater  than  rated  speed  of  the  SRM  drive)  due  to 
reduction  in  the  value  of  (4-  4./).  However,  this  method  can  be  applied  for  estimating  the  low 
speed  or  static  torque  characteristics  of  the  SRM. 
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•  It  may  be  possible  to  fit  the  measured  data  using  non-linear  ARMA  models  characterizing  the 
SRM  non-linear  characteristics.  However,  this  approach  was  not  adopted  since  it  was 
expected  that  it  might  be  more  computational  expensive  to  do  the  model  fitting  in  order  to 
allow  on-line  adaptation  as  compared  to  the  proposed  ANN  solutions.  Future  work  may 
investigate  the  feasibility  of  this  approach. 

•  The  current  models  for  torque  mapping  and  on-line  estimation  (including  those  of  this  report) 
do  not  consider  the  influence  of  losses  (due  to  eddy  currents  at  high  speeds)  on  the 
electromagnetic  torque  output  and  flux  at  a  given  phase  current  and  rotor  position.  It  is 
essential  to  subtract  the  loss  component  of  current  from  the  net  phase  current  to  get  the 
component  contributing  to  electromagnetic  torque  production.  This  requires  further  work  to 
introduce  models  for  the  loss  component  of  the  current  in  order  to  perform  this  correction. 

•  The  studies  of  the  report  indicated  the  feasibility  of  implementing  several  ANNs  to  aid  in 
tackling  different  aspects  of  the  SRM  control  problem.  A  scheme  for  interfacing  the  C6701 
EVM  and  the  C240  EVM  was  also  presented. 

•  At  the  present  time,  the  cost  of  the  C6701  may  make  it  unfeasible  for  low-power  low-cost 
drives.  However,  the  control  benefits  obtained  by  applying  ANN  aids  makes  it  attractive  for 
the  higher-power  drives.  The  feasibility  of  incorporating  high  speeds  ANNs  as  control  aids 
when  they  are  implemented  on  the  C6xxx  architecture  may  also  indicate  the  direction  for 
possible  enhancements  in  the  architecture  of  DSPs  such  as  the  C240,  that  are  specialized  for 
motion  control.  It  might  be  useful  to  incorporate  some  of  the  features  of  express-DSPs  such 
as  the  C6xxx  on  the  motion  control  DSPs  to  allow  the  implementation  of  more  advanced  and 
adaptive  control  techniques  based  on  ANNs.  The  report  has  shown  the  adequacy  of  present 
day  processors  to  implement  the  complete  ANN  based  system.  Further  work  is  required  to 
implement  and  apply  this  ANN  based  system  and  test  the  suggested  communication  scheme 
between  the  two  DSPs  . 
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Appendix  A 

DERIVATIONS  OF  SIZING  EQUATIONS 

A.l  Estimation  of  Power  Losses  in  a  SRM 

The  SRM  is  controlled  by  applying  a  current  pulse  to  a  phase  winding  whose  inductance  variation  as  a 
function  of  rotor  position  is  positive.  Assuming  an  ideal  rectangular  current  pulse,  Figure  A.1  shows  the 
current  waveform,  the  ideal  inductance  variation  of  the  torque  producmg  phase  and  the  stator-pole  flux 
waveform  for  two  different  current  levels.  The  flux  waveforms  in  other  regions  of  the  magnetic  circuit  are 
related  to  the  flux  waveform  in  the  stator  pole  [1].  The  current  pulse  either  occupies  the  full  feasible  region 
for  torque  production  (i.e.,  from  rotor  angle  0,  to  rotor  angle  0j)  or  the  dwell  angle  is  changed  so  that  it 
occupies  a  smaller  region.  This  is  illustrated  by  the  current  pulse  I  from  rotor  angle  0,  to  rotor  angle  0^  in 
Rgure  A.1.  For  the  same  average  torque,  the  amplitude  I  has  to  be  greater  than  that  of  I;  the  two  currents 
are  related  by -0,)  = /2(6), -0^). 

One  can  get  an  idea  of  the  meaning  of  the  difference  (6^-  0,)  by  noting  that  it  is  equal  to  30°  for  the 
three-phase  6/4  SRM.  So  the  number  of  phase  current  pulses  per  revolutions  (i.e.,  in  this  case,  4)  multiplied 
by  the  number  of  phases  (i.e.,  3)  and  [0^-  0,]  (i.e.,  30°)  is  equal  to  360°;  the  angle  corresponding  to  one  full 
revolution  of  the  rotor. 

The  relationships  relating  the  variations  of  the  power  losses  as  function  of  output  power  and 
motor  speed  are  necessary  in  order  to  determine  the  continuous  power  rating  of  the  SRM  in 
Chapter  Two.  Unfortunately,  it  is  not  possible  to  determine  the  exact  variations  unless  the  SRM 
is  already  designed  and  built;  therefore,  an  estimation  of  these  loss  variations  will  be  used  instead. 
The  following  assumptions  are  made  in  estimating  these  variations  of  power  losses  with  the  SRM 
output  power  and  speed: 

1)  Effects  of  magnetic  saturation  are  neglected. 

2)  The  formulation  relating  the  power  output  P^.  and  speed  to  the  losses  is  only  valid  as  long  as 
the  current  pulse  is  restricted  to  the  region  (02-0,)  or,  in  general,  to  the  torque  producing  zone 
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(see  Figure  A.l).  In  an  actual  case,  the  decaying  trailing  edge  of  the  current  pulse  will  be 
present  in  the  region  immediately  after  this  decaying-current  pulse  is  neglected. 

3)  It  is  assumed  that  the  current  pulse  applied  to  the  winding  always  has  a  rectangular-wave 
shape.  This  is  only  true  at  low-speed  operation  of  the  SRM.  For  high-speed  operation,  the 
applied  current  pulse  is  "peaky". 

Such  a  waveform  will  not  lead  to  any  additional  error  in  the  estimation  of  copper  losses 


Figure  A.1  (a)  The  ideal  phase  inductance  and  current,  (b)  The  stator-pole  flux  waveform.  The 
effect  of  reducing  dwell-angle  on  the  current  and  the  stator-pole  flux  is  illustrated  in  (b). 


94 


as  long  as  its  rms  value  is  identical  to  that  of  the  rectangular-wave  shape  case.  The  equations  for 
eddy-current  loss  and  hysteresis  loss  are  only  valid  as  long  as  the  current  waveform  is  flat-topped  (i.e., 
constant  amplitude).  The  error  in  the  estimated  losses  will  not  be  high  as  long  as  the  periods  when 
(di/dt)  is  large  are  small  [2]. 

4)  The  equation  for  the  copper  loss  is  not  valid  for  zero  motor  speed. 

5)  As  per  the  Steinmetz  equation,  the  hysteresis  loss  is  determined  by  (o ,  where  Steinmetz 

found  that  “n”  has  a  value  of  1.6  for  a  wide  variety  of  materials  [3].  For  the  new  ferromagnetic  alloys, 
“n”  may  vary  fiom  1.5  to  2.5.  In  fact,  “n”  is  not  a  constant  except  for  a  limited  range.  By  taking  “n” 
constant,  we  are  assuming  that  the  shape  of  the  hysteresis  curve  remains  unchanged.  This  assumption 
is  necessary  to  be  able  to  simplify  the  formulation  for  the  hysteresis  loss  in  the  absence  of  knowledge 
for  Therefore,  we  decided  to  use  the  hysteresis-loss  equation  W,,  co.  It  should  be  noted 

that  this  Steinmetz  equation  based  approach  has  been  applied  only  to  the  case  of  the  hysteresis  loss  as 
a  function  of  the  peak  flux-density  under  each  condition.  Since  the  complete  Steinmetz  equation  that 
also  accounts  for  eddy-current  losses  [  1]  is  applicable  only  to  the  case  of  sinusoidal  excitation,  we  have 
developed  an  alternate  approach  for  the  eddy-current  loss  based  on  the  basic  laws  governing  eddy- 
current  loss. 

From  the  driving  schedules  and  the  EV  main  specifications,  we  know  the  required  output  power  and 
speed  of  the  motor  drive  at  various  time  instants.  In  order  to  calculate  the  equivalent  continuous  power 
rating  of  the  SRM,  we  will  relate  the  losses  in  the  SRM  at  known  output  power  and  speed  (i.e.,  only 
one  operating  point)  to  the  required  output  powers  and  speeds  at  various  time  instants  in  the  considered 
driving  schedule  (i.e.,  P^  and  (O.).  The  following  analysis  illustrates  this  procedure. 


A.  1 . 1  Variation  of  the  copper  loss 
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Consider  a  SRM  producing  torque  T,  at  speed  to,  and  torque  at  speed  o^.  Since  the  SRM  torque 
varies  as  a  function  of  the  square  of  the  current  under  the  assumption  of  magnetic  linearity,  the  copper  loss 
is  then  directly  proportional  to  the  output  torque.  Thus,  the  copper  loss  at  these  two  points  are  related 


by: 

II. 

Since  7-  =  ,  the  copper  losses  at  two  different  output  powas  and  speeds  are  related  as : 

(0 

ly...  _  Pi  ^2 

H'.2  P2  0), 


(A.1) 


Now,  consider  that  the  torque  required  at  a  particular  speed  is  produced  by  applying  a  current  pulse 
throughout  the  angle  range  (6^-  dj  and  the  same  torque  is  produced  by  applying  a  current  pulse  throughout 
the  smaller  angle  range  (0^-6)  as  illustrated  in  Figure  A.  1 .  The  copper  loss  is  then  given  by: 

•  Current  pulse  throughout  the  angle  range  (B^  -  0,)  with  amplitude  I:  ,  oc  /  ^ . 

•  Current  pulse  with  reduced  dwell  angle  { 6  -  OJ  and  amplitude  I :  yy  «  /  "  ~^i)  oc  /  ^  • 

(02-0.) 

Thus,  the  length  of  the  dwell  angle  does  not  affect  the  copper  loss  for  a  given  torque  value. 

A.  1 .2  Variation  of  the  eddy-current  loss 

The  flux  in  the  stator  pole  is  proportional  to  the  phase  flux  linkages  that  may  be  expressed  as 

where  with  La  being  the  aligned  inductance  and  L„  the  unaligned  inductance. 

02—0) 


Also,  the  flux  density  B  is  proportional  to  the  flux  per  turn  and  also  to  the  flux  linkage  X.  Then,  the  rate 
of  change  of  the  flux  density  can  be  written  as  follows: 

dB  dX  d(LJ  +  KI  {6  -e.) 

-  OC  -  OC  - - - i - 

dt  dt  dt 
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dB 

dt 


oc  KIO) 


The  eddy-current  losses  can  be  approximated  by: 


W, 


W,  « 


Now,  assume  that  a  SRM  develops  torque  T,  at  speed  (o,  and  torque  Tj  at  speed  oij,  then  we  can  write: 


(O, 


=  L-^ILl 

p,  fflj 


(A.2) 


W  j  2  ^2  ®  I 

It  can  be  shown  that  the  eddy-current  loss  remains  unchanged  if  the  same  torque  at  a  particular 
speed  is  produced  by  either  applying  a  current  pulse  throughout  the  angle  range  (G^  -  0,)  or  by  using  a 
reduced  dwell-angle  since  W^oe  K¥a)^ 


A.  1.3  Variation  of  the  hysterisis  loss 

The  hysterisis  loss  is  given  by  Wj,  B]^(0  . 
Since®  oe0  oeA  =  (L  +  L  )I  • 
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If  the  same  torque  at  a  particular  speed  is  developed  by  reducing  the  dwell  angle,  then: 

=  LJ'+ 

then,  w,_  B_  ^  -b.  /■ 


W,  B. 


B,-e,  I 


w,  f 

- ii-  CSC 
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-e, 

Q  2  ~  ^  \ 


Above  rated  speed,  the  dwell-angle  is  changed  such  that 

_  Oir„,e^ 

02  -  0|  (O 
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Hence, 

W,  \  (0 

Based  on  the  above  formulations,  the  efficiency  contours  have  been  determined  and  plotted  in 
Figure  A.2  for  a  SRM  having  an  efficiency  of  85%  at  peak  output  power  and  rated  speed, /  = 
0.333  and  c  -  0.4.  These  contours  are  similar  in  shape  to  those  given  in  [4]. 


A.2  Derivation  of  the  Coefficients  Ki,K2e,  K2h 

Let  us  assume  that  we  know  the  copper  loss  the  hysteresis  loss  the  eddy-current 

loss  of  the  motor  at  the  peak  output  power  and  rated  speed  to 


Figure  A.2  Efficiency  contours  of  a  SRM. 
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A3  Derivation  of  the  Expression  for 

Knowing  the  values  of  K,,  and  at  a  particular  operating  condition  (e.g.,  a  driving 
schedule),  the  total  loss  to  be  dissipated  by  the  electric  motor  can  be  written  using  (A.4),  (A.5) 
and  (A.6)  as  follows; 
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WrouL  = 

w 

.  ^TOTM^  ^  a:,  x(l-c)  +  ATj,  xcx  /  +  /fj^xcxO-  /) 

'^SSMp 

The  losses  can  be  also  written  in  terms  of  the  continuous  power  rating  .  that  would  result 
in  the  same  losses  as  the  operating  condition  under  study  by  using  (A.1),  (A.2)  and  (A.3)  as 
follows: 


By  setting  ^ 


W, 


p 

^  COM  • 


TOTAL 


-Kp  +wj+ 


r  p 

^coM  ,i 


''  emp 


,..]?W.  =  .^{(l_c)  +  cX /}+ 
W  p  ^ 

^^SRMp  ‘emp 


\  -  J 
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^cont  jL(2A\  the  following  equation  results: 


•S|{(1  -  c)  +  c  X  /}+  ^/'Vx  (1  -  /)  = 
K^x(l-c)+  f  +  K^,^xcx(l-  f) 


(2.2) 
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Appendix  B 


B.l  C++  Program  for  Sizing  the  SRM 

This  program  was  developed  in  C++  to  take  advantage  of  the  Object  Oriented  Programming 
(OOP)  features  of  the  language.  The  EV  and  its  electric  motor  drive  are  represented  as  two 
classes  vehicle  and  drive,  respectively,  in  which  the  data  related  to  each  are  encapsulated.  A 
friend  function  interface  (vehicle, drive,float)  enables  calculation  of  the  power  ratings  of  the  drive 
by  implementing  an  interface  between  the  vehicle  and  the  drive  for  a  given  driving  schedule.  The 
advantage  of  such  an  OOP  based  code  is  its  easy  extensibility  to  multiple  drives  for  one  EV  or 
multiple  types  of  drives.  For  example,  this  code  could  also  be  easily  extended  for  sizing  the  SRM 
drive  for  a  Hybrid  EV,  The  program  is  listed  below: 
srmsiz,c 

/*In  this  program,  the  drive  power  requirement  is  calculated  based  on  the  thought  that  for  the  velocity  samples 
available,  the  drive  power  calculation  will  be  at  a  instant  exactly  at  the  center  of  the  period  between  the  two  instant 
where  the  velocity  sample  is  available.  Assuming  a  linear  velocity  profile  between  the  two  instants,  the  force  from  the 
road  load  characteristic  can  be  calculated  corresponding  to  the  velocity  at  this  center  instant  and  the  force  for 
acceleration  is  the  function  of  the  end  velocities  and  the  time-difference*/ 

#include  <stdio.h> 

#include  <iostream.h> 

#include  <math.h> 

//The  parameters  defined  here  are  hard-coded  and  need  to  be  changed  for  each  case 

#define  kd  0.29  //dragging  coefficient 

#define  kr  0.0 1 3  //rolling  coefficient 

#define  mev  1172  //mass  of  electric  vehicle 

#define  g  9.8  //acceleration  due  to  gravity 

#define  Af  2. 1 3  //Vehicle  frontal  Area 

#define  vhw  5.00  //Wind  Velocity 

#define  ad  1.29  //Air  Density 

#define  km  1 .0  //  rotational  inertia  coefficient 

//The  number  of  points  in  the  schedule.  This  is  identified  during 
//pre-processing,  after  reading  the  driving  schedule  file 
int  NumSch; 

float  MAX_RATIO,alpha; 

//A  structure  for  storing  the  time  instant  of  the  schedule  and  the  drive 
//power  required 
typedef  struct 


float  power; 
float  time; 


powerSch; 


101 


//Trapezoidal  integrator.  Takes  the  two  end  points  and  the  step 
//as  input 

inline  float  integr(float  x  1 , float  x2, float  delv) 

{ 

float  y; 

y=0.5  *(x  1  +x2)*del  v; 
retum(y); 

} 


//Class  drive  pre-declared 
class  drive; 

//Class  vehicle  declared  and  immediately  defined 
class  vehicle 
{ 


float  vev;//velocity  of  vehicle 
float  fev;//force  requirement  of  the  vehicle 
float  vevmax;//maximum  velocity  of  vehicle 
float  vevr;//rated  velocity  of  vehicle 
float  ta;//acceleration  time  to  rated  velocity 

public: 

void  inline  calc_force()//function  for  force  computation  given  velocity 

{ 

fev=mev*g*(kr+alpha)+0.5*ad*kd*AP(vev+vhw*0.4470)*(vev+vhw*0.4470); 


void  invalO 

{ 

cout«"Enter  maximum  velocity  of  car: 

cin  »vevmax; 

vevmax=vevmax*0.4470; 

cout«"Enter  rated  velocity  of  car: 

cin»vevr; 

vevr=vevr*0.4470; 

cout«"Enter  acceleration  time  of  car: 
cin»ta; 

} 

void  calc_rat(int  &k) 
k=vevmax/vevr; 

} 

//Outputs  the  road-load  characteristic 
void  calc_curve() 

{ 

FILE  *fp; 

fp=fopen("curve","w”); 
int  i=0; 

float  dk=vevmax/20; 
float  k=0; 

float  cruisePower=0; 
for  (i=0;i<=20;i++) 


vev=k; 

calc_force(); 

cruisePo  wer=fev*k/ 1 000; 

fprintf(fp,"%7.2f  %7.2f  %7.2f\n",(k/0.447),fev,cruisePower); 
k+=dk; 


fclose(fp); 

} 

//These  functions  do  computations  that  require  both  the  vehicle  as 
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//  well  as  the  drive  object 
friend  void  cruisLoss(vehicle, drive); 
friend  void  interface(vehicle, drive  &, float); 
friend  void  schedule(vehicle, drive  &); 


//Class  drive  defined 
class  drive 
{ 

float  vemr;//rated  speed  of  drive 

float  peinr;//rated  power  of  drive  based  on  acceleration  performance 
float  pemnnax;//drive  power  requirement  at  maximum  speed  of  vehicle 
powerSch  *ps;//The  power  requirement  of  drive  as  per  the  schedule 
public: 

void  displayO 

{ 

cout«"Power  Rating=  ’'«pemr«"  kW;  Rated  Speed=  ”«(vemr/0.4470)«"  mph"«endl; 
cout«endl; 

} 

void  showScheduleO 

{ 

int  i=0; 

FILE  *fp; 

fp=fopen("out","w"); 
while  (i<NumSch-l) 

{ 

fprintf(fp,"%7.2f  %7.2f\n",ps[i].time,ps[i].power); 
i++; 

} 

fclose(fp); 

} 

friend  void  cruisLoss(vehicle, drive); 
friend  void  interface(vehicle,  drive  &, float); 
friend  void  schedule( vehicle, drive  &); 


//Function  calculates  the  drive  losses  at  points  between  0-max  speed 
//  and  stores  in  a  file  cruisLoss 
void  cruisLoss(vehicle  c, drive  m) 

{ 

float  k=0; 

float  dk=c.vevmax/200; 
float  p; 

FILE  *fp; 

float  K1,K2E,K2H; 
fp=fopen("cruisLoss " , "  w" ) ; 
while  (k<=c.vevmax) 

{ 

//Remove  next  two  comments  when  you  want  at  constant  torque 

c.vev=k;/*c.vev=c.vevmax;*/ 

c,calc_force(); 

/*c.vev=k;*/ 
p=c.fev*c.vev/l 000; 

if  (c.vev>0)  Kl=(p*m.vemr)/(c.vev*m.pemr); 

K2E=(p*c.vev)/(m.pemr*m.vemr); 

If  ((c.vev>0)&&(p>0)&&(c.vev<=m.vemr))  K2H  =  ((exp(0.8*log(p/m.pemr))) 

*(exp(0.2*(c.vev/m.vemr)))); 

if  ((c,vev>0)&&(p>0)&&(c.vev>m.vemr))  K2H  =  ((exp(0.8*log(p/m.pemr)))  * 

(exp(0.2*(c.vev/m.vemr)))  *  (exp(0.5*log(m.vemr/c.vev)))); 
fprintf(fp,”%8.5f  %8.5f  %8.5f  %8.5f  %8.5f\n",p,(c.vev/0.4470),Kl,K2E,K2H); 
k+=dk; 

} 
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fclose(fp); 


} 

//Function  to  calculate  drive  peak-power  for  a  given  vehicle  performance 
void  interface(vehicle  c,drive  &m,  float  n) 

{ 

float  pmr,dpnir; 
float  vmr,dk; 
float  vmrmax; 
float  acc_time,ta; 
int  k,l,flag; 
float  v,vmin,xl,x2; 
ta=c.ta; 

c.vev=c.vevmax; 

vmnnax=c.  vevmax ; 

c.calc_force(); 

pmr=:c.fev*vmrmax; 

m.pemrmax^pmr; 

vmin=c.  vevmax/n ; 

dpmr=pmr/100; 

flag=0; 

while  (!flag) 

{ 

vmr=vmin; 

dk=c.vevr/1000; 

k=c.vevr/dk; 

acc_time=0.0; 

v=:0.0; 

c.vev=0; 

c.calc_force(); 

x2=vmr/(pmr-c.fev*vmr); 

for  (l=l;l<=k;l++) 


v=v+dk; 
if  (v<=vnir) 

{ 

c.vev=v; 

c,calc_force(); 

xl=x2; 

x2=vmr/(pmr-c.fev*vmr); 
acc„time=acc_time+integr(x  1  ,x2,dk); 


else 

{ 

c.vev=v; 

c.calc_force(); 

xl=x2; 

x2=v/(pmr-c.fev*v); 
acc_time=acc_time+integr(x  1  ,x2,dk); 

} 

} 

acc_time=acc_time*mev*km; 
if  ((acc_time<=ta)&&(acc_time>=:0)) 

{ 

flag=l; 

m.peinn=pmr/ 1 000; 
m.vemr=vmr; 

} 

pmr=pnir+dpmr; 


} 
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//Function  to  calculate  the  constants  Kl,K2e  and  K2h 
void  schedule(vehicle  c, drive  &m) 


char  drivFileName[100]; 

cout«"Please  enter  file-name  for  driving  schedule: "«endl; 
scanf("%s",drivFileName); 

FILE  *fp; 
float  veloc[ 10000]; 
float  time  [10000]; 
float  driveFMax; 
float  fChk; 
char  ch; 

powerSch  ps[  10000]; 
float  FI, F2; 
float  FEV=0; 
float  K1,K2E,K2H; 
int  i=0; 

NumSch=0; 

fp=fopen(drivFileName,"r”); 
while  ((ch=getc(fp))!=EOF) 

{ 

if  (ch=='\n’) 


NumSch++; 

} 


fclose(fp); 

//ps=(powerSch  *)malloc(NumSch*sizeof(powerSch)); 
//veloc=(float  *)malloc(NumSch*sizeof(float)); 
//time=(float  *)malloc(NumSch*sizeof(float)); 
fp=fopen(drivFileName,"r"); 
i=0; 

while  (i<NumSch) 

{ 

fscanf(fp, "  %f  %f\n "  ,&time  [i]  ,& veloc  [i]) ; 
i++; 

} 

for  (i=0;i<NumSch;i++) 

{ 

veloc[i]=veloc[i]  *0.447 ; 


fclose(fp); 

K1=0; 


K2E=0; 

K2H=0; 


for  (i=0; 


i<(NumSch-l);i++) 


Fl=mev*(veloc[i+l]-veloc[i])/(time[i+l]-time[i]); 

c.vev={veloc[i]+veloc[i+l])/2; 

c.calc_force(); 

FEV=c.fev; 

F2=FEV; 

driveFMax=m.pemr*  1 000/m.  vemr; 
if  (c.vev>m.vemr)  driveFMax=m.pemr*1000/c.vev; 

//The  braking  is  always  by  mechanical  means  and  motor  produces 
//no  power  during  this. 

//The  above  is  over-ridden.  Braking  is  always  electrical  now. 

ps  [i]  .po  wer=(F  1  +F2)*c.vev/ 1 000; 

ps[i].time=(time[i]+time[i+l])/2; 

fChk=:Fl+F2; 

if  (fChk<0)  fChk=-fChk; 

if  (fChk>driveFMax) 
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printf("At  Time  %7.2f',ps[i].time); 
printf("  drive  capability  violated\n"); 


} 

if  (veloc[i+l]<veloc[i]) 

{ 

if  (ps[i].power<0)  ps[i].power=-ps[i]. power; 

} 

if  (c.vev>0)  Kl=Kl+(ps[i].power/c.vev); 

K2E=K2E+ps[i].power*c.vev; 

if  ((c.vev>0)&&(ps[i].power>0)&&(c.vev<=m.vemr))  K2H  =  K2H  + 

((exp(0.8*log(ps[i].power/m.pemr)))  *  (exp(0.2*(c.vev/m.vemr)))); 
if  ((c.vev>0)&&(ps[i].power>0)&&(c.vev>m.vemr))  K2H  =  K2H  + 

((exp(0.8*log(ps[i].power/m.pemr)))  *  (exp(0.2*(c.vev/m.vemr)))  *  (exp(0.5*log(m.vemr/c,vev)))); 

} 

m.ps=ps; 

K 1  =K  1  *m.  veinr/(m,pemr*(NumSch- 1 )); 

K2E=K2E/(m.  vemr*m,pemr*(NumSch’  1 )); 

K2H=K2H/(NumSch-l); 

printf("%7.5f  %7.5f  %7.5fm",Kl,K2E,K2H); 

} 


int  main() 

{ 

vehicle  car; 
drive  srm; 
int  k; 
float  i; 

cout«"Enter  the  maximum  ratio:"; 
cin»MAX_RATIO; 

//MAX^RATIO  Max.  Ratio  Maximum/Rated  Speed  upto  which  to  compute 

i=MAX_RATIO; 

cout«" Enter  the  gradient:"; 

cin»alpha; 

//alpha  =  Gradient 

car.invalQ; 

interface(car,srm,i); 

cout  «"Maximum  Ratio  Maximum/Rated  Speed  of  Drive="«i«endl; 

srm.displayO; 

cruisLoss(car,srm); 

schedule(car,srm); 

car.calc_curve(); 

srm.  sho  wScheduleO ; 

} 


B.2  MATLAB™  Code  for  Sizing  Other  Motors 

A.  PMSM 


%Code  for  sizing  the  permanent  magnet  synchronous  machine 
clear  all 
clear, clc 


xdr= . 5 ; 
xqr= . 66 ; 
vr=35 ; 
vmax=7  0  ; 
E=0.8; 
pr=l; 


%the  d-axis  reactance  at  the  rated  drive  speed 
%the  q-axis  reactance  at  the  rated  drive  speed 
%the  electric  vehicle  rated  speed  is  14  mph 
%the  electric  vehicle  maximum  speed  is  70  mph 
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prated=49;%The  peak  power  rating  of  the  machine  is  known  to  be  28  kW  at  ratio 
of  5  and  37  kW  at  ratio  2from  optim.c 
f lag=0 ; 

[ idr , iqr ] =computen ( E , 1 , xdr , xqr , pr ) ; 

ir=sqrt{iqr^2+idr'^2)  ; 

kfeer= ( (E+xdr*idr) /E) ^2 ; 

kfehr=:{  (E+xdr*idr)  /E)^1.6; 

kcu=0 ; 

kfee=0; 

kfeh=0; 

C=0 . 4; 

f=0.3333; 

%%%%%%%For  driving  schedules%%%%%%%%%%%%%%%%%%%%%%%%%%%UNCOMMENT  DESIRED 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%SECTION 

fid=fopen( 'pshwy.ouf , 'r' ) ; 

[data, count] =fscanf {fid, ' %f ' ) ; 
fclose(fid) ; 
div= count/3 ; 

%%%%%%%For  cruising  at  maximum  velocity%%%%%%%%%%%%%%%%UNCOMMENT  DESIRED 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%SECTION 

%count=3 ; 

%data(l)=:l; 

%data (2 ) =20 ; 

%data(3)=55*0.4470; 

%count=3 ; 

%div=count/3 ; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
for  i=l : (count/3 ) 

k=data( (i-1) *3+3 ) /O . 4470/vr; 
if  k<=l 
if  k~=0 

[id, iq] =computen (E*k, k, k*xdr , k*xqr , (data ( (i-1) *3+2 ) /prated) ) ; 
i=sqrt  (id""2+iq'^2)  /ir; 

end 

if  k==0 

div=div-l; 

end 

else 

% [id, iq] =computen (E*k, 1+ ( (k-1) /4) , k*xdr , k*xqr , (data ( ( i-1) *3+2 ) /prated) ) ; 
[id, iq] =computen (E*k, 1 , k*xdr , k*xqr , (data ( ( i-1) *3+2 ) /prated) ) ; 
i=sqrt (id^2+iq^2 ) /ir; 

end 

if  k-'=0 

kcu=kcu+i'^2  ; 

end 

if  k-'=0 

kfee=kfee+  (  (  (  (E+xdr*id)  /E)  "2)*k'^2)/kfeer; 
kfeh=kfeh+(  (  (  (E+xdr*id)  /E)  "'1 . 6)  *k) /kfehr ; 

end 

end 

kcu=kcu/div 
kf ee=kf ee/div 
kfeh=kfeh/div 

w=kcu* (1-c) +kfee*c*f+kfeh*c* ( 1-f ) 
s=((w-c)/(l-c))'^.5 


function  [id, iq] =computen(Eb, V,Xd,Xq, P) 
id=0  ; 

delid=0.01; 

iq=P/ {Eb+ (Xd-Xq) *id) ; 
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f =iq'"2+  (  (Xd/Xq)  "2  )  *  (  ( id+  (Eb/Xd)  )  ^2 )  -  (V/Xq)  "2  ; 
while  abs(f)>0.01 
id=:id-delid; 
iq=P/ (Eb+ (Xd-Xq) *id) ; 

f=iq-2+{  (Xd/Xq)-2)*  (  (id+ (Eb/Xd)  )  "2 )  -  (V/Xq)  ^^2  ; 

end 


B.  IM 

%Code  for  sizing  the  induction  machine 
clear  all 
clear, clc 

vr=35;  %the  electric  vehicle  rated  speed  is  14  mph 
vmax=70;  %the  electric  vehicle  maximum  speed  is  70  mph 
pr=l; 

prated=49 ;%The  peak  power  rating  of  the  machine  is  known  to  be  28  kW  at  ratio 

%%%%%%%%%%%of  5  and  37  kW  at  a  ratio  of  2 from  optim.c 

kcu=0 ; 

kf ee=0 ; 

kfeh=0; 

c=0.4; 

f=0.3333; 

%%%%%%For  driving  schedules%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
fid=fopen( 'psftp.out ’ , 'r ' ) ; 

[data, count] =fscanf (fid, ' %f ' ) ; 
f close (fid) ; 
div=count/3 ; 

%%%%%%%For  cruising  at  maximum  velocity%%%%%%%%%%%%%%%% 

%count=3 ; 

%data(l)=:l; 

%data(2)=19; 

%data(3)=70*0.4470; 

%count=3; 

%div=count/3 ; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
for  i=l : (count/3 ) 

k=data( (i-l)*3+3)/0.4470/vr; 
p=data ( (i-1) *3+2) /prated; 
if  k==0 

div=div-l ; 

end 

if  k<=:l 
if  k-=0 

kcu=kcu+ (p/k) ^2 ; 
end 
else 

kcu=kcu+p'^2  ; 

end 

if  k<=l 
if  k-'=0 

kf  ee=kf  ee+k''2  ; 
kfeh=kfeh+k; 

end 

else 

kf ee=kfee+l ; 

kfeh=kfeh+( ( (1/k) . 6) *k)  ; 

end 

end 

kcu=kcu/div 
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kfee=kfee/div 

kfeh=kfeh/div 

w=kcu* (1-c) +kfee*c*f+kfeh*c* (1-f ) 
s=({w-c)/{l-c))^.5 


C  DC  Motor 

%Code  for  sizing  the  dc  machine 
clear  all 
clear , clc 

vr=35;  %the  electric  vehicle  rated  speed  is  14  mph 
vmax=70;  %the  electric  vehicle  maximum  speed  is  70  mph 
pr=l; 

prated=49;%The  peak  power  rating  of  the  machine  is  known  to  be  28  kw  form 

optim.c  at  ratio  of  5 

kcu=0; 

%%%%%%%For  driving  schedules%%%%%%%%%%%%%%%%%%%%%%%%%% 
fid=fopen ( 'ps ftp. out ' , 'r'); 

[data, count] =fscanf ( fid, ' %f ' ) ; 

f close (fid) ; 

div=count/3; 

%%%%%%%For  cruising  at  maximum  velocity%%%%%%%%%%%%%%%% 

%count=3 ; 

%data(l)=l; 

%data (2 ) =20 ; 

%data(3)=55*0.4470; 

%count=3 ; 

%div=count/3 ; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
for  i=l : (count/3 ) 

k=data{ ( i-1 ) *3+3 ) /O . 4470/vr ; 
p=data{ (i-l) *3+2) /prated; 
if  k==0 

div=div-l; 

end 

if  k<=l 
if  k-=0 

kcu=kcu+ (p/k)  ^"2  ; 
end 
else 

kcu=kcu+p^2 ; 

end 

end 

kcu=kcu/div 
w=kcu 
s=kcu'" .  5 
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Appendix  C 
ANSYS  CODE 


C.l.  6/4  SRM 

'This  file  calculates  an  electromagnetic  solution  of  the  SRM 

! Written  by:  Shyam  S.R. 

!on:  11/8/99 


! variables  of  the  motor  geometry: 


sor=0 . 2321 

smr=0 . 1775 

sir=0.1510 

ror=0 . 1490 

rmr=0 . 1043 

sr=0.0468 

spw=30 

rpw=32 

theta=-30 

nsp=6 

nrp=nsp-2 

cww=360/nsp-spw 

i agw= 3 6 0 / nrp-rpw 

mdepth=0 .2981 

current=258 

turns =16 


!  stator  outside  radius 
!  stator  middle  radius 
!  stator  inside  radius 
! rotor  outside  radius 
!  rotor  middle  radius 
!  shaft  radius 
! stator  pole  width 
•rotor  pole  width 

[position  of  the  rotor  with  respect  to  the  stator  (0  -  -45) 

[number  of  stator  poles 

[number  of  rotor  poles 

[copper  winding  width 

[ inner  airgap  width 

[depth  of  the  motor 

[current  applied  to  the  windings 

[niimber  of  windings 


[in  the  following,  l=phase  a  2=phases  a  and  b  3=phases  a  and  c 

phase=l  [which  phases  are  excited 

[calculate  the  angle  from  the  origin  to  the  top  of  the 
[small  angle  on  the  pole  of  the  rotor 


[numerator  =  num 

num=tan(0*3 . 14159265/180) * (ror-rmr) +ror*sin (rpw/2 *3 . 14159265/180) 
phi=rpw/2 [asin (num/rmr) *180/3.14159265 

[enter  the  preprocessor  of  the  program 

/prep7 

[define  the  element  types  to  be  used 

et, l,planel3 , 0 

[define  the  units  (metric) 


emunit,mks 

[associate  material  properties  (permeabilities)  with  element  types 

mp,murx,l,l  [air 

mp,murx,2,l  [steel 

mp , murx ,3,1  [ copper 

[activate  a  data  table  to  analyze  the  rotor  and  stator 
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Imaterial  properties  and  assign  those  properties  a  material  number 
TB,BH,4, ,33 

TBPT,defi,771030.21,3 .1 
TBPT, def i , 691452 .74,3.0 
TBPT,defi, 611875.74,2.9 
TBPT, def i , 532297 .80,2.8 
TBPT,def 1,452720.33, 2. 7 
TBPT, def i, 373142 .86,2.6 
TBPT, def 1,293565.39, 2. 5 
TBPT,defi, 213987 .91,2 .4 
TBPT, def i, 134410 .44,2 .3 
TBPT, defi, 77534. 21, 2. 2 
TBPT, def 1,43639.5,2.1 
TBPT,defi, 310*79.6, 2 
TBPT, defi, 190*79.6, 1.9 
TBPT, def 1,100*79.6,1.7875 
TBPT,defi, 80*79.6, 1.747917 
TBPT,defi, 56*79.6, 1.7 
TBPT, def 1,40*79.6, 1.658333 
TBPT, defi, 20*79 .6,1.6041667 
TBPT, def 1,10*79.6, 1.552083 
TBPT, def 1,8*79.6,1. 526042 
TBPT , defi ,6.4*79.6,1.5 
TBPT, defi, 5*79.6, 1.46875 
TBPT, def 1,4*79.6,1.4479167 
TBPT, def 1,3*79.6,1.395833 
TBPT, def 1,2 .1*79.6, 1.3 
TBPT, defi, 1.5*79. 6, 1.1 
TBPT, def 1 , 1 .1*79 .6,0.9 
TBPT, def 1 ,0.85*79.6,0.7 
TBPT, defi, 0.68*79.6, 0.5 
TBPT, defi, 0.6*79.6, .3854167 
TBPT , defi ,0.54*79.6,0.3 
TBPT, defi, 0.4*79.6, .1666667 
TBPT, def 1,0.33*79. 6, 0.1 
!TBPLOT,BH, 4 

•define  the  resistivity  of  the  copper  windings 
res=l . 512 

! associate  the  resistivity  with  the  material  number  of  copper 
mp, rsvz, 3 , res 

! activate  cylindrical  cooordinate  system 
csys , 1 

! define  nodes  used  in  flux  calculations 

n, 1 , sor , -spw/2 
n, 2 , sor, 0 
n, 3 , sor, spw/2 

n,  4,  smr,  -spw/2 
n , 5 , smr , 0 
n,  6,  smr,  spw/2 

n, 7 , sor, 180-spw/2 

n, 8, sor, 180 

n, 9, sor, -180+spw/2 
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n, 10 , smr , 180“Spw/2 

n, 11, smr, 180 

n,  12 ,  siur ,  -180+spw/2 

n, 13 , ror, rpw/2 
n, 14, ror, 0 
n, 15, ror, “rpw/2 

n, 16 , ror, 180“rpw/2 

n, 17 , ror, 180 

n, 18, ror, -180+rpw/2 

•define  the  keypoints  of  the  SRM 

•keypoints  of  the  stator  poles  and  copper  windings 

k, , sir, spw/2 

k,  ,  smr ,  spw/2 

k, , sir, spw/2+cww/2 

k, , smr, spw/2+cww/2 

k, , sir , spw/2+cww 

k, , smr, spw/2+cww 

ksel , s , , , 1 , 6 

kgen,nsp, all, , , , 360/nsp 

•keypoints  of  the  rotor  and  shaft 

k,  ,  ntir ,  theta-phi 
k, , ror, theta-rpw/2 
k, , ror, theta+rpw/2 

k,  , rmr , theta+phi 
ksel,s, , ,nsp*6+l,nsp*6+4 
kgen, nrp, all , , , , 360/nrp 

Hines  for  the  rotor  poles 

csys , 0 

l ,  6*nsp+l , 6*nsp+2 
1, 6*nsp+3 , 6*nsp+4 
1, 6*nsp+5 , 6*nsp+6 
1 , 6*nsp+7 , 6*nsp+8 
1 , 6*nsp+9 , 6*nsp+10 
1, 6*nsp+ll, 6*nsp+12 
1, 6*nsp+13 , 6*nsp+14 
1 , 6*nsp+15 , 6*nsp+16 
csys , 1 

! select  all  keypoints  of  the  SRM 

ksel , all 
kplot,all 

[generate  areas  of  the  SRM 

[areas  of  the  stator  ring 

pcirc, smr, sor , 0 , spw 

agen, 2*nsp, all, , , ,360/ (2*nsp) 

[areas  of  the  airgap  ring 

pcirc, ror , sir, 45+theta, 135+theta 
agen , nrp , 2  *nsp+l ,,,,360 /nrp 
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‘areas  of  the  rotor  ring 

pcirc, sr , rmr , theta, 360/nrp+theta 
agen,nrp, 2*nsp+nrp+l, , , ,360/nrp 

!area  of  the  shaft 

pcirc, , sr, theta, 360+theta 

! areas  of  the  stator  poles 

a, 1, 2 , nsp*6 ,nsp*6-l 

agen, nsp, 2*nsp+2*nrp+2 , , , , 360 /nsp 

! areas  of  the  copper  windings 


a, 1,2, 4, 3 
a, 3 , 4, 6, 5 

agen , nsp , 3  *nsp+2  *nrp+2 , 3  *nsp+2  *nrp+3 ,,,36  0 /nsp 

! areas  of  the  inner  airgap  area 

a, nsp* 6+3 , nsp* 6+4 , nsp* 6+5 , nsp* 6+6 
agen, nrp, 5*nsp+2*nrp+2 , , , , 360/nrp 

! areas  of  the  rotor  poles 

a, nsp* 6+1 , nsp* 6+2 , nsp* 6+3 , nsp* 6+4 
agen , nrp , 5  *nsp+3  *nrp+2 ,,,,360 /nrp 

! merge  coincident  keypoints,  lines,  etc, 

ksel, all 
nximmrg ,  kp 

!glue  all  correct  areas 

asel , s, loc, X, smr , sor 
aglue, all 

asel, s, loc,x, sir,  smr 
aglue , all 

asel , s , loc , X , rmr ,  smr 
Isel, s, loc,x,  smr 
asll ,u 

Isel, s, loc, X, rmr 
asll, u 
aglue, all 

asel , s, loc, X, rmr , ror 
Isel, s, loc,x, sir 
asll , u 
aglue, all 

asel , s , loc , x, 0 , rmr 
aglue, all 

! group  commom  areas  of  the  SRM  together  under  a 
! component  name. 

! stator  ring 

Isel, s, loc,x, sor 
asll 

cm,  stator, area 
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aatt, 4, ,  1 


{windings 
!phase=a  (1) 

* i f , phase , eq , 1 , then 

{positive  copper  windings  of  interest 

asel, s, loc,x, sir, smr 
asel , r , loc,y, spw/2+cww/4 
asel , a, loc ,y, 180“Spw/2-cww/4 
asel,r, loc,x, sir, smr 
cm,pwind,  area 
e  -t,3,,l 

{negative  copper  windings  of  interest 

asel, s, loc,x, sir, smr 
asel, r, loc,y, -spw/2-cww/4 
asel , a, loc,y, -180+spw/2+cww/4 
asel , r , loc, x, sir, smr 
cm,nwind, area 
aatt, 3 , , 1 

{non  interest  copper  windings 

asel, s, loc,x, sir, smr 

asel, r, loc, y, -180+spw/2+ . 25*cww 

*do,i,l,nsp 

n=-180+spw/2+.25*cww+ (360/nsp) *i 
asel, a, loc,y,n,n 
*  enddo 

asel, a, loc,y, -180+spw/2+ . 75*cww 
*do, i, 1, nsp 

n=-180+spw/2+. 75*cww+ (360/nsp) *i 
asel , a, loc,y, n, n 
* enddo 

asel,r, loc,x, sir, smr 
cmsel , u,pwind 
cmsel , u, nwind 
cm,  owind,  area 
aatt, 3 , , 1 

{phase=ab  (2) 

*elseif, phase, eq, 2 

{positive  copper  windings  of  interest 

asel, s, loc,x, sir, smr 
asel, r, loc,y, spw/2+cww/4 
asel, a, loc,y, 180“Spw/2“Cww/4 
asel, a, loc,y, 1 . 5*spw+l . 25*cww 
asel , a, loc, y, ~180+spw/2+ . 75*cww 
asel, r, loc, x, sir , smr 
cm,pwind,area 
aatt, 3 , , 1 

{negative  copper  windings  of  interest 
asel, s, loc, X, sir, smr 
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asel, r , loc,y, -spw/2~cww/4 

asel, a, loc,y, -180+spw/2+cww/4 

asel, a, loc,y, spw/2+ . 75*cww 

asel , a, loc,y, -180+1 . 5*spw+l . 25*cww 

asel , r, loc, x, sir, smr 

cm,  nwind,  area 

aatt, 3 , ,  1 

!non  interest  copper  windings 

asel, s, loc,x, sir, smr 

asel , r , loc,y, -180+spw/2+ . 25*cww 

*do, i , l,nsp 

n=-180+spw/2+.25*cww+ (360/nsp) *i 

asel, a, loc,y,n,n 

*enddo 

asel , a, loc,y, -180+spw/2+ . 75*cww 
*do, i, l,nsp 

n=-180+spw/2+ . 75*cww+ (360/nsp) *i 

asel , a, loc, y, n, n 

*enddo 

asel, r, loc,x, sir, smr 
cmsel  ,u,pwind 
cmsel,u, nwind 
cm, owind, area 
aatt, 3 , , 1 

!phase=ac  (3) 

*elseif , phase, eq, 3 

! positive  copper  windings  of  interest 

asel, s, loc,x, sir, smr 
asel,r, loc,y, spw/2+cww/4 
asel, a, loc ,y, 180-spw/2-cww/4 
asel , a, loc,y, 180-spw/2“ . 75*cww 
asel , a, loc,y, -1 . 5*spw-l . 25*cww 
asel , r , loc, x, sir,  smr 
cm,pwind,area 
aatt, 3 , , 1 

! negative  copper  windings  of  interest 

asel , s , loc, X, sir, smr 
asel , r , loc,y, -spw/2-cww/4 
asel, a, loc,y, -180+spw/2+cww/4 
asel, a, loc,y, 180-1 . 5*spw-l . 25*cww 
asel, a, loc,y, - . 5*spw- .75*cww 
asel, r, loc,x, sir, smr 
cm,  nwind,  area 
aatt, 3 , , 1 

!non  interest  copper  windings 

asel, s, loc,x, sir, smr 

asel, r, loc,y, -180+spw/2+ . 25*cww 

*do, i, l,nsp 

n=-180+spw/2+.25*cww+ (360/nsp) *i 

asel , a, loc, y, n, n 

*enddo 

asel, a, loc,y, -180+spw/2+ . 75*cww 
*do, i, l,nsp 


n=-180+spw/2+ . 75*cww+ (360/nsp) *i 

asel, a, loC/y,n,n 

*enddo 

asel,r, loc,x, sir, smr 
cmsel,u,pwind 
cmsel,u,nwind 
cm, owind, area 
aatt, 3 , , 1 

*endif 

! calculate  the  area  of  the  positive  windings  and 
•assign  the  value  to  the  variable  poswind 

cmsel, s,pwind 
asum 

*get, poswind, area, ,area 


•calculate  the  area  of  the  negative  windings  and 
! assign  the  value  to  the  variable  negwind 

cmsel , s , nwind 
asum 

*get, negwind, area, ,area 

! stator  poles: 

asel,s, loc,x,sir,  smr 
cmsel , u,pwind 
cmsel,u, nwind 
cmsel , u, owind 
cmsel , u , s  tator 
cm, spole, area 
aatt, 4, , 1 

! shaft 

asel , s , loc,  X,  0 
cm, shaft, area 
aatt , 2 , , 1 

lairgap  ring 

Isel , s , loc, X, sir 
asll 

cmsel , u, spole 
cmsel,u, nwind 
cmsel, u,pwind 
cmsel,u, owind 
cm,  airgap, area 
aatt, 1, , 1 

! rotor  ring 
Isel, s , loc, X, sr 
asll 

cmsel, u,  shaft 
cm,  rotor , area 
aatt, 4, , 1 

•all  inner  airgap  areas 
asel , s , loc , x, rmr , ror 
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asel, r , loc,y, ”180+rpw/2+iagw/2+theta 
*do, i, l,nrp 

n=-180+rpw/2+iagw/2+ (360/nrp) *i+theta 

asel,a, loc,y,n,n 

*enddo 

asel ,  r,  loc,X/  ror 
cmsel , u, airgap 
citi;  inag,  area 
aatt , 1, ,1 

! rotor  poles 
Isel,  s,  loc,x,nnr 
asll 

cmsel, u, rotor 
cmsel, u, inag 
cm, rpole, area 
aatt, 4, , 1 

Isize  lines  for  meshing 

cmsel, s, airgap 
Isla 

Isel , r , loc,x, sir 
lesize, all, , , 160 
Isla 

lsel,r, loc,x,ror 
lesize, all, , , 160 

cmsel, s, stator 
Isla 

lsel,r, loc,x, sor 
lesize, all , , , 8 
Isla 

Isel , r , loc, X, smr 
lesize, all, , , 8 
Isla 

lesize, all , , , 3 , . 8 

cmsel, s, spole 
Isla 

Isel , r , loc, X, sir 
lesize, all , , , 18 
Isla 

Isel , r , loc,x, smr 
lesize, all, , , 18 

cmsel , s , owind 
cmsel,a,nwind 
cmsel, a,pwind 
Isla 

Isel, r , loc, X, sir 
lesize, all, , , 9 
Isla 

Isel , r , loc, X, smr 
lesize, all , , , 9 

Isel, s, loc, y, -180+spw/2+cww 
*do, i, l,nsp 

n=-180+spw/2+cww+ (spw+cww) *i 

Isel , a, loc,y,n, n 

*enddo 

Isel, r , loc, X, sir, smr 
lesize, all , , , 16 , . 6 


Isel , s , loc ,y, “180+spw/2 
*do, i, l,nsp 

n=“-180+spw/2+ (spw+cww)  *i 
Isel, a, loc,y,n, n 
*enddo 

lsel,r,  loc,x,  sir, sitir 
lesize, all , , ,16,1.66667 

cmsel,  s,  owind 
cmsel , a, nwind 
cmsel , a,pwind 
Isla 

lesize, all, , , 16, . 6 

cmsel , s, rpole 
Isla 

Isel, r, loc,x, ror 
lesize, all, , ,20 
Isla 

Isel, r, loc,x, rmr 
lesize, all, , ,20 
Isla 

lesize, all , , , 12 

cmsel, s, inag 
Isla 

Isel, r, loc, X, ror 
lesize, all, , ,20 
Isla 

Isel, r, loc, X, rmr 
lesize, all, , ,20 

cmsel , s , rotor 
Isla 

Isel, r, loc,x, rmr 
lesize, all, , , 12 
Isla 

Isel , r, loc,  X,  sr 
lesize, all, , , 12 
Isla 

lesize, all , , , 5 , . 8 

allsel 
amesh, all 

! constraint  equations 

! stator  ring  to  poles  and  windings 

allsel 
esel,none 
cmsel , s , stator 
Isla 

Isel, r, loc,x, smr 

nsll , s 

esln 

cmsel, s, spole 
cmsel , a, owind 
cmsel, a, nwind 
cmsel, a,pwind 
Isla 

Isel, r, loc,x, smr 


118 


nsll ,  a 
ceintf 

lairgap  to  stator  poles  and  windings 

allsel 
esel,none 
cmsel, s, spole 
cmsel , a, owind 
cmsel,  a,nwind 
cmsel , a,pwind 
Isla 

Isel , r, loc,X/ sir 

nsll, s 

esln 

cmsel, s, airgap 
Isla 

Isel, r, loc,x, sir 

nsll , a 

ceintf 

lairgap  to  rotor  poles  and  inner  airgap  region 

allsel 

esel,none 

cmsel, s, rpole 

cmsel , a, inag 

Isla 

Isel , r, loc, X, ror 

nsll, s 

esln 

cmsel , s , airgap 
Isla 

lsel,r, loc, X, ror 

nsll, a 

ceintf 

! rotor  poles  and  inner  airgap  region  to  rotor  yoke 

allsel 

esel , none 

cmsel, s, rotor 

Isla 

Isel, r, loc,x,rmr 

nsll, s 

esln 

cmsel , s , rpole 
cmsel, a, inag 
Isla 

Isel , r, loc,x,  rmr 

nsll , a 

ceintf 

finish  !the  preprocessor  is  finished 

/solution 

! the  solution  sequence  is  to  solve  the  problem  over  two  load  steps 
! the  first  load  step  obtains  an  approximate  solution 
! the  second  load  step  obtains  a  converged  solution 
! the  analysis  is  static  by  default 
!the  Newton-Raphson  option  is  by  default 

! unselect  all  components 
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cmsel , none 


! select  the  components  of  the  SRM  that  a  current  density  will  be  applied 
! to  and  then  select  all  elements  of  those  components.  After  this,  the 
!bfe  command  applies  the  current  density  to  the  selected  elements. 

!phase=a  (1) 

* i f , phase , eq , 1 , then 

cmsel , s , nwind 
esla, s 

bfe, all, js, 3 , “Current* turns*2/negwind  !“Z  direction 

cmsel , none 
cmsel , s , pwind 
esla, s 

bfe, all, js, 3 , current*turns*2/poswind  !+Z  direction 

!phase=ab  (2) 

*elseif, phase, eq, 2 

cmsel, s, nwind 
esla, s 

bfe, all, js, 3 , -current* turns *4 /negwind  !-Z  direction 

cmsel , none 
cmsel,  s, pwind 
esla, s 

bfe,all, js, 3, current* turns *4 /poswind  !+Z  direction 

!phase=ac  (3) 

*elseif, phase, eq, 3 

cmsel , s, nwind 
esla, s 

bfe, all, js, 3 , -current* turns *4 /negwind  ! -Z  direction 

cmsel , none 
cmsel , s, pwind 
esla, s 

bfe, all, js, 3 , current* turns *4 /poswind  !+Z  direction 

*endif 


!far-field  boundary  conditions 
allsel 

nsel , s , loc , x, sor 
d, all , az , 0 


! select  all  entities  of 
allsel 

! first  load  step 

neqit ,  1 
nsubst, 5 
outres , al 1 , none 
solve 

! second  load  step 


the  SRM  for  the  solution 


!  1  equilibrium  iteration 
!  5  substeps 

Ido  not  save  results  data 
!  solve  the  first  load  step 
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neqit, 20 
nsubst , 1 
cnvtol, csg, , le-7 
outres , all, last 
solve 


!up  to  20  equilibrium  iterations 
!  1  substep 

!  convergence  criteria 
Isave  converged  solution 
!  solve  the  second  load  step 


finish 


!  the  solution  is  finished 


! start  the  postprocessing 
/postl 

! return  results  in  cylindrical  coordinate  system 
rsys ,  1 

•calculate  the  stored  energy  of  the  srm 

!argl=material  #  of  first  material 
!arg2 -material  #  of  last  material 
!arg3= increment  of  the  material  # 

! units  are  per  length 

allsei 


/nopr 

parsav, all , sasimac 
esel, s,mat, ,1,4,1 
etable, sene, sene 
/out, scratch 
ssum 
/out 

*get, energy, ssum, , item, sene 
/gopr 

menergy= energy *mdepth 

esel, s, mat, ,1,4,1 
senergy, 1 , 0 
c  o  en= c_eng  *mdep  th 

esel, s, mat, ,1,4,1 
senergy, 0 , 0 
ener=s_eng*mdepth 

lnkage= (ener+coen) /current 

! calculate  the  airgap  torque 

allsei 

cmsel , s , airgap 
esla 

torqc2d, ror+ ( (sir-ror) /2) , 100 
agtorque=torque*mdepth 

parsav, scalar , parameters 

Iplots 

!flux  lines 

allsei 

!plf2d,54 
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!  azr 

cmsel , s , rpole 
cmsel , a, rotor 
cmsel/ a, shaft 
nsla 

nsel / a, loc, X, 0 , rmr 
plnsol / a, z 


C.2.  8/6  SRM 


The  code  for  constructing  the  geometry  and  executing  the  run  is  almost  the  same  as  for  the 
three-phase  6/4  SRM.  The  difference  is  in  the  initial  assignment  of  dimensions,  angle,  number  of 
stator  and  rotor  poles,  the  selection  of  winding  areas,  meshing  and  assignment  of  negative  or 
positive  current  sheets.  By  current  sheet,  we  mean  the  product  of  coil  turns  and  coil  current  that 
will  be  encountered  in  the  cross-section  of  the  coil  in  the  2D  FEA  model. 


A.  Assignment  of  Dimensions  and  Other  Values 


sor=0 , 1851 

smr=0 . 1524 

sir=0.1195 

ror=0.1183 

nnr=0 .0828 

sr  =0.0477 

spw=22.5 

rpw=24 . 5 

nsp=8 

nrp=nsp-2 

cww=360/nsp-spw 

i agw= 3 6 0 / nrp - r pw 

mdepth=0 .2366 

arg2=131.1517 

arg3=131.1517 

theta=“22 . 5 

turns=25 


! stator  outside  radius 
! stator  middle  radius 
! stator  inside  radius 
'.rotor  outside  radius 
! rotor  middle  radius 
'shaft  radius 
! stator  pole  width 
! rotor  pole  width 
{number  of  stator  poles 
Iniomber  of  rotor  poles 
{copper  winding  width 
!  inner  airgap  width 
{depth  of  the  motor 

{current  working— | current  in  two  adjacent  phases 
'current  leading- [under  short-flux-path  mode 
{angle  of  rotor  pole  wrt  stator  pole 
{number  of  windings 


B,  Selection  of  Winding  A  reas 


{windings 

{positive  copper  windings  of  interest 

asel, S/ loc,x, sir, smr 
asel/ r, 1oc,y, spw/2+cww/4 
asel/a, loC/Y/ 180-spw/2-cww/4 
asel, r, loc,x, sir, smr 
cm, pwind2 , area 
aatt, 3 , ,1 

asel , s , loc , x , s ir , smr 

asel, r, loc,y, spw/2+ . 75*cww 

asel, a, loc,y, -180+1 . 5*spw+l . 25*cww 

asel, r, loc,x, sir,  smr 

cm, pwind3 , area 
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aatt , 3 , , 1 


! negative  copper  windings  of  interest 

asel,  s,  loc,x,  sir,  snir 
asel,  r,  loc,y,  “Spw/2-cww/4 
asel , a, loc,y, -180+spw/2+cww/4 
asel, r, loc, x, sir, smr 
cm, nwind2 , area 
aatt, 3, ,1 

asel, s, loc,x, sir,  smr 

asel,r, loc,y, 1.5*spw+1.25*cww 

asel, a, loc,y, -180+spw/2+ , 75*cww 

asel , r , loc,x, sir, smr 

cm, nwind3 , area 

aatt, 3 , , 1 

!non  interest  copper  windings 

asel , s, loc, X, sir, smr 

asel, r, loc,y, “180+spw/2+ . 25*cww 

*do,i,l,nsp 

n=-’180+spw/2  + .  25*cww+  (360/nsp)  *i 

asel, a, loc,y,n,n 

*enddo 

asel, a, loc,y, -180+spw/2+ . 75*cww 
*do, i , l,nsp 

n=-180+spw/2+. 75*cww+ (360/nsp) *i 

asel , a, loc, y, n, n 

*enddo 

asel, r, loc,x, sir, smr 
cmsel , u,pwind2 
cmsel , u , pwind3 
cmsel , u , nwind2 
cmsel , u , nwind3 
cm, owind, area 
aatt, 3 , , 1 

! calculate  the  area  of  the  positive  windings  and 
! assign  the  value  to  the  variable  poswind 

cmsel , s,pwind2 
asum 

*get , poswind , area , , area 


! calculate  the  area  of  the  negative  windings  and 
! assign  the  value  to  the  variable  negwind 

cmsel , s , nwind2 
asum 

*get, negwind, area, ,area 


C  Meshing 

!size  lines  for  meshing 

cmsel , s , stator 
Isla 

lsel,r, loc,x, sor 
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lesize, all , , 10 
Isla 

Isel, r, loc,x, smr 
lesize^all, , 10 

cmsel , s , spole 
Isla 

Isel, r, loC;X, sir 
lesize, all , , 2 
Isla 

Isel,  r,  loc,x,  smr 
lesize, all, , 2 

cmsel,  s,  owind 
cmsel , a , nwind2 
cmsel , a , nwind3 
cmsel , a , pwind2 
cmsel , a , pwind3 
Isla 

lsel,r, loc,x, sir 
lesize, all, , 2 
Isla 

Isel , r , loc, X, smr 
lesize, all , , 2 

Isel, s, loc,y, “180+spw/2+cww 
*do, i, l,nsp 

n=’-180+spw/2+cww+ (spw+cww)  *i 

Isel, a, loc,y,n,n 

*enddo 

Isel , r, loc, X, sir, smr 
lesize, all, , , 12 ,  . 6 

Isel, s, loc,y, -180+spw/2 
*do, i , 1 , nsp 

n=-180+spw/2+ (spw+cww) *i 
Isel, a, loc, y, n, n 
*enddo 

Isel , r , loc, X, sir,  smr 
lesize, all, , ,12,1,6667 

Isel, s, loc,y, “180+spw/2+cww/2 
*do, i, l,nsp 

n=:-180+spw/2+cww/2+  (spw+cww)  *i 

Isel, a, loc,y,n,n 

*enddo 

Isel, r, loc,x, smr, sor 
lesize, all, , , 3 

Isel, s , loc,y, -180 
*do, i , 1, nsp 
n=-- 180+  (spw+cww)  *i 
Isel, a, loc,y,n,n 
*enddo 

Isel, r, loc,x, sir, smr 
lesize, all, , , 12 ,  .  6 

cmsel,  s, owind 
cmsel , a, nwind2 
cmsel, a,nwind3 
cmsel , a,pwind2 
cmsel , a,pwind3 
Isla 
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lesize, all ,  ,  ,\2,  .6 


cmsel , s , rpole 
Isla 

Isel, r, loc,x,ror 
lesize, all , , 3 
Isla 

lsel,r,  loc,x,nnr 
lesize, all, , 3 

cmsel, s, rpole 
Isla 

Isel , u, loc,x, ror 
lsel,u,  loc,x,nnr 
k=l 

*do, i, l,nrp 
Isel , u, line, k, k 
k=k+2  *nrp 
*enddo 

lesize, all , , ,12,1.25 

cmsel, s, rpole 
Isla 

lsel,u, loc,x,ror 
Isel , u, loc, X, rmr 
k=2 

*do, i, l,nrp 
lsel,u, line, k, k 
k=k+2*nrp 
*enddo 

lesize, all, , , 12 ,  .8 

cmsel, s , inag 
Isla 

Isel , r, loc, X, ror 
lesize, all, , 3 
Isla 

Isel , r, loc, X, rmr 
lesize, all , , 3 

cmsel , s , inag 
Isla 

Isel, u, loc, X, ror 
Isel, u, loc, X, rmr 
k=l 

*do, i, l,nrp 
lsel,u, line, k,k 
k=k+2*nrp 
*enddo 

lesize, all, ,,12,1.1 

cmsel, s, inag 
Isla 

Isel , u, loc,x, ror 
Isel , u, loc, X, rmr 
k=2 

*do, i, l,nrp 
Isel , u, line, k, k 
k=k+2*nrp 
*enddo 

lesize, all, , , 12 ,  .  9 
cmsel , s , rotor 
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Isla 

Isel ,r, loc, X,  rmr 
lesize, all, , 5 
Isla 

Isel , r , loc, X, sr 
lesi2e,all, , 5 

cms el, s, rotor 
Isla 

lsel,u, loc,x, rmr 
lsel,u, loc,x, sr 
lesize, all , , , 8 

cmsel, s,airgap 
Isla 

Isel , r , loc, X, sir 
lesize, all, , .25 
Isla 

Isel, r , loc , X, ror 
lesize, all , ,0.25 
Isla 

lsel,u, loc,x, sir 
Isel , u, loc, X, ror 
lesize, all , , , 2 

allsel 
amesh, all 

D.  Assignment  of  Current  Sheets 


! unselect  all  components 


cmsel , none 

! select  the  components  of  the  SRM  that  a  current  density  will  be  applied 
!tO  and  then  select  all  elements  of  those  components.  After  this,  the 
!bfe  command  applies  the  current  density  to  the  selected  elements. 


cmsel , s,nwind2 
esla, s 

bfe, all, js, 3 , -arg2*turns*2/negwind  !-Z  direction 
cmsel , none 
cmsel, s,pwind2 
esla, s 

bfe,all, js, 3 , arg2*turns*2/poswind  !+Z  direction 


cmsel , s , nwind3 
esla, s 

bfe, all, js, 3 , -arg3*turns*2/negwind  !-Z  direction 
cmsel , none 
cmsel , s,pwind3 
esla, s 

bfe, all, js, 3 , arg3*turns*2/poswind  !+Z  direction 


C.3.  10/8  SRM 

The  code  is  similar  to  that  for  the  four-phase  8/6  SRM.  The  main  difference  is  in  the 
assignment  of  values  and  meshing  section. 
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A.  Assignment  of  Dimensions  and  Other  Values 


sor=0 . 1856 

smr=0 . 1594 

sir=0.1195 

ror=0.1183 

rmr=0.0828 

sr  =0.054 

spw=18 

rpw=2  0 

nsp=10 

argl=-18 

arg2=131.153 

arg3=131.153 

nrp=nsp“2 

cww=360/nsp-spw 

iagw=3  60 /nrp-rpw 

mdepth=0 .2366 

turns=25 


!  stator  outside  radius 
!  stator  middle  radius 
!  stator  inside  radius 
!  rotor  outside  radius 
!  rotor  middle  radius 
!  shaft  radius 
!  stator  pole  width 
!  rotor  pole  width 
! number  of  stator  poles 
!  Theta 

!  Current  working 
[Current  leading 
[number  of  rotor  poles 
[copper  winding  width 
[  inner  airgap  width 
[depth  of  the  motor 
[number  of  windings 


B.  Meshing 

[size  lines  for  meshing 

cmsel , s , airgap 
Isla 

Isel, r, loc,x, sir 
lesize, all , ,0.25 
Isla 

Isel, r, loc,x, ror 
lesize, all, ,0.25 
Isla 

lesize, all, , , 2 

cmsel, s, stator 
Isla 

Isel , r , loc, X, sor 
lesize, all, , 10 
Isla 

Isel, r, loc,x, smr 
lesize, all , , 10 

cmsel , s , spole 
Isla 

Isel , r , loc, X, sir 
lesize, all, , 2 
Isla 

Isel, r, loc,x, smr 
lesize, all, , 2 

cmsel , s , owind 
cmsel , a, nwind2 
cmsel , a,nwind3 
cmsel , a , pwind2 
cmsel , a , pwind3 
Isla 

Isel, r, loc,x, sir 
lesize, all, , 2 
Isla 

Isel, r, loc,x, smr 
lesize, all, , 2 


127 


Isel ,  s,  loc,y,  “180+spw/2+cww 
*do,  i ,  1,  nsp 

n=-180+spw/2+cww+ (spw+cww)  *i 

Isel, a; loc,y,n,n 

*enddo 

Isel , r , loc, X, sir, smr 
lesize, all, , , 12 , , 6 

Isel, s, loc, y, -180+spw/2 
*do, i, l,nsp 

n=-180+spw/2+ (spw+cww) *i 
Isel , a, loc, y, n, n 
*enddo 

Isel , r, loc,x, sir, smr 
lesize, all, , ,12,1.6667 

Isel, s, loc,y, -180+spw/2+cww/2 
*do, i, l,nsp 

n=-180+spw/2+cww/2+ (spw+cww) *i 

Isel, a, loc,y,n,n 

*enddo 

Isel , r, loc, X, smr, sor 
lesize, all , , , 3 

Isel, s, loc, y, -180 
*do, i , 1 , nsp 
n=- 180+ (spw+cww) *i 
Isel, a, loc,y,n,n 
*enddo 

Isel, r, loc,x, sir, smr 
lesize, all, , , 12 , . 6 

cmsel, s, owind 
cmsel , a , nwind2 
cmsel, a,nwind3 
cmsel, a, pwind2 
cmsel , a , pwind3 
Isla 

lesize, all, , , 12 ,  .  6 

cmsel, s, rpole 
Isla 

Isel , r, loc,x, ror 
lesize, all , , 3 
Isla 

Isel, r, loc,x, rmr 
lesize, all, , 3 

cmsel , s, rpole 
Isla 

lsel,u, loc,x, ror 
lsel,u, loc, X,  rmr 
k=l 

*do, i, l,nrp 
lsel,u, line,k,k 
k=k+2*nrp 
*enddo 

lesize, all, , ,12,1.25 

cmsel, s, rpole 
Isla 

Isel , u, loc, X, ror 


128 


Isel, u, loc, X,  rmr 
k=2 

*do, i, l,nrp 
Isel, u, line, k, k 
k=k+2*nrp 
*enddo 

lesize,all, , , 12  ,  . 8 

cmsel , s , inag 
Isla 

Isel, r, loc,x, ror 
lesize, all , , 3 
Isla 

Isel, r, loc,x, rmr 
lesize, all, , 3 

cmsel , s, inag 
Isla 

Isel, u, loc, X, ror 
lsel,u, loc, X, rmr 
k=l 

*do, i, l,nrp 
lsel,u, line,k,k 
k=k+2  *nrp 
*enddo 

lesize, all, ,,12,1.1 

cmsel, s, inag 
Isla 

lsel,u, loc,x, ror 
Isel , u, loc, X, rmr 
k=2 

*do, i, l,nrp 
lsel,u, line,k,k 
k=k+2*nrp 
*enddo 

lesize, all, , , 12 , . 9 

cmsel , s , rotor 
Isla 

Isel, r, loc,x, rmr 
lesize, all , , 5 
Isla 

Isel , r, loc,x, sr 
lesize, all, , 5 

cmsel , s, rotor 
Isla 

lsel,u, loc,x, rmr 
Isel , u, loc, X, sr 
lesize, all, , , 8 

allsel 
amesh, all 


C.4.  12/8  SRM 

The  main  difference  in  the  code  for  the  single-phase-excited  three-phase  12/8  SRM  is  that 
quadrature  current  sheets  carry  same  direction  of  current  resulting  in  local  short  paths  for  the  flux. 
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The  assignment,  meshing,  winding  areas  selection  and  the  current  sheet  application  code  are 
listed  below. 


A.  Assignment  of  Dimensions  and  Other  Values 


!This  file  calculates  an  electromagnetic  solution  of  the  SRM 
'This  is  the  file  for  the  12/8  SRM  using  single-phase  excitation 
! Quadrature  current  sheets  carry  same  direction  current.  This 
'corresponds  to  single  phase  excitation  case  of  6/4  SRM. 

!The  self  and  mutual  inductances  in  these  quadrature  phases  together 
•contribute  to  the  net  torque. 

! Written  by:  Shyam  S.R. 

!on:  11/6/99 


! variables  of  the  motor  geometry: 


sor=0 . 1876 

smr=0 . 1657 

sir=:0.1195 

ror=0.1183 

rmr=0 .0828 

sr=0 . 0583 

spw=15 

rpw=17 

theta=-15 

22.5) 

nsp=12 

nrp=nsp-4 

cww= 3 6 0 / nsp - spw 

iagw= 3 6 0 / nr p - r pw 

mr  ipth^O  .2366 

t  ^rent=131 . 1527 

turns =2 5 


I  stator  outside  radius 
!  stator  middle  radius 
!  stator  inside  radius 
!  rotor  outside  radius 
!  rotor  middle  radius 
!  shaft  radius 
!  stator  pole  width 
! rotor  pole  width 

! position  of  the  rotor  with  respect  to  the  stator  (0  -  - 

Inumber  of  stator  poles 
Inumber  of  rotor  poles 
I  copper  winding  width 
I  inner  airgap  width 
! depth  of  the  motor 
I  current  applied  to  the  windings 
Inumber  of  windings 


B.  Selection  of  Winding  Areas 
I windings 
!phase=a  (1) 

* i f , phase , eq , 1 , then 

I positive  copper  windings  of  interest 

asel/ s, loc, X, sir, smr 
asel , r , loc, y, spw/2+cww/4 
asel , a, loc,y, 90-spw/2-cww/4 
asel , a, loc,y, 180+spw/2+cww/4 
asel, a, loc,y, -90-spw/2-cww/4 
asel, r , loc, x, sir, smr 
cm ,  pwind ,  area 
aatt , 3 , , 1 

I  negative  copper  windings  of  interest 

asel, s, loc, X, sir, smr 
asel , r , loc,y, “Spw/2-cww/4 
asel , a, loc,y, 90+spw/2+cww/4 
asel, a, loc,y, 180-spw/2-cww/4 
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asel, a, loc,y/ -90+spw/2+cww/4 
asel, r, loc,x, sir, smr 
cm, nwind, area 
aatt, 3 , , 1 

!non  interest  copper  windings 

asel, s, loc,x, sir, smr 

asel, r, loc,y, -180+spw/2+ . 25*cww 

*do, i , 1 , nsp 

n=-180+spw/2+.25*cww+ (3 60 /nsp) *i 

asel, a, loc,y,n,n 

*enddo 

asel, a, loc,y, -180+spw/2+ . 75*cww 
*do, i, l,nsp 

n=-180+spw/2+ . 75*cww+ (360/nsp) *i 

asel, a, loc,y,n,n 

*enddo 

asel, r, loc,x, sir, smr 
cmsel , u,pwind 
cmsel,u, nwind 
cm,  owind,  area 
aatt, 3 , , 1 

!phase=ab  (2) 

*elseif ,phase, eq, 2 

‘positive  copper  windings  of  interest 

asel, s, loc,x, sir, smr 
asel, r, loc,y, spw/2+cww/4 
asel , a, loc, y, 180-spw/2-cww/4 
asel , a, loc,y, 1 . 5*spw+l . 25*cww 
asel, a, loc,y, -180+spw/2+ , 75*cww 
asel,r, loc,x, sir, smr 
cm,pwind,  area 
aatt, 3 , , 1 

! negative  copper  windings  of  interest 

asel , s , loc, X, sir, smr 

asel ,  r ,  loc,y,  -spw/2-*cww/4 

asel, a, loc,y, -180+spw/2+cww/4 

asel , a, loc,y, spw/2+ . 75*cww 

asel, a, loc,y, -180+1. 5*spw+l , 25*cww 

asel , r , loc, x, sir, smr 

cm,  nwind,  area 

aatt, 3 , , 1 

!non  interest  copper  windings 

asel, s, loc, X, sir, smr 

asel, r , loc,y, -180+spw/2+ . 25*cww 

*do, i, l,nsp 

n=-180+spw/2+ . 25 *cww+ (360/nsp) *i 

asel, a, loc,y,n,n 

*enddo 

asel, a, loc,y, -180+spw/2+ . 75*cww 
*do, i, l,nsp 

n=-180+spw/2+.75*cww+ (360/nsp) *i 

asel, a, loc,y,n,n 

*enddo 
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asel, r, loc,x, sir, smr 
cmsel,u,pwind 
cmsel,u, nwind 
cm, owind, area 
aatt , 3 , , 1 

!phase=ac  (3) 

*elseif, phase, eq, 3 

! positive  copper  windings  of  interest 

asel , s, loc, X, sir, smr 
asel, r, loc,y, spw/2+cww/4 
asel , a, loc,y, 180-spw/2-cww/4 
asel, a, loc,y, 180-spw/2- . 75*cww 
asel, a, loc,y, -1. 5*spw-l . 25*cww 
asel , r, loc, x, sir, smr 
cm, pwind, area 
aatt, 3 , , 1 

! negative  copper  windings  of  interest 

asel, s, loc,x, sir, smr 
asel , r , loc,y, -spw/2-cww/4 
asel, a, loc,y, -180+spw/2+cww/4 
asel, a,  loc,y,  180-1 . 5*spw*-l . 25*cww 
asel, a, loc,y, - . 5*spw- . 75*cww 
asel, r, loc,x, sir, smr 
cm,  nwind,  area 
aatt, 3 , , 1 

!non  interest  copper  windings 

asel, s, loc, X, sir, smr 

asel, r, loc,y, -180+spw/2+ . 25*cww 

*do, i, l,nsp 

n=-180+spw/2+ . 25*cww+ {360/nsp) *i 

asel, a, loc,y,n,n 

*enddo 

asel , a, loc,y, -180+spw/2+ . 75*cww 
*do, i, l,nsp 

n=-180+spw/2+ . 75*cww+ (360/nsp) *i 

asel, a, loc,y,n,n 

*enddo 

asel, r , loc, x, sir , smr 
cmsel , u, pwind 
cmsel,u, nwind 
cm,  owind,  area 
aatt, 3 , , 1 

*endif 

! calculate  the  area  of  the  positive  windings  and 
! assign  the  value  to  the  variable  poswind 

cmsel, s, pwind 
asum 

*get, poswind, area, ,area 


!  calculate  the  area  of  the  negative  windings  and 
!  assign  the  value  to  the  variable  negwind 
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cmsel , s , nwind 
asum 

*get , negwind , area , , area 


C  Meshing 

!size  lines  for  meshing 

cmsel , s, spole 
Isla 

Isel, r, loc,x, sir 
lesize,all, , , 6 
Isla 

Isel, r, loc,x, smr 
lesize, all, , , 6 

cmsel , s , owind 
cmsel , a, nwind 
cmsel,a,pwind 
Isla 

Isel , r , loc, X, sir 
lesize, all, , , 6 
Isla 

Isel, r, loc,x, smr 
lesize, all , , , 6 

cmsel, s, stator 
Isla 

Isel , r, loc,x, sor 
lesize, all, , , 3 
Isla 

Isel , r , loc, X, smr 
lesize, all, , , 3 

Isel, s, loc,y, 0 
*do, i , 1, 2*nsp 
n=spw*i 

Isel, a, loc,y, n, n 
*enddo 

Isel , r , loc, X, smr, sor 
lesize, all, , , 3  ,  . 8 

Isel , s , loc,y, -180+spw/2+cww 
*do, i, l,nsp 

n=’-180+spw/2+cww+ (spw+cww)  * 

Isel , a, loc,y,  n,n 

*enddo 

Isel, r, loc,x, sir,  smr 
lesize, all, , , 16,  .6 

Isel, s, loc,y, -180+spw/2 
*do, i , 1 , nsp 

n=-180+spw/2+ (spw+cww) *i 
Isel , a, loc,y, n, n 
*enddo 

lsel,r, loc,x, sir, smr 
lesize, all, , ,16,1.66667 

cmsel , s , owind 
cmsel ,  a,  nwind 
cmsel ,  a,pwind 


Isla 

lesize^all, / , 16,  . 6 

cmsel, s,rpole 
Isla 

Isel , r , loc, X, ror 
lesize, all, , , 10 
Isla 

Isel , r , loc, X, rmr 
lesize,all, , , 10 
Isla 

lesize, all, , , 12 

cmsel, s, inag 
Isla 

Isel , r, loc, X, ror 
lesize, all, , , 10 
Isla 

Isel, r, loc,x, rmr 
lesize, all, , , 2 

cmsel, s, rotor 
Isla 

Isel , r , loc, X, rmr 
lesize, all, , , 6 
Isla 

Isel, r, loc,x, sr 
lesize, all , , , 6 
Isla 

lesize, all , , , 5 , . 8 

cmsel, s, airgap 
Isla 

Isel , r , loc, X, sir 
lesize, all , , . 25 
Isla 

Isel, r, loc,x, ror 
lesize, all , ,0.25 
Isla 

Isel , u, loc,x, sir 
Isel , u, loc,x, ror 
lesize, all,  ,  ,  2 

allsel 
amesh, all 

D.  Assignment  of  Current  Sheets 
! unselect  all  components 
cmsel , none 

! select  the  components  of  the  SRM  that  a  current  density  will  be  applied 
!to  and  then  select  all  elements  of  those  components.  After  this,  the 
!bfe  command  applies  the  current  density  to  the  selected  elements. 

!phase=a  (1) 

* i f , phase , eq , 1 , then 

cmsel , s , nwind 
esla, s 

bfe, all , js , 3 , -current*turns*4/negwind  !-Z  direction 


134 


! +Z  direction 


cmsel,none 
cinsel,s,pwind 
esla, s 

bfe, all, js, 3 , current* turns* 4 /poswind 

!phase=ab  (2) 

*elseif , phase, eq, 2 

cmsel,  s ,  nwind 
esla, s 

bfe, all, js, 3 , -current* turns *4 /negwind  ! -Z  direction 
cmsel , none 
cmsel, s,pwind 
esla, s 

bfe, all, js, 3 , current*turns*4/poswind  ! +Z  direction 

!phase=ac  (3) 

*elseif , phase, eq, 3 

cmsel, s, nwind 
esla, s 

bfe, all , js , 3 , -current*turns*4/negwind  !-Z  direction 

cmsel , none 
cmsel, s,pwind 
esla, s 

bfe, all , js , 3 , current*turns*4/poswind  ! +Z  direction 
*endif 


C.5.  FRSRM 


This  code  differs  in  the  assignment  of  variables,  winding  area  selection  and  assignment  of 
current  sheets.  These  different  sections  of  the  code  are  listed  below. 


A.  Assignment  of  Dimensions  and  Other  Values 

!This  file  calculates  an  electromagnetic  solution  of  the  SRM 
!This  is  the  file  for  the  12/8  SRM  using  fractional  windings 
!The  correct  areas  are  selected  to  get  the  windings  for  the  phase  a,b,c 
!so  that  body-force  current  load  can  be  applied  to  that  area  to  get 
! f rsrm  operation  as  desired. 

!The  areas  are  pwnda,nwnda - >  Phase  a 

!  pwndb,nwndb - >  Phase  b 

!  pwndc,nwndc - >  Phase  c 

'Written  by:  Shyam  S.R, 

!on:  11/6/99 


'variables  of  the  motor  geometry: 


sor=0 . 1786 
smr=0 . 1567 
sir=0.1195 
ror=0.1183 
rmr=0 . 0828 
sr=0.0583 
spw=15 
rpw=17 


!  stator  outside  radius 
!  stator  middle  radius 
!  stator  inside  radius 
!  rotor  outside  radius 
!  rotor  middle  radius 
!  shaft  radius 
!  stator  pole  width 
!  rotor  pole  width 
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theta=~7 . 5 
-  -22.5) 
nsp=12 
nrp=nsp“4 
cww=360/nsp-spw 
iagw=360/nrp-rpw 
mdepth=0 .2366 
currenta=136 . 6174 
currentb=136 . 6174 
currentc=-136 .6174 
turns =12 


! position  of  the  rotor  with  respect  to  the  stator  (0 

! number  of  stator  poles 
! number  of  rotor  poles 
!  copper  winding  width 
!  inner  airgap  width 
! depth  of  the  motor 

! current  applied  to  the  phase  a  -Make  0  for  Quasi 
! current  applied  to  the  phase  b 
! current  applied  to  the  phase  c 
Inumber  of  windings 


B.  Selection  of  Winding  A  reas 

! windings 
!phase=a 

! negative  copper  windings  of  interest 

asel, s, loc,x, sir,  smr 
asel, r, loc,y, spw/2+cww/4 
asel , a, loc,y, 60-spw/2-cww/4 
asel, a, loc,y, 180+spw/2+cww/4 
asel, a, loc,y, 240-spw/2-cww/4 
asel,r, loc,x, sir, smr 
cm,  nwinda ,  area 
aatt, 3 , , 1 

{positive  copper  windings  of  interest 

asel, s, loc,x, sir, smr 

asel , r, loc, y, -30-spw/2-cww/4 

asel, a, loc,y, -90+spw/2+cww/4 

asel, a, loc,y, 150-spw/2-cww/4 

asel,  a,  loc,y,  904*spw/2-fcww/4 

asel, r, loc,x, sir, smr 

cm , pwinda , area 

aatt, 3 , , 1 

!phase=b 

{negative  copper  windings  of  interest 

asel, s, loc,x, sir, smr 
asel, r, loc,y, -spw/2-cww/4 
asel, a, loc,y, 120+spw/2+cww/4 
asel, a, loc,y, 180-spw/2-cww/4 
asel, a, loc, y, -60+spw/2+cww/4 
asel, r, loc,x, sir, smr 
cm,  nwindb,  area 
aatt, 3 , , 1 

{positive  copper  windings  of  interest 

asel, s, loc,x, sir, smr 
asel, r, loc,y, 30+spw/2+cww/4 
asel, a, loc,y, 90-spw/2-cww/4 
asel , a, loc, y, -150+spw/2+cww/4 
asel, a, loc,y, -90-spw/2-cww/4 
asel, r, loc,x, sir, smr 
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cm, pwindb , area 
aatt , 3 , , 1 

!phase=c 

! negative  copper  windings  of  interest 

asel, s, loc, X, sir ,  smr 

asel , r , loc,y, 60+spw/2+cww/4 

asel , a, loc,y, 120-spw/2-cww/4 

asel, a, loc, y, -120+spw/2+cww/4 

asel, a, loc,y, -60-spw/2~cww/4 

asel , r , loc, x, sir,  smr 

cm, nwindc , area 

aatt, 3 , , 1 

! positive  copper  windings  of  interest 

asel, s, loc,x, sir, smr 

asel , r , loc, y, 30-spw/2-cww/4 

asel, a, loc,y, -30+spw/2+cww/4 

asel , a, loc,y, 150+spw/2+cww/4 

asel, a, loc,y, -150-spw/2-cww/4 

asel, r, loc,x, sir, smr 

cm ,  pwindc ,  area 

aatt, 3 , , 1 


'calculate  the  area  of  the  positive  windings  and 
! assign  the  value  to  the  variable  poswind 

cmsel , s , pwinda 
asum 

*get,poswinda, area, ,area 


! calculate  the  area  of  the  negative  windings  and 
! assign  the  value  to  the  variable  negwind 

cmsel ,  s ,  nwinda 
asum 

*get , negwinda , area , , area 

! calculate  the  area  of  the  positive  windings  and 
! assign  the  value  to  the  variable  poswind 

cmsel , s , pwindb 
asum 

*get,poswindb, area, ,area 


! calculate  the  area  of  the  negative  windings  and 
{assign  the  value  to  the  variable  negwind 

cmsel , s , nwindb 
asum 

*get , negwindb, area, ,area 

{calculate  the  area  of  the  positive  windings  and 
{assign  the  value  to  the  variable  poswind 

cmsel , s , pwindc 
asum 
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get , poswindc , area , , area 


! calculate  the  area  of  the  negative  windings  and 
! assign  the  value  to  the  variable  negwind 

cmsel, s,nwindc 
asum 

*get , negwindc , area , , area 


C.  Assignment  of  Current  Sheets 
1 unselect  all  components 
cmsel,none 

! select  the  components  of  the  SRM  that  a  current  density  will  be 
applied 

!to  and  then  select  all  elements  of  those  components.  After  this,  the 
!bfe  command  applies  the  current  density  to  the  selected  elements. 

!phase=a 

cmsel , s , nwinda 
esla, s 

bfe, all, js, 3 , “Currenta*turns*4/negwinda  ! -Z  direction 
cmsel , none 
cmsel, s,pwinda 
esla, s 

bfe, all, js, 3 , currenta*turns*4/poswinda  ! +Z  direction 

!phase=b 

cmsel , s , nwindb 
esla, s 

bfe, all,  js, 3 , -currentb*turns*4/negwindb  ! -Z  direction 
cmsel , none 
cmsel , s , pwindb 
esla, s 

bfe, all , js, 3 , currentb*turns*4/poswindb  ! +Z  direction 

!phase=c 

cmsel , s , nwindc 
esla, s 

bfe, all, js, 3 , -currentc*turns *4 /negwindc  ! -Z  direction 
cmsel , none 
cmsel, s,pwindc 
esla, s 

bfe, all, js, 3 , currentc*turns*4/poswindc  ! +Z  direction 


C.6.  Full-Pitch  SRM 
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Again,  this  code  differs  in  the  assignment  of  variables,  winding  area  selection  and  assignment  of 
current  sheets.  These  different  sections  of  the  code  are  listed  below. 


A.  Assignment  of  Dimensions  and  Other  Values 


•This  file  calculates  an  electromagnetic  solution  of  the  SRM 
IThis  is  the  file  for  the  12/8  SRM  using  fractional  windings 
!The  correct  areas  are  selected  to  get  the  windings  for  the  phase  a,b,c 
Iso  that  body-force  current  load  can  be  applied  to  that  area  to  get 
!  f rsrm  operation  as  desired. 

IThe  areas  are  pwnda,nwnda - >  Phase  a 

!  pwndb,nwndb - >  Phase  b 

!  pwndc,nwndc - >  Phase  c 

•Written  by:  Shyam  S.R. 

•on:  11/6/99 


•variables  of  the  motor  geometry: 


sor=0.1721 

smr=0 . 1502 

sir=0.1195 

ror=0.1183 

rmr=0 . 0828 

sr=0 . 0583 

spw=15 

rpw=17 

theta=-15 

-  -22.5) 

nsp=12 

nrp=nsp-4 

cww=360/nsp-spw 

i agw=  3  6  0 / nr p - rpw 

mdepth=0 .2366 

current=136 . 6174 

turns=12 


!  stator  outside  radius 
!  stator  middle  radius 
!  stator  inside  radius 
!  rotor  outside  radius 
•rotor  middle  radius 
!  shaft  radius 
!  stator  pole  width 
!  rotor  pole  width 

! position  of  the  rotor  with  respect  to  the  stator  (0 

•number  of  stator  poles 
Inumber  of  rotor  poles 
I  copper  winding  width 
•  inner  airgap  width 
I  depth  of  the  motor 
I  current  applied  to  the  phase  a 
Inumber  of  turns  per  pole 


B.  Selection  of  Winding  A  reas 

! windings 

! positive  copper  windings  of  interest 

asel, s, loc, X, sir, smr 
asel, r, loc,y, spw/2+cww/4 
asel , a, loc, y, 30-spw/2-cww/4 
asel , a, loc, y, 90-spw/2-cww/4 
asel , a, loc, y, 60+spw/2+cww/4 
asel , a, loc,y, 180+spw/2+cww/4 
asel , a, loc,y, 210-spw/2-cww/4 
asel , a, loc,y, 270-spw/2-cww/4 
asel , a, loc,y, 240+spw/2+cww/4 
asel, r, loc,x, sir, smr 
cm , pwind , area 
aatt, 3 , , 1 

! negative  copper  windings  of  interest 

asel , s, loc, X, sir, smr 
asel, r, loc,y, -spw/2-cww/4 
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asel,a,  loc,y,  -30+spw/2+cww/4 
asel , a, loc,y, 90+spw/2+cww/4 
asel,a, loc,y, 120-spw/2-cww/4 
asel, a, loc,y, 180-spw/2-cww/4 
asel , a, loc, y, 150+spw/2+cww/4 
asel , a, loc,y, -90+spw/2+cww/4 
asel,  a,  loc,y,  -60-spw/2“"Cww/4 
asel, r, loc,x, sir, smr 
cm, nwind, area 
aatt , 3 , , 1 

I  non  interest  copper  windings 

asel, s, loc,x, sir, smr 

asel,r, loc,y, -180+spw/2+ . 25*cww 

*do, i , 1, nsp 

n=-180+spw/2+ . 25*cww+ (3 60 /nsp) *i 

asel , a, loc,y, n, n 

*enddo 

asel , a, loc, y, -180+spw/2+ . 75*cww 
*do, i, l,nsp 

n=-180+spw/2+.75*cww+ (360/nsp) *i 

asel, a, loc,y,n,n 

*enddo 

asel, r, loc,x, sir,  smr 
cmsel,  u,pwind 
cmsel,u, nwind 
cm, owind, area 
aatt, 3 , , 1 


! calculate  the  area  of  the  positive  windings  and 
! assign  the  value  to  the  variable  poswind 

cmsel , s , pwind 
asum 

*get, poswind, area, ,area 

! calculate  the  area  of  the  negative  windings  and 
! assign  the  value  to  the  variable  negwind 

cmsel , s , nwind 
asum 

*get , negwind , area , , area 


C  Assignment  of  Current  Sheets 
! unselect  all  components 
cmsel , none 

’select  the  components  of  the  SRM  that  a  current  density  will  be 
applied 

!to  and  then  select  all  elements  of  those  components.  After  this,  the 
!bfe  command  applies  the  current  density  to  the  selected  elements. 

cmsel, s, nwind 
esla,  s 

bfe, all, js, 3 , -current* turns *8 /negwind  ! -Z  direction 
cmsel , none 
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cmsel, s,pwind 
esla,  s 

bfe, all , js, 3 , current* turns* 8 /poswind 


l+Z  direction 


C.7.  Writing  Code  for  Automated  Looping 

Using  ANSYS™,  it  is  possible  to  set  up  the  code  to  perform  automated  looping  with 
variation  of  file  variables  in  each  run.  The  code  should  be  re-written  as  a  .mac  file  and  is  run 
under  the  control  of  another  file  that  invokes  the  .mac  file.  The  variables  that  must  be  varied  are 
re-named  as  argl,  arg2, . .  argl8. 

For  example,  the  6/4  ANSYS  code  file  is  re-written  as  a  .mac  file  srminduc.mac  to  determine 
the  flux-linkage  at  different  rotor  positions  and  phase  current  values.  At  the  start  of  the  .mac  file, 
the  “/clear,nostart”  operation  is  performed.  A  file  is  opened  in  the  controlling  file  using  the 
command  *CFOPEN  and  the  calculated  values  during  each  run  are  written  into  it  using  the 
*VWRrrE  command.  In  the  .mac  file,  argl  is  the  rotor  position  and  arg2  and  arg3  are  the 
currents  in  the  working  and  leading  phases,  respectively. 


Controlling  File  srmc 

*CFOPEN, linkage 8 6, out, /hoine/ssr2/ 
arg3=0 


*DO,argl, -30,0,2 
*DO,arg2, 13 .5,135,13.5 

srminduc , argl , arg2 , arg3  ! srminduc .mac  invoked 

*ENDDO 

*ENDDO 

!This  file  enables  calculation  of  the  solution  for  8/6  SRM  characteristics 
!to  determine  the  coenergy  and  flux-linkages  which  can  be  used  to  get  a 
! multi -phase  model  of  the  SRM 

!At  the  end  of  calculation  for  each  argl  value,  *VWRITE  command  is  used 
! to  write  argl , arg2 , arg3 , coenergy  and  torque  values  to  file  linkagesfp.out 


IWritten  by  Shyam  S.R.  on  11/11/99 


Skeleton  of  .mac  file  srminduc.mac 
/clear, nostart 
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!This  file  calculates  an  electromagnetic  solution  of  the  SRM 

! Written  by:  Shyam  S.R. 

!on:  11/11/99 


! variables  of  the  motor  geometry: 


sor=0.1851 
smr=0 . 1524 
sir=0.1195 
ror=0.1183 
rmr=0.0828 
sr  =0.0477 
spw=22 . 5 
rpw=24 . 5 
nsp=8 
nrp=nsp-2 
cww=360/nsp-spw 
iagw=3 60 /nrp-rpw 
mdepth=0 .2366 
turns=25 


!  stator  outside  radius 
!  stator  middle  radius 
!  stator  inside  radius 
! rotor  outside  radius 
!  rotor  middle  radius 
!  shaft  radius 
!  stator  pole  width 
!  rotor  pole  width 
Inumber  of  stator  poles 
Inumber  of  rotor  poles 

I  copper  winding  width 
I  inner  airgap  width 
I  depth  of  the  motor 
'number  of  windings 


I  calculate  the  angle  from  the  origin  to  the  top  of  the 
Ismail  angle  on  the  pole  of  the  rotor 


I  numerator  =  num 

num=tan( 5*3. 14159265/ 180) * (ror-rmr) +ror*sin (rpw/2 *3 . 14159265/180) 
phi=asin (num/rmr) *180/3.14159265 

I  enter  the  preprocessor  of  the  program 

/prep7 

'define  the  element  types  to  be  used 
et, l,planel3, 0 


Ikeypoints  of  the  stator  poles  and  copper  windings 

k, , sir , spw/2 

k, , smr, spw/2 

k, , sir, spw/2+cww/2 

k,  ,  smr,  spw/2+cww/2 

k, , sir , spw/2 +CWW 

k, , smr, spw/2+cww 

ksel, s, ,,1,6 

kgen, nsp, all, , , , 360/nsp 

Ikeypoints  of  the  rotor  and  shaft 

k, , rmr , argl-phi 
k, ,ror, argl -rpw/2 
k, , ror, argl+rpw/2 
k, ,  rmr, argl+phi 
ksel,s, , , nsp*6+l , nsp*6+4 
kgen, nrp, all, , , , 360/nrp 

Hines  for  the  rotor  poles 


142 


csys , 0 

1, 6*nsp+l, 6*nsp+2 
1 , 6*nsp+3 , 6*nsp+4 
1, 6*nsp+5, 6*nsp+6 
1, 6*nsp+7, 6*nsp+8 
1, 6*nsp+9, 6*nsp+10 
1, 6*nsp+ll, 6*nsp+12 
1, 6*nsp+13, 6*nsp+14 
1, 6*nsp+15; 6*nsp+16 
1, 6*nsp+17, 6*nsp+18 
1, 6*nsp+19, 6*nsp+20 
1, 6*nsp+21, 6*nsp+22 
1 , 6*nsp+23 , 6*nsp+24 
csys , 1 

! select  all  keypoints  of  the  SRM 

ksel, all 
kplot, all 

! generate  areas  of  the  SRM 

! areas  of  the  stator  ring 

pcirc, smr, sor, 0, spw 
agen,2*nsp,all, , , ,360/ {2*nsp) 


cmsel,none 

! select  the  components  of  the  SRM  that  a  current  density  will  be  applied 
!to  and  then  select  all  elements  of  those  components.  After  this,  the 
!bfe  command  applies  the  current  density  to  the  selected  elements. 

cmsel,  s,nwind2 
esla, s 

bfe, all, js, 3 , -arg2*turns*2/negwind  !-Z  direction 
cmsel, none 
cmsel , s ,pwind2 
esla, s 

bfe, all, js, 3 ,arg2*turns*2/poswind  !+Z  direction 

cmsel, s,nwind3 
esla, s 

bfe,all, js,3,-arg3*turns*2/negwind  ! -Z  direction 
cmsel, none 
cmsel , s , pwind3 
esla, s 

bfe,all, js,3,arg3*turns*2/poswind  !+Z  direction 

!far-field  boundary  conditions 
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/nopr 

parsav,  all ,  sasiitiac 
esel, s,mat, , 1,4;  1 
etable, sene, sene 
/out, scratch 
ssum 
/out 

*get, energy, ssum, , item, sene 
/gopr 

menergy=energy*mdepth 

esel, s,mat, ,1,4,1 
s energy, 1 , 0 
c  o  en = c_eng  *mdep  th 

esel, s,mat, ,1,4,1 
s energy, 0 , 0 
ener=s_eng*mdepth 

! calculate  the  airgap  torque 

allsel 

cmsel, s, airgap 
esla 

torqc2d,  ror+  {  (sir-ror )  /2 )  ,  100 
agt=torque*mdepth 

*VWRITE , argl , arg2 , arg3 , agt , coen , ener 
( F6 . 1 , F6 . 1 , F8 . 1 , F8 . 2 , F9 . 5 , F9 . 5 ) 


finish 
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APPENDIX  D 


SIMULINK  MODELS 

D.l  SIMULINK  Model  for  the  SRM  using  ANN 

This  Section  describes  the  changes  to  the  SIMULINK™  implementation  of  the  SRM  single¬ 
phase  model  in  [1]  when  using  the  ANN  to  account  for  the  multi-phase  excitation.  The  ANN 
suggested  in  this  Report  map  the  rotor  position  with  respect  to  a  phase,  the  current  in  the  same 
phase  (called  the  working  phase)  and  the  one  ahead  of  it  in  the  direction  of  rotation  (called  the 
leading  phase)  to  the  electromagnetic  torque  and  mutual  flux  interaction  functions.  The  overall 
system  and  the  SIMULINK  blocks  for  torque  and  mutual  flux  computations  are  shown  in  Figures 
D.  1,  D.2  and  D.3  for  the  simulation  of  a  four-phase  8/6  SRM. 


R_P 

Figure  D.  1  SIMULINK™  model  of  the  four-phase  8/6  SRM. 
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Figure  B.3  The  ANN-based  mutual  flux  computation  block. 
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In  the  development  of  this  simulation  approach,  the  assumption  is  that  more  than  two  phases 
do  not  operate  simultaneously.  The  torque  computation  block  uses  selective  switching  at  the 
output  to  select  the  phase  which  is  currently  the  working  phase  and  derive  the  torque  output  with 
respect  to  that  phase.  In  Figure  D.3,  FMl,  FM2,  FM3  and  FM4  denote  respectively  the  mutual 
flux  component  of  the  net  flux  for  each  phase.  Hence,  this  component  is  used  to  correct  the  phase 
flux  for  the  mutual  interactions. 

D.2  SIMULINK™  Model  for  On-Line  Torque  Estimation 

The  simulation  of  the  on-line  torque  estimation  system  is  based  on  measurements  of  SRM 
SIMULINK  model  variables.  The  phase  current  ik,  phase  rotor  position  thk  and  the  terminal 
voltage  vk  are  measured  quantities.  The  SRM  model  of  [1]  has  one  time  step  of  delay  between  the 
terminal  voltage  and  actual  phase  current.  To  account  for  this,  a  delayed  value  of  the  terminal 
voltage  is  considered  in  the  estimation  system.  The  shift  register  is  used  to  store  sampled  values 
over  several  time  steps  as  determined  by  a  buffer  size  specified  in  the  initialization  stage.  Several 
values  of  the  buffer  size  were  tried  out  to  see  the  performance  of  the  estimation  system  when  it  is 
trained  to  greater  accuracy  near  the  current  operating  point  of  the  SRM.  However,  it  was  found 
that  a  buffer  size  of  1  results  in  the  best  performance  since  the  ANN  is  able  to  train  to  represent 
the  system  completely.  Figure  D.4  shows  the  portion  of  the  model  that  is  used  for  on-line 
estimation.  The  MATLAB  function  used  for  the  ANN  and  its  inversion  is  also  included. 


MATLAB  Function  for  ANN  Training  and  Inversion  in  the  Adaptive  Torque  Estimation 
Simulation: 

function  [out] =randcots2ialt ( inik, inikl, in thk, inthkl, invl) 
global  net  netl  bufsize  vw  n\iinln 
T=  5.0e-5  ;  %the  sampling  time 

Inp=4;  %number  of  inputs 

totB=numIn+l ; 
for  i=l:bufsize 
P(l, i) =inik(i) ; 

P{2, i) =invl (i) ; 

P{3, i) =inthk(i) ; 

P(4, i) =inthk{i) -inthkl (i) ; 
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Figure  D.4  On-line  torque  estimation  system 


Tr (i) =inik(i) -inikl (i) ; 

end 

net=train (net , P, Tr) ; 

B=cell2niat  (net  .b)  ; 

test (1) =inik (buf size) ; 

test (2 ) =invl (buf size) ; 

test (3) =inthk(bufsize) ; 

test (4) =inthk(bufsize) -inthkl (buf size) ; 

%Detennining  the  network  output  for  the  last  case 
iw=sini (net, test ') ;  %  Rate  of  current  rise 

fl=200000; 
if  vw>5 

min=vW”5 ; 
inax=vw+5 ; 
elseif  vw  <  1 
min=l ; 
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max=6; 

else 

min=vw 
max=vw+5 ; 
end 

for  i=min: 0 . 25 :max 

test ( 1 ) =inik (buf size) ; 
test(2)=i-l; 
test (3 ) =inthk (bufsize) ; 
test (4) =inthk(bufsize) -inthkl (bufsize) ; 
f=sim(net, test ' ) ; 
if  abs(f)<fl 
vw=i-l ; 
fl^abs(f) ; 

end 

end 

Tr=vw*inik (bufsize) *180*T/ (inthk (bufsize) -inthkl (bufsize) ) /pi; 
out= [iw, vw, Tr] ; 
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